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Preface 


The  purpose  of  this  study  was  to  develop  a  computer  program  that 
could  model  thin  film  growth  by  vapor  deposition.  Although  other  com¬ 
puter  programs  have  been  written  to  do  the  same,  this  one  is  unique  in 
that  it  incorporates  several  deposition  parameters  both  new  and  old. 

Some  of  these  include  the  variation  of  the  angle  of  incidence  about  some 
main  angle,  the  mobility  of  molecules,  the  doping  of  impurities  into  the 
microstructure,  the  creation  of  multilayered  films  and  the  ability  to 
vary  the  substrate.  Like  many  other  models  this  is  a  two-dimensional 
simulation. 

The  program  is  written  in  FORTRAN  77  and  was  designed  to  be  used  on 
a  VAX  computer  under  a  VMS  operating  system. 

Throughout  the  program  development,  a  process  which  took  over  five 
months,  it  was  neccessary  to  rewrite  sections  of  the  program  in  order  to 
get  rid  of  bugs  and/or  make  improvements.  Although  the  resulting  pro¬ 
gram  is  fairly  sophisticated  in  nature,  it  is  by  no  means  the  ultimate 
in  the  modeling  of  vapor  deposition.  More  programs  should  be  written 
in  the  future  to  provide  a  better  understanding  of  thin  films  and  the 
effect  vapor  deposition  has  on  their  growth. 

While  producing  the  computer  model  and  writing  this  thesis  I  was 
very  fortunate  to  be  surrounded  by  several  kind  and  dedicated  individ¬ 
uals.  I  am  deeply  indebted  to  my  faculty  advisor,  Major  John  Wharton, 
for  having  faith  in  me  and  for  his  assistance  in  times  of  need.  I  also 
wish  to  thank  Dr.  David  Lee  of  the  math  department  for  his  assistance  in 
the  early  stages  of  program  development.  A  word  of  thanks  is  also  owed 
to  the  system  manager  Frank  Bakos  for  his  help  in  answering  my  many 
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Abstract 

A  VAX-11/785  computer  was  used  to  simulate  the  two-dimensional 
growth  of  thin  films  produced  by  vapor  deposition.  In  this  model  mole¬ 
cules  and  impurities  were  represented  by  three  different  sized  disks. 

In  order  to  simulate  varying  deposition  conditions  and  evaporants,  sev¬ 
eral  variable  parameters  were  introduced.  Among  these  parameters  were 
the  variation  of  the  deposition  angle  about  some  main  angle,  the  mobil¬ 
ity  of  the  disks  upon  collision,  the  ability  to  introduce  impurities 
into  the  microstructure,  the  simulation  of  multilayered  coatings  and 
the  ability  to  introduce  imperfections  into  the  substrate. 

The  results  obtained  by  this  model  show  that  disks  can  be  used  to 
simulate  some  of  the  main  features  exhibited  by  vapor  deposited  films. 
Among  these  features  are  the  formation  of  columns  and  their  compliance 
with  the  "tangent  rule",  and  the  disappearance  of  this  structure  in  the 
case  of  large  disk  mobility.  Another  feature  found  to  be  exhibited  in 
the  modeled  films  is  that  under  certian  conditions,  impurities  and  sub¬ 
strate  imperfections  can  produce  large  voids  and/or  nodules.  Other 
characteristics  found  in  the  simulated  films  include  pores  which  could 
allow  water  absorption,  and  increased  packing  density  for  films  produced 
with  angle  variations  along  with  a  moderate  amount  of  disk  mobility. 


AN  IMPROVED  MODEL  OF  THIN  FILM  GROWTH 


I.  Introduction 

Thin  films  are  widely  used  for  various  scientific  and  industrial  ap¬ 
plications,  most  notably  in  microelectronics  and  optics.  In  microelec¬ 
tronics,  their  use  is  predominately  in  the  manufacturing  of  integrated 
circuits  as  in  the  isolation  of  circuit  elements.  In  optics,  however, 
thin  films  have  been  used  for  a  variety  of  purposes  such  as  antireflec¬ 
tion  coatings,  mirrors,  sun  protection  coatings  as  well  as  band-pass  fil¬ 
ters.  The  fabrication  technique  used  in  producing  thin  films  varies  but 
the  most  commonly  used  technique  in  optics  is  vapor  deposition. 

Since  the  advent  of  using  vacuum  chambers  for  the  vapor  deposition 
of  thin  films,  the  process  by  which  the  mircostructure  is  formed  has  not 
been  completely  understood.  The  amorphous  and  sometimes  crystalline 
thin  films  produced  by  this  type  of  deposition  are  usually  characterized 
by  voids,  pores  and  other  vacancies.  However,  the  most  dominant  feature 
of  these  films  is  their  columnar  shape.  More  specifically,  the  micro¬ 
structure  consist  of  an  array  of  somewhat  parallel  columnar  regions  of 
high  density  material  surrounded  by  a  network  of  material  of  lower  den¬ 
sity.  It  is  also  a  characteristic  of  this  microstructure  that  the  dia¬ 
meters  of  the  individual  columns  are  approximately  uniform.  This  struc¬ 
ture  has  been  observed  by  microfractography  (13,18-21,29,40,41),  by 
transmission  electron  microscopy  (6,29,36-38,43,48),  as  well  as  by  both 
small  angle  electron  (29,38,46),  and  x-ray  scattering  (4,34,45).  Exper- 


iments  have  shown  that  this  columnar  structure  can  also  change  some  of 
the  properties  of  the  material  when  compared  to  those  of  the  bulk  mater¬ 
ial.  Some  of  the  changes  that  have  been  noted  include  magnetic  (19-21, 
29,33,40,41,43),  optical  (7,8),  electrical  (8,22,42)  and  mechanical  (41) 
properties. 

In  an  effort  to  try  to  understand  the  growth  patterns  of  the  colum¬ 
nar  structure  and  how  this  might  affect  the  properties  of  the  film,  the 
deposition  process  was  first  simulated  by  use  of  nucleation  models  (1,2, 
9,25,42).  Since  then  researchers  have  used  computer  models  to  simulate 
the  vertical  structure  of  thin  films  (3,5,10,11,15,16,24,26,28,30,31,46, 
47,49).  The  results  of  these  latest  models  confirm  the  fact  that  shadow¬ 
ing  plays  a  major  role  in  the  growth  of  films  deposited  at  oblique  inci¬ 
dence.  However,  the  columns  comprising  the  mircrostructure  produced  by 
the  models  are  approximately  one  to  two  orders  of  magnitude  too  small  in 
diameter.  Other  characteristics  that  past  models  fail  to  properly  simu¬ 
late  include  the  tendency  of  the  models  to  produce  films  with  low  pack¬ 
ing  densities,  and  a  failure  to  produce  films  with  a  crystalline  struc¬ 
ture. 


Objective 

The  objective  of  this  research  was  to  improve  upon  a  previous  two- 
dimensional  computer  model  (47)  of  thin  film  vapor  deposition;  and  in 
doing  so,  incorporating  into  the  new  model  parameters  which  could  be 
changed  at  will  by  the  user  of  the  program.  This  allows  the  user  to 
then  examine  the  simulated  film  and  make  observations  concerning  the 
effect  a  particular  parameter  has  made.  The  first  parameter  to  be  in¬ 
troduced  into  the  model  was  the  variation  of  the  angle  of  deposition 
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about  some  main  angle.  Other  parameters  introduced  into  the  new  model 
include  allowing  a  certain  percentage  of  molecules  to  undergo  more  than 
normal  mobility,  the  ability  to  simulate  the  doping  of  impurities  into 
the  microstructure  as  well  as  being  able  to  modify  the  substrate  to  mo¬ 
del  imperfections.  The  final  parameter  to  be  put  into  the  program  was 
the  ability  to  deposit  two  different  sized  molecules  and  in  doing  so 
simulate  multilayered  coatings. 

General  Approach 

The  approach  used  to  accomplish  the  objective  was  to  divide  the 
problem  into  three  parts.  First,  after  the  problem  had  been  defined  and 
goals  set,  a  comprehensive  literature  search  was  made  from  available  ma¬ 
terial.  In  particular,  any  papers  found  containing  computer  simulations 
of  vapor  deposition  were  especially  scrutinized.  The  second  step  in  sol¬ 
ving  the  problem  consisted  of  developing  the  computer  program  used  to  mo¬ 
del  the  film  growth.  This  included  making  any  simplifying  assumptions, 
developing  the  mathematics  and  the  model,  as  well  as  testing  the  pro¬ 
gram.  The  last  step  consisted  of  running  the  program  and  then  analyz¬ 
ing  the  resulting  films.  Throughout  this  last  step  the  previously  men¬ 
tioned  parameters  were  varied  and  any  significant  changes  to  the  film 
noted . 

Sequence  of  Presentation 

The  results  of  this  study  will  be  presented  in  the  five  remaining 
chapters  and  three  appendices.  In  chapter  two,  some  of  the  current  know¬ 
ledge  pertaining  to  the  new  model  on  thin  films  will  be  reviewed.  In 
chapter  three,  the  basic  requirements  needed  to  produce  a  computer  model 


and  how  these  are  incorporated  into  the  new  model  will  be  discussed. 
Chapter  four  will  present  the  new  model  and  how  it  was  derived  including 
its  organization  and  implementation.  The  fifth  chapter  will  show  some 
of  the  resulting  films  produced  by  the  computer  simulation  as  well  as 
the  analysis  of  these  films.  A  summary  of  conclusions  and  recommenda¬ 
tions  will  be  presented  in  the  sixth  chapter. 

The  final  section  of  this  thesis  contains  three  appendices  which 
are  provided  for  follow  on  studies.  Appendix  A  contains  several  mathe¬ 
matical  derivations  used  in  the  development  of  the  program,  while  Ap¬ 
pendix  B  contains  the  actual  Fortran  code.  The  last  appendix,  which  is 
extremely  important  for  follow  on  studies,  contains  the  documentation  of 
the  computer  program,  a  brief  explanation  of  how  the  program  is  organ¬ 
ized,  and  a  list  of  the  variables  used  along  with  a  description  of  their 


purpose 


Although  there  are  several  ways  of  depositing  thin  films  onto  a  sub¬ 
strate,  the  most  commonly  used  technique  for  producing  optical  compon¬ 
ents  is  by  physical  vapor  deposition.  The  reason  for  using  this  tech¬ 
nique  is  that  it  tends  to  produce  a  more  uniform  coating  of  the  sub¬ 
strates.  In  most  cases  this  process  is  accomplished  inside  a  vacuum 
chamber  where  the  bulk  material  is  either  heated  resisti'vely  or  by  an 
electron  beam  which  is  directed  into  the  material  (23).  Upon  heating, 
the  vapor  thus  produced  coats  the  substrates  placed  inside  the  chamber 
by  condensation. 


Columnar  Structure 

As  mentioned  previously  in  chapter  one,  the  most  dominant  feature 
exhibited  by  vapor  deposited  films  is  their  columnar  structure.  It  is 
important  to  realize,  however,  that  not  all  vapor  deposited  films  ex¬ 
hibit  this  columnar  structure  due  to  the  mobility  of  the  molecules. 

This  process  will  be  discussed  in  a  future  section. 

Tangent  Rule.  Of  all  the  parameters  which  affect  the  columnar 
structure  exhibited  by  some  films,  the  most  influential  is  the  angle  of 
deposition  (36,37,39,44).  In  1966  Nieuwenhuizen  and  Haanstra  (39)  in¬ 
troduced  an  empirically  derived  expression  which  has  become  known  as  the 
"tangent  rule."  This  expression  relates  the  angle  of  growth  of  the  film 
to  the  angle  of  deposition  in  the  following  manner: 

2  tan  B  -  tan  A  (1) 

where  B  is  the  angle  of  growth  of  the  film  as  measured  from  the  sub- 
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strate  normal,  and  A  is  the  angle  of  deposition  also  measured  relative 
to  the  substrate  normal.  Figure  1  shows  a  graph  of  the  tangent  rule 
for  deposition  angles  from  0  to  90  degrees  (49).  It  is  true  that  not 
all  films  exhibiting  columnar  structure  follow  the  rule  exactly,  but  for 
the  most  part  the  tangent  rule  is  still  obeyed. 

Shadowing.  One  of  the  main  reasons  that  vapor  deposited  films  pro¬ 
duce  a  columnar  structure  and  that  most  films  follow  the  tangent  rule 
is  due  to  a  phenomenon  called  "shadowing."  Put  simply,  shadowing  is  an 
effect  which  limits  the  eventual  placement  of  the  incident  molecules  in 
the  microstructure  due  to  previously  deposited  molecules  blocking  their 
path.  Under  actual  deposition  conditions,  a  molecule’s  trajectory  is  a 
straight  line  due  to  the  fact  that  the  process  takes  place  in  a  vacuum. 
In  other  words,  the  molecule's  mean  free  path  is  very  large  implying 
that  the  molecules  first  collision  will  be  with  the  substrate  or  pre¬ 
viously  deposited  film.  If  one  allows  disks  to  represent  molecules,  the 
idea  of  shadowing  can  be  understood  in  figure  2.  If  disk  X  represents 
a  previously  deposited  molecule,  it  can  be  thought  of  as  a  barrier  which 
other  molecules  may  encounter.  For  instance,  if  disk  Y  had  a  trajectory 
to  the  left  of  disk  X  at  an  angle  A,  the  closest  that  it  could  impact 
the  substrate  would  be  some  distance  away.  In  other  words,  previously 
deposited  molecules  or  even  imperfections  in  the  substrate  shield  or  sha¬ 
dow  unoccupied  sites,  thus  creating  voids.  If  the  mobility  of  the  mole¬ 
cules  after  condensation  do  not  fill  up  the  voids,  the  void  structure  is 
maintained  and  grows  with  subsequent  deposition  of  molecules.  As  im¬ 
plied  in  figure  2,  the  larger  the  angle  of  deposition  the  larger  the  re- 


Shape  of  the  Columnar  Structure.  The  angle  of  deposition  also  has 


an  effect  on  the  shape  of  the  individual  columns.  For  normal  deposition 
the  cross  sections  of  the  columns  tend  to  be  circular.  However,  films 
deposited  at  oblique  incidence  yield  columns  with  cross  sections  that 
are  elliptical  in  shape.  Contrary  to  what  one  might  think,  the  axis  of 
elongation  occurs  in  a  plane  perpendicular  to  the  plane  of  incidence 
(30),  and  is  a  direct  result  of  shadowing.  A  figure  exhibiting  this 
elongation  can  be  seen  in  figure  3.  As  the  angle  of  deposition  in¬ 
creases  this  columnar  shape  becomes  more  elliptical,  and  the  void  re¬ 
gions  between  the  columns  become  thinner  along  a  direction  parallel  to 
the  incident  plane. 

Another  characteristic  of  columnar  shape  that  is  gaining  more  sup¬ 
port  is  that  the  individual  columns  are  made  up  of  even  smaller  columns, 
usually  referred  to  as  dendrites  (10,31).  These  dendrites,  which  may 
only  be  a  few  molecules  in  diameter,  are  thought  to  come  together  in  the 
early  stages  of  growth  and  intertwine,  thus  forming  the  relatively  large 
columns  seen  by  transmission  electron  microscopy. 


Density.  Because  the  amount  of  voids  is  related  to  the  angle  of  in 
cidence,  the  resulting  density  of  the  film  is  also  affected.  As  expec¬ 
ted,  the  density  of  the  films  decreases  with  an  increasing  angle  of  depo' 
sition.  The  density  of  a  film  can  be  defined  in  several  ways,  but  the 
most  commonly  referred  to  density  is  the  packing  density.  Packing  den¬ 
sity  is  usually  defined  as  follows: 

Volume  of  the  solid  part  of  the  film  (2) 

^  *  Total  volume  of  the  film  (solid  +  voids) 

where  p  is  the  packing  density. 

Angle  Variation.  One  of  the  many  ways  that  people  have  gone  about 
affecting  the  packing  density  and  producing  more  uniform  coatings  is  by 
rotating  the  substrates.  By  rotating  the  substrates  the  angle  of  depo¬ 
sition  is  varied  over  time,  thereby  smoothing  out  the  film  and  making 
it  more  uniform.  This  also  tends  to  increase  the  overall  density  of 
the  film.  As  to  date,  there  have  been  few  if  any  studies  looking  at  the 
amount  of  variation  as  well  as  the  speed  of  rotation,  and  how  this  af¬ 
fects  the  structure  of  the  film. 

Mobility 

In  order  to  further  affect  the  density  of  the  films,  that  is  trying 
to  produce  more  dense  films,  many  people  have  tried  to  increase  the  mo¬ 
bility  of  the  molecules  as  they  are  being  deposited.  There  are  basi¬ 
cally  two  reasons  why  one  would  want  to  affect  the  mobility  of  the  inci¬ 
dent  molecules  as  they  are  being  deposited.  One,  as  just  mentioned, 


would  be  to  try  to  produce  a  film  with  a  higher  packing  density.  The 
reason  this  is  desirable  is  that  if  the  density  of  the  film  can  approach 


that  of  the  bulk  material,  the  properties  of  the  film  will  tend  to  be 
more  stable,  especially  the  optical  ones.  The  second  reason  for  wanting 
to  affect  the  mobility  of  the  molecules  is  to  help  cut  down  on  the 
amount  of  water  absorption.  As  more  water  is  absorbed  into  the  mircro- 
structure  the  properties  of  the  film  change,  and  in  the  case  of  optical 
properties  wavelength  shifts  can  occur  (31). 

Substrate  Temperature.  It  is  well  known  that  when  the  ratio  of  the 
substrate  temperature  to  that  of  the  bulk  material  (the  material  being 
deposited  as  a  film)  is  extremely  high,  the  resulting  mircrostructure  is 
without  columnar  structure  (16).  This  is  due  to  the  fact  that  the  mobil¬ 
ity  of  the  molecules  is  large  and  most,  if  not  all,  of  the  voids  are 
filled.  The  high  substrate  temperature  also  has  the  effect  that  it  in¬ 
creases  the  density  of  the  film.  The  problem  with  using  substrate  tem¬ 
perature  to  increase  the  density  of  films  is  that  most  films  are  com¬ 
posed  of  several  layers  and  the  high  temperature  would  destroy  the  lay¬ 
ered  structure. 

Ion  Bombardment.  The  use  of  ions  to  increase  the  mobility  of  the 
molecules  making  up  the  film  has  two  primary  purposes.  First,  ion  bom¬ 
bardment  has  been  used  to  help  anneal  the  film  to  the  substrate.  Se¬ 
cond,  the  use  of  bombarding  ions  has  been  used  to  help  increase  the  over¬ 
all  packing  density  of  the  film.  In  both  cases,  authors  (32)  have  exper¬ 
imented  with  varying  current  densities  in  hopes  of  producing  films  with 
greater  substrate  adhesion  and  packing  densities. 

Ultra-Violet  Radiation.  Another  technique  used  to  increase  the  mo¬ 
bility  of  the  molecules  under  deposition  has  been  the  use  of  UV  radia¬ 
tion  (14,49).  Although  this  technique  has  not  had  as  much  success  as 


ion  bombardment,  the  packing  density  of  films  have  been  increased  some¬ 
what. 

Nodular  Defects 

Although  there  have  been  attempts  to  try  to  improve  the  microstruc¬ 
ture  of  thin  films  by  such  processes  as  ion  bombardment  and  UV  irradia¬ 
tion,  defects  within  the  film  still  arise  due  to  growth  asperities  (11- 
13,17,27,31).  In  thin  film  circles  these  defects  are  usually  referred 
to  as  nodules.  Put  simply,  a  nodule  is  a  structure  protruding  from  the 
thin  film  surface  in  a  domelike  manner.  A  cross  section  of  a  nodule 
shows  that  their  shape  is  like  that  of  an  inverted  cone,  either  parabo¬ 
lic  or  linear,  increasing  their  diameter  above  the  tip  of  the  cone.  See 
figure  4  for  an  example  of  the  general  shape  of  a  nodule.  Micrographs 
of  multi1ayer  stacks  (11,13,17)  show  that  the  individual  layering  se¬ 
quence  is  maintained  within  the  stack  even  though  a  nodule  may  be  pre¬ 
sent.  This  evidence  has  led  some  to  believe  that  nodules  are  growth  de¬ 
fects  and  not  massive  spatter  particles  or  empty  bubbles.  Therefore  a 


Figure  4.  Shape  of  Nodule 


nodule  can  be  thought  of  as  a  reproduction  of  an  Individual  minute  par¬ 
ticulate.  The  problem  with  nodules  occurring  in  films  is  basically  two¬ 
fold.  First,  the  non-uniform  structure  that  these  nodules  present  af¬ 
fect  the  properties  of  the  film.  Second,  nodules  have  a  tendancy  to  se¬ 
parate  from  the  film,  leaving  behind  a  hole  which  can  also  cause  prob- 
blems  with  some  of  the  film  properties. 

Substrate.  One  cause  of  nodules  is  surface  asperities.  Even 
though  most  substrates  are  thought  to  be  "flat",  there  usually  exists 
some  defects  with  the  substrate.  These  may  arise  due  to  surface  rough¬ 
ness  caused  by  insufficient  polishing,  or  may  even  be  caused  by  small 
polishing  grains  left  over  due  to  the  polishing  process  (10).  Whatever 
may  be  the  cause  of  these  small  defects,  there  is  much  evidence  leading 
to  the  belief  that  surface  asperities  can  and  do  cause  the  formation  of 
nodules. 

Impurities.  Another  cause  of  nodular  defects  within  films  is  due 
to  impurities  which  get  trapped  within  the  microstructure  during  deposi¬ 
tion.  These  impurities  may  vary  from  hydrogen  molecules  from  oil  diffu¬ 
sion  pumps  as  well  as  nitrogen,  oxygen  or  water  vapor  which  may  still  re¬ 
side  in  the  vacuum  chamber  (10).  Another  cause  of  impurities  may  even 
be  large  spatterings  of  the  material  being  used  to  produce  the  film  coat¬ 
ing.  The  effect  that  these  impurities  have  within  the  microstructure  is 
to  provide  the  deposition  material  an  asperity  to  grow  on. 

Intentional  doping  of  impurities  into  films  have  shown  a  reduction 
in  the  tensile  strength.  This  is  due  to  the  fact  that  the  doped  impuri¬ 
ties  form  grain  boundaries  which  affect  the  surface  energy.  This  in 
turn  influences  the  internal  mechanical  stress  of  the  thin  films  (11). 


Multilayers 

The  use  of  multilayers  in  thin  film  optics  is  primarily  due  to  the 
interference  effects  caused  by  layers  which  are  a  quarter  of  a  wave¬ 
length  in  thickness.  Although  most  of  the  properties  caused  by  this  phe¬ 
nomenon  are  known,  the  growth  patterns  and  the  microstructure  making  up 
these  layers  is  still  not  completely  understood. 

By  the  use  of  electron  mircrographs,  the  layered  structured  within 
multilayered  films  has  been  seen  along  with  a  columnar  microstructure 
making  up  the  individual  layers  (10,13,17).  However,  these  same  mircro¬ 
graphs  have  shown  that  there  are  obvious  discontinuities  at  the  inter¬ 
faces  between  these  layers.  This  effect  has  caused  some  researchers  to 
investigate  the  water  absorption  problems  caused  by  these  discontinu¬ 


ities  at  the  boundaries  (31). 


III.  Modelim 


The  purpose  of  this  stud;  was  to  improve  upon  a  previous  computer 
model  of  thin  film  growth  by  vapor  deposition.  In  order  to  accomplish 
this  task  several  concepts  as  to  how  the  model  would  work  had  to  be 
thought  out.  Some  of  these  included  how  the  molecules  would  be  repre¬ 
sented,  the  movement  of  the  molecules  after  collision  and  the  final 
placement  of  the  molecules  within  the  film.  In  an  effort  to  gain  an 
insight  on  how  these  questions  could  be  answered,  past  models  in  the 
literature  were  reviewed. 

Past  Models 

Although  thin  films  by  vapor  deposition  have  been  around  for  many 
years,  the  process  by  which  films  grow  is  not  clearly  understood.  There¬ 
fore,  since  computers  have  been  used  increasingly  in  the  scientific  com¬ 
munity,  people  have  used  computers  to  model  thin  film  growth  (1-3, 5, 9- 
11,15,16,24-26,28,30,31,42,46,47,49).  These  models  fall  into  two  main 
categories,  with  the  first  being  nucleation  models  and  the  second  belong¬ 
ing  to  a  group  hereafter  referred  to  as  the  columnar  category.  Although 
nucleation  theory  may  be  a  relevant  concept,  this  report  deals  with  the 
columnar  nature  of  thin  film  growth,  and  therefore  nucleation  models 
will  not  be  discussed. 

Of  all  the  columnar  growth  models  of  interest,  the  most  notable  is 
the  one  produced  by  Dirks  and  Leamy  (5).  In  this  model,  a  two-dimen¬ 
sional  simulation,  both  the  aspects  of  shadowing  and  atomic  relaxation 
are  discussed  along  with  some  aspects  of  three  dimensional  models. 

Mobility.  Whatever  is  used  to  model  the  actual  molecules  making 


up  Che  film  (disks  are  usually  used  in  two  dimensions  and  spheres  in 
three  dimensions),  the  amount  of  mobility  given  to  the  "molecules"  is 
the  most  important  aspect  in  the  development  of  the  simulation.  Dirks 
and  Leamy  showed  that  in  order  to  approach  the  proper  modeling  of  vapor 
deposited  films,  the  molecules  (i.e.  disks  or  spheres)  must  possess  some 
sort  of  limited  mobility.  This  can  be  explained  by  looking  at  the  two 
possible  extremes.  In  the  case  of  no  mobility,  the  simulated  "film" 
produced  by  Dirks  and  Leamy  was  extremely  porous  and  showed  no  columnar 
structure,  but  was  rather  composed  of  a  chain-like  structure.  The  other 
extreme  would  be  when  the  molecules  possessed  infinite  mobility,  meaning 
that  the  molecules  would  continue  to  move  until  they  found  the  point  in 
the  film  which  possessed  the  lowest  potential  energy.  This  of  course 
would  produce  a  very  dense  film,  hexagonal  closed  packed  in  two  dimen¬ 
sions,  and  was  not  even  considered  by  Dirks  and  Leamy.  Therefore  they 
assumed  a  "limited  mobility"  which  is  somewhere  in  between  the  two  ex¬ 
tremes.  More  specifically,  in  the  case  of  two  dimensions  the  incident 
disks  making  up  the  evaporant  were  allowed  to  relax  into  the  nearest 
"saddle"  or  "pocket"  made  by  two  disks.  In  the  case  of  three  dimen¬ 
sions  this  pocket  was  composed  of  three  spheres.  The  films  produced  by 
this  limited  mobility  were  composed  of  columns  three  to  five  disks  in 
diameter.  Since  then  other  models  have  used  this  same  limited  mobility 
in  their  simulations  (10,11,16,24,26,28,31,46,47,49). 

Defects.  Other  than  trying  to  produce  columnar  structures  within 
simulated  films,  computer  models  have  been  used  to  understand  how  nod¬ 
ules  are  formed  within  the  mircrostructure  and  how  this  affects  the 
growth  of  the  film.  The  most  well  known  works  concerning  the  modeling 


of  defects  have  been  by  Karl  Guenther  (11,15).  His  models  have  been 
able  to  produce  nodular  growths  in  films  by  both  substrate  imperfections 
and  the  doping  of  impurities  within  the  microstructure.  As  far  as  sub¬ 
strate  imperfections  are  concerned,  this  has  been  accomplished  by  modify¬ 
ing  the  substrates  in  such  a  way  as  to  produce  bumps  and  other  surface 
asperities.  When  it  comes  to  simulating  the  doping  of  impurities  into 
the  film,  large  spatterings  of  approximately  three  times  the  size  of  the 
simulated  molecules  have  been  used. 

Form  Birefringence.  In  addition  to  simulating  columnar  structure, 
three  dimensional  models  have  been  used  to  study  form  birefringence  (26). 
Put  simply,  films  deposited  at  oblique  incidence  produce  variations  in 
the  index  of  refraction  due  to  the  density  of  the  films  being  different 
in  various  directions,  a  phenomenon  caused  by  the  columnar  structure. 

Assumptions 

After  a  study  of  the  past  models  had  been  reviewed,  it  was  then  ne¬ 
cessary  to  make  some  simplifying  assumptions  in  order  to  produce  the  new 
model.  These  are  provided  below  along  with  any  reasons  as  to  why  the 
assumptions  were  made: 


1.  The  deposition  of  thin  films  can  be  represented  by  a  two-dimen¬ 
sional  computer  model.  Although  a  three-dimensional  model  may  be  more 
applicable  or  realistic,  it  is  still  possible  to  learn  new  things  from  a 
two-dimensional  model. 

2.  The  incident  particles  (atoms  or  molecules)  making  up  the  evaporant 
were  assumed  to  be  disks  of  a  particular  diameter.  Of  the  three  differ¬ 
ent  size  disks  in  the  model,  the  smaller  two  were  used  to  simulate  the 


evaporant.  The  largest  disk,  being  approximately  eight  times  the  size 
of  the  smallest,  was  used  to  simulate  an  impurity  being  trapped  in  the 
mircrostructure.  Although  the  size  of  the  molecules  may  differ  from  one 
to  the  next  in  reality,  the  order  of  magnitudes  are  approximately  cor¬ 
rect  for  use  in  the  model. 

3.  Each  disk,  independent  of  size,  was  assumed  to  travel  on  a  straight 
line  at  an  angle  A  relative  to  the  substrate  normal.  It  continued  on 
this  trajectory  until  coming  into  contact  with  one  of  the  already  depos¬ 
ited  disk  or  the  substrate.  This  is  in  good  agreement  with  what  ac¬ 
tually  happens  in  vapor  deposition.  The  reason  being  that  the  deposi¬ 
tion  takes  place  in  a  vacuum  plant,  and  there  are  few  if  any  molecules 
removed  from  their  trajectory. 

4.  Each  individual  disk  was  assumed  to  be  deposited  serially.  Al¬ 
though  the  actual  deposition  rate  of  the  molecules  may  be  such  that  par¬ 
allel  activities  take  place  (i.e.  two  molecules  arriving  at  the  sub¬ 
strate  at  the  same  time),  the  probability  of  these  two  events  being 
close  enough  to  affect  one  another  is  rare. 

5.  It  was  assumed  that  after  an  incident  disk  came  into  contact  with 
an  already  deposited  disk,  it  still  had  a  limited  mobility.  This  mobil¬ 
ity  allowed  the  disks  to  come  to  rest  in  one  of  two  ways.  First,  the 
incident  disk  remained  in  contact  with  the  first  disk  and  then  relaxed 
into  the  nearest  "pocket”  until  it  made  contact  with  another  previously 
deposited  disk.  Second,  a  set  percentage  of  the  incident  particles  were 
given  an  extra  mobility  allowing  them  to  move  to  a  pocket  farther  away 
but  still  maintaining  contact  with  the  first  disk.  The  reason  these 
types  of  mobility  were  introduced  into  the  model  was  to  find  the  amount 


of  limited  mobility  needed  to  produce  realistic  films,  as  discussed  in 
the  previous  chapter. 


Further  discussion  of  why  some  of  these  assumptions  were  made  and 
how  they  were  used  in  the  model  can  be  found  in  the  following  sections. 

Disk  Size  and  Shape 

After  reviewing  the  literature  it  became  evident  that  the  most  com¬ 
mon  shape  used  to  represent  molecules  was  a  sphere  for  three-dimensional 
models  and  circles/disks  for  two  dimensions.  Even  though  these  are  rela¬ 
tively  simple  shapes  compared  to  molecules  actually  used  in  vapor  depo¬ 
sition,  they  are  relatively  easy  to  represent  mathematically.  It  is  for 
this  reason,  and  the  fact  that  the  old  model  used  disks,  that  this  model 
also  used  disks  to  represent  the  molecules. 

The  size  of  the  disk  used  to  represent  the  molecules  in  the  film 
was  determined  by  looking  at  the  two-dimensional  arrays  used  to  hold  the 
disks.  More  specifically,  three  arrays  (each  300  wide  and  200  in 
height),  were  used  to  represent  the  film,  accounting  for  60000  unit 
cells.  The  first  array  was  used  to  store  occupancy  data  of  the  x-y 
field  (i.e.  whether  or  not  a  unit  cell  was  occupied).  The  second  and 
third  arrays  were  used  to  store  the  x  coordinates  and  y  coordinates  of 
the  disks  respectively.  In  order  to  avoid  the  problem  of  two  or  more 
disks  occupying  a  single  cell,  the  diameter  of  the  smallest  disk  was 
equal  to  the  square  root  of  two  (see  figure  5).  Since  the  model  was  to 
simulate  multilayered  coatings,  a  second  disk  of  diameter  equal  to  two 
times  the  square  root  of  two  was  used  to  simulate  a  larger  molecule  or 
one  with  a  greater  interaction  distance.  Finally,  in  order  to  simulate 


impurities  a  third  disk  was  chosen  to  be  eight  times  the  square  root  of 
two  in  diameter  or  eight  times  the  size  of  the  smaller  disk.  See  figure 
6  for  a  comparative  sketch  of  the  three  different  disk  sizes. 

Disks  are  represented  in  the  computer  by  the  radius  of  the  individ¬ 
ual  disks  and  its  center  coordinates.  Once  the  address  of  a  cell  has 
been  determined  by  truncating  its  coordinates,  the  information  is  stored 
in  the  arrays.  In  the  occupation  array  a  zero  represents  that  the  cell 
is  unoccupied;  while  a  one,  two  or  eight  specify  the  size  of  disk  with 
its  center  located  inside  the  particular  cell. 

Contact  Circles 

As  an  incident  disk  approaches  those  that  have  been  previously  de¬ 
posited  the  position  of  the  incident  disk  when  it  first  collides,  as 
well  as  its  final  resting  position,  needs  to  be  determined.  In  order  to 
accomplish  this  task  the  idea  of  contact  circles  from  the  old  model  was 
invoked.  Put  simply,  the  contact  circle  of  an  individual  disk  is  a  cir¬ 
cle  about  the  disk,  whose  radius  is  such  that  the  edge  of  the  circle  pas¬ 
ses  through  the  center  of  the  neighboring  disk  when  they  are  in  contact. 
An  example  of  two  contact  circles  can  be  seen  in  figure  7  where  the 
smaller  disk  is  in  incident  on  the  larger  one.  Implied  is  the  fact  that 
it  is  possible  for  one  disk  to  have  more  than  one  size  contact  circle 
when  it  is  surrounded  by  different  size  disks. 

The  idea  of  contact  circles  is  first  used  in  the  model  in  determ- 
ing  the  position  of  an  incident  disk  after  collision.  As  an  incident 
disk  approaches  an  already  deposited  disk  or  the  substrate,  which  is  com¬ 
posed  of  disks  (a  matter  discussed  in  chapter  4),  the  position  of  the 
disk  immediately  upon  collision  can  be  determined  mathematically.  More 
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equations  used  to  find  the  x  and  y  intercepts  are 


[(1/m  )  +  1] 


X  -  Y/m  +  a 


where  the  significance  of  the  variables  as  well  as  the  derivation  of  the 
equations  can  be  found  in  appendix  A. 

The  second  place  where  the  idea  of  contact  circles  comes  into  play 
is  determining  the  final  resting  point  of  the  incident  disk.  As  stated 
above,  the  incident  disk  will  come  to  rest  in  a  "pocket"  formed  by  the 
disk  in  which  it  first  makes  contact,  known  as  the  collision  disk,  and  a 
neighboring  or  rest  disk.  In  order  to  calculate  the  position  of  the  cen 
ter  of  the  disk  after  it  comes  to  rest,  the  contact  circles  of  the  col¬ 
lision  and  rest  disk  (whose  radii  are  determined  by  the  radius  of  the  in 
cident  disk  and  their  own  respectively)  are  used  as  in  figure  9.  By 
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Figure  9.  Rest  Point  Contact  Circles 


finding  the  intersection  of  the  two  contact  circles  the  rest  point  can 
be  determined.  The  mathematical  formulas  used  to  calculate  the  rest 
coordinate  are 


Y  -  [-a*k  +/-  (a2*k2  -  c*b)]/c 
X  =  [-a*h  +/-  (a2*h2  -  c*b)]/c 


where  the  significance  of  the  variables  and  the  derivations  of  the  eq¬ 
uations  can  be  found  in  appendix  A.  Because  these  formulas  actually 
specify  two  distinct  points,  the  point  which  allows  the  disk  to  move 
the  least  is  chosen  as  the  rest  point,  except  in  the  case  of  added  mo¬ 
bility;  a  point  which  will  be  discussed  in  the  next  section. 


Mobility 

As  stated  earlier  in  the  section  concerning  past  models,  the  mobil¬ 
ity  of  the  disks  after  initial  collision  should  be  limited.  In  an  ef¬ 
fort  to  try  to  understand  how  much  mobility  the  molecules  should  have 
after  collision,  the  new  model  allows  the  user  to  look  at  three  differ¬ 
ent  mobilities.  The  first  type  of  mobility  that  the  disks  can  undergo 
is  no  mobility  at  all.  When  disks  undergo  no  mobility,  the  resting 
point  of  the  disk  is  determined  by  the  collision  point  coordinates  pre¬ 
viously  mentioned.  Therefore,  whenever  the  disk  first  comes  into  con¬ 
tact  with  another  disk,  the  deposition  process  for  that  disk  is  com¬ 
plete.  As  shown  in  figure  10  the  resulting  film  is  very  loosely 
packed.  The  second  type  of  mobility  possible  is  what  will  be  referred 
to  as  normal  mobility.  Under  normal  mobility,  the  final  resting  point 
of  the  incident  disk  is  determined  by  the  disk  coming  to  rest  in  the 
nearest  pocket.  As  shown  in  figure  11,  normal  mobility  implies  that 
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each  individual  disk  is  in  contact  with  at  least  two  other  disks 


The  third  type  of  mobility,  which  will  be  referred  to  as  added  or 
extra  mobility,  was  introduced  as  one  of  the  parameters  in  the  new  mo¬ 
del.  The  reason  for  introducing  this  parameter  is  to  try  to  produce 
films  with  greater  packing  densities  and  larger  individual  columns.  Un¬ 
der  extra  mobility,  the  incident  disk  comes  to  rest  in  the  pocket  which 
allows  the  disks  greater  movement  before  coming  to  rest.  As  shown  in 
figure  12,  normally  the  incident  disk  will  come  to  rest  at  point  A.  How¬ 
ever,  when  extra  mobility  is  introduced  the  incident  disk  will  come  to 
rest  at  point  B  making  the  column  wider  than  before.  More  specifically, 
under  extra  mobility  the  incident  disk  still  remains  in  contact  with 
the  collision  disk  but  just  rolls  to  the  next  available  pocket  produced 
by  surrounding  disks. 


Figure  12.  Extra  Disk  Mobility 
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Program  Implementation 


In  order  to  simulate  the  vapor  deposition  of  thin  films  by  use  of  a 
computer,  it  was  first  necessary  to  decide  what  the  program  was  to  mo¬ 
del.  That  is,  what  specific  tasks  was  the  program  going  to  be  asked  to 
do.  In  order  to  model  the  deposition  process  properly,  the  following 
major  tasks  needed  to  be  incorporated  into  the  program.  One,  the  pro¬ 
gram  had  to  be  able  to  produce  a  substrate,  and  to  simulate  an  actual 
substrate  it  was  essential  that  it  be  able  to  produce  imperfections. 

Two,  the  program  needed  to  deposit  a  film  onto  the  substrate,  and  in  do¬ 
ing  so  the  deposition  process  needed  to  contain  parameters  that  could  be 
changed  to  simulate  different  deposition  conditions.  Three,  in  order  to 
analyze  the  results,  the  program  needed  to  have  an  analysis  section  for 
calculating  such  variables  as  packing  density  and  angle  of  growth  by  the 
film.  And  lastly,  the  results  of  the  deposition  needed  to  be  plotted, 
so  a  fourth  task  of  the  program  was  to  make  possible  the  plotting  of  the 
deposited  film  and  substrate. 

The  program  listed  in  appendix  B  is  able  to  take  care  of  the  four 
tasks  by  means  of  four  subroutines.  Each  subroutine  is  connected  to  the 
other  by  a  main  interface  loop.  The  names  of  the  four  subroutines  in 
the  program  are  simply  the  substrate,  deposition,  analysis  and  moving 
subroutines. 

Substrate  Subroutine 


The  substrate  subroutine  is  divided  into  five  areas  allowing  for 
the  initialization,  creation,  storing,  recalling  and  transferring  of  sub¬ 
strates.  The  movement  between  these  sections  is  controlled  by  an  inter- 


face  loop  at  the  beginning  of  the  subroutine. 

In  an  effort  to  keep  the  user  from  destroying  valuable  information, 
the  subroutine  uses  a  set  of  buffer  arrays  for  manufacturing  the  sub¬ 
strate.  These  arrays  (namely  xc,  yc  and  subbuf),  each  of  which  is  300 
in  width  and  10  in  height,  are  used  to  create  the  substrate  before  be¬ 
ing  transferred  to  the  main  arrays.  More  specifically,  the  information 
contained  in  these  arrays  hold  information  concerning  the  x  coordinate, 
y  coordinate  and  occupation  (i.e.  whether  or  not  a  cell  is  occupied)  of 
the  3000  cells  used  to  produce  the  substrate. 

As  in  the  previous  model  the  substrate  is  composed  of  disks,  each 
of  which  is  the  square  root  of  two  in  diameter.  In  order  to  produce  the 
substrate  these  disks  are  laid  down  either  one  on  top  of  the  other  or 
side  by  side  so  that  each  disk  just  comes  into  contact  with  its  neigh¬ 
bor.  See  figure  13  for  an  example  of  how  these  disks  can  be  stacked. 


Because  the  arrays  are  300  wide  and  the  diameter  of  each  disk  is  the 
square  root  of  two,  a  total  of  213  disks  can  be  laid  across  the  bottom 
of  the  array  thus  producing  a  "flat”  substrate.  See  figure  14  for  an 
example.  Although  this  may  not  appear  to  be  flat  due  to  the  shape  of 
the  disks,  it  must  be  remembered  that  some  optical  substrates  are  no¬ 
thing  more  than  cleaved  NaCl  which  is  not  "flat"  on  an  atomic  scale. 

In  order  to  produce  a  more  complicated  substrate,  that  is  one  which 
contains  imperfections,  the  remaining  parts  of  the  buffers  are  used. 
Since  the  arrays  are  10  high,  the  number  of  disks  that  can  be  stacked 
one  upon  another  comes  to  7,  implying  that  1491  disks  can  be  used  to 
create  the  substrate.  The  fact  that  there  are  213  disks  across  the  bot¬ 
tom  of  the  arrays  is  used  to  envision  213  columns  within  the  substrate, 
with  each  column  having  a  maximum  possible  height  of  7  disks.  By  spec¬ 
ifying  the  height  in  each  column  it  is  then  possible  to  produce  sub¬ 
strates,  as  seen  in  figures  15  through  20. 

Since  the  program  was  intended  to  be  user  friendly,  it  is  not  re¬ 
quired  that  the  user  enter  the  height  in  each  column  to  produce  a  sub¬ 
strate.  In  fact,  after  initializing  the  buffers  and  entering  into  the 
creation  part  of  the  subroutine,  a  flat  substrate  like  the  one  in  fig¬ 
ure  14  is  automatically  produced.  Furthermore,  if  the  user  desires  a 
more  complicated  substrate  and  this  substrate  contains  a  section  where 
the  height  remains  constant  over  several  columns,  the  user  is  spared  the 
trouble  of  entering  the  same  height  over  and  over. 

After  entering  the  height  in  a  particular  column,  the  program  asks 
whether  or  not  the  height  is  to  be  continued  and  if  so  to  what  column 
this  should  occur.  An  example  of  how  the  creation  of  a  substrate  may 


take  place  is  given  below  with  the  prompts  and  appropriate  answers. 

Do  you  wish  to  create  more  than  just  a  horizontal  line  for  the 
substrate? 

If  yes  -  Enter  1 
Otherwise  -  Hit  Return 

1 

There  are  213  columns  columns  across  the  bottom  of  the  substrate. 
Each  one  of  these  must  have  a  height  specified  to  it.  The  follow 
ing  steps  will  help  in  that  process. 

You  are  in  column  1 

Enter  height  (must  be  7  or  less) 

1 

Do  you  wish  this  height  (1)  to  continue? 

If  yes  -  Enter  1 
Otherwise  -  Hit  Return 

1 

To  what  column  should  this  continue? 

(must  be  213  or  less) 

100 

You  are  in  column  101 

Enter  height  (must  be  7  or  less) 

7 

Do  you  wish  this  height  (5)  to  continue? 

If  yes  -  Enter  1 
Otherwise  -  Hit  Return 

1 

To  what  column  should  this  continue? 

(must  be  213  or  less) 

105 

You  are  in  column  106 

Enter  height  (must  be  7  or  less) 


Do  you  wish  this  height  (1)  to  continue? 

If  yes  -  Enter  1 
Otherwise  -  Hit  Return 

1 

To  what  column  should  this  continue? 

(must  be  213  or  less) 

213 

The  substrate  is  now  complete ! ! 

Ready  for  transfer  to  film  matrices. 

A  sample  of  the  substrate  just  created  by  the  above  commands  can  be 
seen  in  figure  15. 

Once  the  substrate  has  been  produced,  whether  it  is  just  flat  or 
some  imperfection  has  been  introduced,  the  user  may  wish  to  store  the 
values  of  the  substrate  for  future  use.  This  is  accomplished  by  the 
third  section  of  the  subroutine.  The  fourth  section  is  the  counter¬ 
part  of  the  previous  section  and  allows  the  user  to  recall  substrates 
so  that  they  do  not  have  to  be  reproduced.  The  final  section  of  the 
subroutine  takes  care  of  transferring  the  values  of  the  buffer  arrays 
into  the  main  arrays  used  in  the  deposition  of  the  film.  Once  this  has 
occurred  the  user  may  enter  the  deposition  subroutine  and  prepare  to  de¬ 
posit  the  film  onto  the  substrate. 

Deposition  Subroutine 

The  deposition  subroutine  is  divided  into  three  sections,  allowing 
for  the  initialization  of  the  main  arrays,  the  setting  of  the  deposition 
variables  and  lastly  the  deposition  of  the  film  onto  the  substrate.  As 
in  the  substrate  subroutine,  movement  between  these  sections  is  con¬ 
trolled  by  an  interface  loop  located  at  the  beginning  of  the  subroutine. 
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After  the  initialization  of  the  main  arrays  takes  place,  the  user 
of  the  program  is  then  able  to  set  the  deposition  variables  which  con- 
control  the  deposition  process.  The  first  parameter  specified  by  the 
user  is  the  number  of  disks  to  be  deposited.  In  order  to  look  at  chan¬ 
ges  in  the  overall  packing  density  of  the  film,  the  user  should  enter 
more  disks  than  desired.  The  reason  for  this  is  that  deposition  of  the 
film  continues  until  either  all  the  disks  have  been  used  or  until  the 
arrays  become  full.  Therefore,  by  specifying  more  disks  than  the  arrays 
can  hold,  the  effect  that  other  parameters  have  on  the  packing  density 
can  be  measured. 

The  second  parameter  that  the  user  has  control  of  is  the  main  angle 
of  deposition.  The  angle  entered  by  the  user  is  the  angle  as  measured 
from  the  substrate  normal.  Therefore  if  zero  degrees  were  entered,  de¬ 
position  would  occur  at  normal  incidence.  In  order  to  simulate  the  ro¬ 
tation  of  substrates,  the  user  also  has  the  ability  to  vary  the  angle 
of  incidence  as  well  as  the  speed  at  which  the  variation  takes  place. 
This  is  accomplished  by  the  user  entering  the  angle  which  corresponds  to 
the  total  amount  of  variation  as  in  figure  21.  If  the  main  angle  of 


incidence  is  at  35  degrees  and  the  amount  of  variation  is  specified  to 
be  10  degrees,  then  deposition  will  occur  between  30  and  40  degrees. 

The  speed  at  which  the  angle  changes  is  determined  by  the  size  of  the  an¬ 
gle  increment,  implying  that  the  smaller  the  angle  increment  the  slower 
the  angle  changes  over  time.  An  example  of  how  these  variables  are  en¬ 
tered  occurs  below  with  the  prompts  and  possible  answers  given. 

At  what  angle  do  you  wish  the  film  to  be  deposited? 

Please  enter  angle  in  degrees  -  Enter  in  real  format. 

27.0 

Do  you  wish  to  vary  the  angle  of  incidence? 

If  yes  -  Enter  1 
Otherwise  -  Hit  Return 

1 

How  many  degrees  (full  angle)  do  you  wish  to  vary  the  angle? 

Enter  degrees  in  real  format. 

6.0 

What  size  angle  do  you  wish  to  increment  the  angle  of  incidence? 

Enter  degrees  in  real  format. 

.01 

With  the  above  answers  given  to  the  prompts,  the  main  angle  of  deposi¬ 
tion  will  occur  at  27  degrees  and  will  vary  between  24  and  30  degrees 
with  the  angle  being  changed  0.01  degrees  after  every  disk. 

The  next  set  of  parameters  that  the  user  is  able  to  specify  con¬ 
cerns  the  amount  of  mobility  the  disks  have  after  making  initial  contact 
with  a  previously  deposited  disk  or  the  substrate.  The  user  has  the 
choice  of  three  different  disk  mobilities.  First,  the  user  may  option 
to  let  the  disk  have  no  mobility.  This  of  course  would  produce  a  film 


with  a  very  low  packing  density  like  that  of  Dirks  and  Leamy  (5),  and 
does  not  simulate  real  films.  Second,  normal  mobility  may  be  chosen  in 
which  the  disks  come  to  rest  in  a  "pocket”  produced  by  two  other  disks. 
The  final  option  the  user  may  choose  is  allowing  a  predetermined  percen¬ 
tage  of  disks  to  undergo  an  extra  mobility  as  discussed  in  chapter 
three.  The  percentage  is  determined  by  the  user  entering  the  number  of 
disks  to  be  deposited  before  an  extra  mobility  iteration  occurs.  There¬ 
fore  if  the  user  wishes  20  percent  of  the  disks  to  undergo  an  extra  mo¬ 
bility,  the  number  5  is  entered. 

The  fourth  parameter  the  user  enters  determines  whether  or  not  a 
single  impurity,  the  size  of  which  was  discussed  in  chapter  three,  is  to 
be  deposited  within  the  mircrostructure.  If  the  user  does  decide  to  im¬ 
plant  an  impurity,  the  placement  of  the  impurity  is  determined  randomly 
by  the  computer. 

The  last  set  of  parameters  the  user  has  control  over  is  the  size  of 
the  disk  to  be  deposited.  The  user  may  decide  to  deposit  only  one  size 
disk,  either  small  or  large,  or  he  may  decide  to  create  a  multilayered 
film.  If  the  user  decides  upon  creating  a  multilayered  film,  the  pro¬ 
gram  then  prompts  the  user  which  size  disk  is  to  be  deposited  first  and 
then  the  depth  of  the  layers  in  terms  of  the  film  arrays. 

After  the  deposition  variables  are  set,  the  deposition  of  the  film 
onto  the  designed  substrate  may  be  invoked.  It  should  be  noted  that  the 
five  major  parameters,  number  of  disks,  the  angle  of  deposition,  mobil¬ 
ity  of  the  disks,  deposition  of  impurities  and  the  size  of  disks,  are 
all  independent  of  one  another.  Therefore,  any  combination  of  these 
parameters  may  take  place. 


Deposition  of  film.  The  deposition  of  the  film  onto  the  substrate 
is  the  heart  of  the  entire  program.  It  is  this  section  of  the  program 
that  all  other  sections  were  designed  around.  The  section  itself  is  di¬ 
vided  into  three  areas  which  are  responsible  for  calculating  collision 
partners  and  coordinates,  rest  point  coordinates  and,  when  needed,  the 
rest  point  coordinates  for  disks  which  undergo  an  extra  mobility. 

Collision  Coordinates.  When  an  incident  disk  approaches  a  pre¬ 
viously  deposited  disk  or  the  substrate,  it  is  necessary  to  know  with 
what  disk  it  makes  contact  and  the  position  of  the  incident  disk  when 
this  occurs.  The  reason  for  wanting  to  know  the  position  of  the  inci¬ 
dent  disk  upon  collision  comes  about  when  finding  the  proper  rest  disk. 
In  order  to  accomplish  this  task  the  collision  corridor  method  used  in 
the  previous  thesis  (47)  was  utilized  with  some  modifications  since  mul¬ 
tiple  sized  disks  were  used  in  the  new  model. 

Put  simply,  the  collision  corridor  method  uses  parallel  streamers 
to  search  the  area  around  an  incident  disk  trajectory  for  possible  col¬ 
lision  partners.  Each  individual  streamer  scans  until  reaching  a  cell 
not  previously  encountered  by  another  streamer.  Once  an  occupied  cell 
is  found  the  y  coordinate  of  the  incident  disk  is  calculated  by  the  al¬ 
gorithm  in  chapter  three.  Once  each  streamer  has  preformed  its  scan, 
the  collision  partner  which  yields  the  largest  incident  disk  y  coordi¬ 
nate  is  chosen  as  the  collision  disk.  The  x  coordinate  of  the  incident 
disk  is  then  found  by  using  the  y  value  just  obtained. 

The  number  of  streamers  necessary  to  define  the  collision  corri¬ 
dor  is  dependent  upon  two  parameters.  First,  the  spacing  between  the 
individual  streamers  needs  to  be  such  that  no  matter  what  the  angle  of 


deposition,  no  unit  cell  can  be  placed  between  them.  Second,  the  width 
of  the  corridor  has  to  be  large  enough  that  all  unit  cells  capable  of 
containing  collisional  partners  are  checked.  This  second  parameter  is 
therefore  dependent  on  the  size  of  the  incident  disk,  as  well  as  the 
size  of  the  previously  deposited  disks.  As  an  example,  the  case  which 
calls  for  the  least  number  of  streamers  occurs  when  only  the  smallest 
disks  are  being  deposited.  In  this  case  for  a  collision  to  occur  the 
center  of  the  incident  disk  needs  to  be  one  disk  diameter  away  from  the 
center  of  another  disk;  namely,  one  square  root  of  two.  The  minimum 
number  of  streamers  necessary  to  check  around  the  incident  disk  in  this 
particular  example  is  four  (as  shown  in  figure  22).  When  larger  disks 
are  deposited  the  number  of  streamers  increases,  as  in  the  case  when  16 
streamers  are  used  when  an  impurity  is  incident  on  large  disks. 


TRAJECTORY 


STREAMER 


UNIT  CELL 


Figure  22.  Example  of  Streamers 


Rest  Point  Coordinates.  Once  a  collision  disk  is  determined  and 
the  position  of  the  incident  disk  upon  collision  known,  the  final  rest¬ 
ing  place  of  the  disk  can  be  determined.  As  mentioned  previously,  the 
user  of  the  program  has  three  options  as  far  as  the  mobility  of  the 
disks  are  concerned.  With  no  relaxation  of  the  disks,  it  is  implied 
that  the  collision  coordinates  of  the  incident  disk  become  the  rest 
point  coordinates.  For  the  case  of  extra  mobility,  this  subject  will  be 
presented  in  the  next  area. 

When  normal  relaxation  of  disks  occurs,  the  final  resting  place  is 
determined  by  the  mathematical  derivation  explained  in  chapter  three. 
That  is,  find  a  third  disk  or  rest  disk  with  which  the  incident  disk 
rolls  into  contact  with,  while  still  maintaining  contact  with  the  col¬ 
lision  disk.  The  program  algorithm  used  to  accomplish  this  task  in¬ 
volves  using  an  array  search  around  the  collision  disk.  The  rest  disk 
is  then  found  by  determining  which  disk  allows  the  incident  disk  to 
move  the  shortest  distance  before  coming  to  rest.  If  no  further  relax¬ 
ation  is  to  occur  the  rest  position  of  the  incident  disk  is  then  re¬ 
corded  in  the  main  arrays  by  truncating  the  final  position  coordinates. 

Extra  Mobility  Coordinates.  In  order  to  simulate  some  disks  under¬ 
going  more  than  "normal  mobility",  the  third  area  in  the  deposition  sec¬ 
tion  calculates  the  rest  point  coordinates  of  these  disks.  In  essence, 
this  part  of  the  section  finds  the  rest  coordinates  for  the  incident 
disks,  which  allows  it  to  move  to  the  "next"  shortage  distance  in  compar 
ision  to  that  of  normal  mobility.  Like  normal  mobility,  the  incident 
disk  is  allowed  to  remain  in  contact  with  the  collision  disk. 

Although  the  algorithm  for  normal  mobility  may  be  simple,  the  task 


of  finding  an  alternate  rest  place  is  more  difficult.  Since  the  inci¬ 
dent  disk  is  to  maintain  contact  vith  the  collision  disk  and  the  rest¬ 
ing  place  for  normal  mobility  is  to  be  ignored,  there  exists  only  one 
other  possible  resting  position.  Due  to  the  way  the  algorithm  was  de¬ 
rived,  the  second  resting  place  may  or  may  not  correspond  to  the  posi¬ 
tion  which  allows  the  incident  disk  to  move  the  next  shortest  distance 
As  shown  in  figure  23  (a  condition  immediately  after  contact  has  been 
made),  position  1  is  where  the  incident  disk  would  come  to  rest  under 
the  condition  of  normal  mobility.  Although  positions  2  through  5  cor¬ 
respond  to  places  where  the  algorithm  would  designate  as  possible  rest 
points,  these  are  either  occupied  or  physically  impossible  due  to  sur¬ 
rounding  disks.  Position  6  on  the  other  hand  is  the  second  of  the  two 
possible  resting  points.  In  order  to  arrive  at  this  decision  using  a 
computer,  the  extra  mobility  routine  uses  a  set  of  seven  stored  posi- 


tions  and  then  tests  whether  each  is  occupied  or  physically  plausible. 
Immediately  after  the  testing  of  the  seven  different  positions,  an¬ 
other  test  looks  at  the  y  coordinate  of  the  extra  mobility  position  in 
comparison  with  the  y  coordinate  determined  for  normal  relaxation.  If 
the  extra  mobility  position  is  significantly  higher  than  that  for  normal 
relaxation  the  position  is  considered  invalid.  The  reason  for  doing 
this  is  to  ensure  that  the  resulting  film  becomes  more  dense  rather  than 
less  dense.  In  the  case  where  none  of  the  seven  positions  are  found  to 
be  suitable  and/or  the  height  position  test  fails,  the  final  position  of 
the  disk  is  decided  by  default  to  be  the  position  determined  by  normal 
relaxation.  The  next  disk,  in  order  to  make  up  for  the  failure  of  the 
test,  is  then  given  extra  mobility.  This  continues  until  a  suitable 
disk  is  found. 

Periodic  Boundaries.  In  order  to  save  on  space  within  the  memory 
of  the  computer,  periodic  boundaries  were  incorporated  into  the  model. 
More  specifically,  this  condition  allowed  incident  disks  which  would 
normally  disappear  off  the  left  side  to  reappear  on  the  right  and  visa 
versa.  The  result  of  introducing  periodic  boundaries  can  be  seen  in 
in  the  next  chapter  where  the  figures  give  the  impression  that  one  sec¬ 
tion  of  the  film  is  being  viewed. 


Analysis  Subroutine 

After  a  film  has  been  produced  it  is  necessary  to  characterize  its 
properties.  In  this  subroutine  two  properties  of  the  simulated  film  are 
analyzed  by  two  separate  sections,  namely  the  packing  density  and  angle 
of  growth  of  the  film.  The  movement  between  these  sections  is  con¬ 
trolled  by  an  interface  loop  as  in  the  first  two  subroutines. 


Packing  Density  Analysis.  As  mentioned  in  chapter  one  when  the  den¬ 
sity  of  the  film  increases,  the  properties  of  the  film  approach  those  of 
the  bulk  material.  This  is  desirable  because  the  optical  properties  of 
the  film  become  more  stable.  An  example  would  be  the  reduction  of  water 
absorption  due  to  increased  packing  density,  which  in  turn  reduces  the 
amount  of  wavelength  shift.  In  order  to  monitor  the  density  of  the 
film,  the  first  section  of  this  subroutine  is  dedicated  to  that  purpose. 

Reviewing  from  chapter  two,  the  packing  density  of  the  film  is  sim¬ 
ply: 


volume  of  the  solid  part  of  the  film  (2) 

^  “  total  volume  of  the  film  (solid  +  voids) 

where  p  is  the  packing  density.  This  implies  that  the  value  of  1.00  is 
the  largest  value  the  packing  density  can  achieve.  In  order  to  make  the 
program  more  user  friendly  the  value  of  p  is  multiplied  by  100  making  it 
a  percentage. 

In  order  to  calculate  the  packing  density  of  the  film,  where  only 
one  size  disk  has  been  deposited,  the  number  of  cells  which  are  occupied 
are  counted  and  then  divided  by  the  total  number  of  cells.  This  number 
is  then  multiplied  by  the  area  of  a  disk  relative  to  a  cell,  and  then 
multiplied  by  100.  However,  in  the  case  where  different  sized  disks  are 
used,  separate  counters  are  needed  for  each  size  disk.  After  calcula¬ 
ting  the  individual  percentages  for  each  size  disk,  the  percentages  are 
then  added  to  obtain  the  total  packing  density.  With  the  three  differ¬ 
ent  sized  disks  used  in  the  program  the  individual  packing  densities 


were  calculated  as  follows: 


psmall  =  (number  of  cells  with  small  disk/60000)*l. 570796*100  (7) 

plarge  =  (number  of  cells  with  large  disk/60000)*6. 283185*100  (8) 

pimpur  =  (number  of  cells  with  impurities/60000)*100. 5309*100  (9) 

where  psmall,  plarge  and  pimpur  are  the  respective  packing  densities  of 

the  three  different  disks. 

Because  films  produced  by  vapor  deposition  tend  to  be  amorphous  in 
nature,  this  implies  that  the  packing  density  throughout  the  film  may 
not  be  uniform.  In  an  effort  to  look  for  spatial  variations  within  the 
simulated  films  the  density  algorithm  was  modified.  Not  only  is  it  pos¬ 
sible  to  find  the  total  packing  density  of  the  film,  but  spatial  varia¬ 
tions  in  both  the  horizontal  and  vertical  directions  are  possible. 
Therefore,  the  user  of  the  program  may  look  at  individual  areas  within 
the  film.  This  was  accomplished  by  simply  specifying  the  upper  and  lo¬ 
wer  heights  for  vertical  variations,  and  left  and  right  boundaries  for 
horizontal  variations.  The  values  entered  correspond  to  the  addresses 
of  the  unit  cells  of  the  arrays. 

Angle  of  Growth.  The  angle  of  columnar  growth  by  the  film  is  re¬ 
sponsible  for  many  characteristics  as  noted  in  chapter  one,  and  is  depen¬ 
dent  on  the  deposition  angle  via  the  tangent  rule  in  many  films.  In  or¬ 
der  to  calculate  the  angle  of  growth  the  algorithm  used  in  the  previous 
model  (47)  was  invoked. 

This  algorithm  calculates  the  angle  of  growth  by  looking  at  varia¬ 
tions  in  the  density  of  the  film  at  varying  angles.  In  doing  so,  the  an¬ 
gle  of  growth  can  be  found  at  different  heights  within  the  film,  amount¬ 
ing  to  vertical  variations  within  the  film. 
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Moving  Subroutine 

The  moving  subroutine  was  developed  in  order  to  store  the  values  of 
the  main  arrays  making  up  the  film  for  future  analysis  and/or  plotting. 
This  is  accomplished  by  dividing  the  subroutine  into  four  sections. 

Again,  as  in  the  previous  subroutines,  the  movement  between  the  sections 
is  controlled  by  an  interface  loop. 

The  first  two  sections  of  the  subroutine  provide  a  somewhat  inef¬ 
ficient  way  of  storing  and  retrieving  values  in  that  all  values  of  the 
main  arrays,  both  zero  and  non-zero,  are  moved  around.  Although  this 
is  not  memory  efficient  it  may  prove  to  be  useful  for  certain  types  of 
graphing  devices. 

The  third  section  of  the  subroutine,  like  the  first,  stores  the  val¬ 
ues  of  the  main  arrays  into  files;  however,  in  this  case  only  the  non¬ 
zero  values  are  stored.  More  specifically,  this  section  produces  two 
files  for  storing  values  for  plotting.  The  first  is  a  master  file  con¬ 
taining  data  on  all  the  disks  independent  of  size,  and  holds  information 
concerning  the  occupancy  of  the  unit  cells  as  well  as  the  x  and  y  coordi¬ 
nates  of  the  individual  disks.  The  second  file,  on  the  other  hand,  con- 
contains  similar  data  regarding  only  the  larger  sized  disks. 

The  fourth  section  of  the  subroutine,  somewhat  like  the  second, 
reads  back  into  the  program  the  values  stored  in  a  file  by  section 
three.  What  makes  this  section  different  from  the  second  is  that  only 
the  data  concerning  the  occupancy  of  the  unit  cells  is  read  back.  This 
in  turn  allows  for  only  an  analysis  of  the  film,  namely  density  calcula¬ 
tions  and  determination  of  the  angle  of  growth. 
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As  with  any  program,  it  is  necessary  to  analyze  the  results  and 
then  make  conclusions  based  on  the  analysis.  The  purpose  of  this  chap¬ 
ter  is  to  look  at  the  resulting  films  and  then  try  to  analyze  the  films 
based  on  the  parameters  under  which  the  films  were  deposited.  Through¬ 
out  the  last  stages  of  program  development  and  finally  during  the  run¬ 
ning  of  the  program  to  obtain  results,  approximately  35  films  were 
printed  out  and  the  results  examined.  However,  to  conserve  on  space 
only  13  of  the  most  interesting  films  will  be  discussed  in  this  chap¬ 
ter. 

In  order  to  provide  a  control  upon  which  the  films  could  be  com¬ 
pared,  5  films  were  deposited  at  various  angles.  With  these  films,  none 
of  the  deposition  parameters  were  varied,  therefore  allowing  future 
films  to  be  compared  against  them.  Three  of  these  films  can  be  seen  in 
figures  24  through  26.  The  other  films  provided  in  this  report  (figures 
27  through  36),  all  have  at  least  one  parameter  that  was  introduced.  In 
order  to  refresh  the  reader  these  parameters  are  the  temporal  variation 
of  the  angle  of  deposition,  the  mobility  of  the  disks  upon  collision, 
the  doping  of  impurities  into  the  microstructure,  the  ability  to  intro¬ 
duce  imperfections  into  the  substrate  and  the  simulation  of  multilayered 
coatings. 

Controls 

By  viewing  the  three  control  films  (figures  24  -  26),  some  of  the 
main  characteristics  exhibited  by  other  columnar  models  can  be  seen.  In 
figure  24,  a  film  deposited  under  normal  incidence,  columns  3  to  5  disks 


thick  can  be  seen  branching  in  both  the  left  and  right  directions  due 
to  the  normal  incidence.  When  the  angle  of  deposition  is  increased 
to  25  degrees  as  in  figure  25,  little  if  any  branching  takes  place  that 
is  headed  toward  the  left.  In  figure  26,  the  film  is  definitely  leaning 
towards  the  right  due  to  the  angle  of  incidence  being  50  degrees.  It 
can  also  be  seen  that  the  overall  density  of  the  film  is  less  than  that 
for  normal  incidence,  and  that  the  spacing  between  the  columns  is  sig¬ 
nificantly  larger  than  that  for  either  normal  incidence  or  at  25  de¬ 
grees.  Comparisons  of  the  overall  densities  of  the  films  show  that 
the  largest  density  occurs  for  the  film  deposited  at  normal  incidence, 
and  that  this  density  decreases  for  increasing  angles  of  deposition. 

The  amount  of  decrease  from  normal  incidence  to  50  degrees  is  approx¬ 
imately  10  percent.  The  angle  of  growth  by  the  control  films  is  in 
very  good  agreement  with  what  the  tangent  rule  predicts.  For  example, 
at  50  degrees  the  angle  of  growth  by  the  film  was  found  to  be  approx¬ 
imately  29  degrees,  where  the  tangent  rule  predicts  30.7  degrees. 

All  of  the  above  results  correspond  very  well  to  results  that  other 
models  have  produced.  Therefore,  it  can  be  assumed  that  the  new  model 
is  not  predicting  anything  radical  and  that  other  results  presented  in 
this  chapter  should  be  considered  as  reliable. 


Angle  Variations 

Films  produced  with  angle  variations  showed  no  consistently  repro¬ 
ducible  results.  However,  there  were  some  characteristics  that  could  be 
gleaned  from  the  resulting  films.  First  of  all,  films  produced  with  an¬ 
gle  variations  were  able  to  produce  columns  that  were  on  the  average  a 
little  wider  than  those  with  no  variations,  but  the  spacing  between  the 


columns  increased,  which  resulted  in  decreased  density.  The  most  inter¬ 
esting  results  came  when  angle  variations  were  carried  out  with  a  reason¬ 
able  amount  of  extra  mobility  on  the  part  of  the  disks.  As  shown  in  fig¬ 
ure  27,  a  multilayered  film  with  disk  mobility  occurring  on  5  percent  of 
the  disks,  the  dendrites  formed  are  wider  than  normal  but  the  columnar 
structure  is  fairly  random.  When  this  mobility  is  brought  together  with 
angle  variation  as  in  figure  28,  large  columns  are  formed  from  smaller 
dendrites  which  begins  to  model  actual  films.  In  this  figure,  9  percent 
of  the  disks  have  undergone  mobility  and  are  deposited  with  5  degrees  of 
angle  variation. 

Mobility 

One  interesting  result  obtained  during  this  research  occurred  when 
extra  mobility  was  given  to  a  percentage  of  the  incident  disks  (figures 
27,28,29,30,31).  As  noted  above,  larger  columns  were  formed  when  the  ex¬ 
tra  mobility  of  disks  were  introduced  with  angle  variations.  However, 
this  alone  was  not  the  only  significant  result. 

As  the  percentage  of  disks  undergoing  extra  mobility  was  increased, 
three  distinct  phenomenon  were  found  to  exist.  First,  the  packing  den¬ 
sity  of  the  films  was  increased  (figure  29).  Although  this  comes  as  no 
surprise,  the  amount  of  increase  for  films  deposited  at  oblique  inci¬ 
dence  was  fairly  substantial  when  compared  to  films  deposited  at  normal 
incidence.  In  figure  30,  a  film  deposited  at  40  degrees  with  20.1  per¬ 
cent  mobility,  the  total  density  was  found  to  be  68.1  percent.  When 
this  is  compared  to  the  control  at  50  degrees  (figure  25)  with  a  packing 
density  of  56.3  percent  this  corresponds  to  a  11.8  percent  increase. 

The  film  in  figure  31  (deposited  at  normal  incidence),  with  a  mobility 


of  21.6  percent,  only  showed  a  0.8  percent  increase  in  the  packing  den- 
density,  compared  to  the  control  deposited  at  the  same  angle  of  deposi¬ 
tion.  This  would  seem  to  imply  that  mobility  has  little,  if  any,  ef¬ 
fect  for  films  deposited  at  normal  incidence,  as  far  as  increasing  pack¬ 
ing  density  is  concerned. 

The  second  phenomenon  that  occurred  as  a  result  of  disk  mobility 
was  the  disappearance  of  a  visible  columnar  structure.  In  figures  27,28 
and  29  where  the  mobility  of  disks  was  kept  between  5  and  10  percent,  a 
visible  columnar  structure  is  still  maintained.  However,  when  the  mobil 
ity  of  the  disks  is  increased  to  around  20  percent  and  greater  as  in 
figures  30  and  31,  the  columnar  structure  starts  to  break  down  and  event 
ually  goes  away.  This  disappearance  is  in  agreement  with  what  Nakahara 
had  found  for  films  with  large  structural  relaxation,  such  as  in  copper 
and  permalloy  films  (35). 

The  third  phenomenon  found  in  films  deposited  with  mobility  was 
that  the  films  possessed  a  more  uniform  packing  density  throughout  the 
thickness  of  the  film.  A  density  analysis  of  the  control  film  at  normal 
incidence  showed  that  the  density  of  the  film  tends  to  decrease  the  far¬ 
ther  away  one  gets  from  the  substrate.  When  the  same  analysis  was  done 
on  the  film  in  figure  31,  the  packing  density  did  not  exhibit  a  consis¬ 
tent  drop  and  the  density  near  the  top  of  the  film  was  significantly 
greater  than  that  for  no  mobility.  This  was  found  to  be  true  for  all 
films  deposited  with  extra  mobility.  Although  no  major  study  has  been 
made  concerning  the  uniformity  of  density  throughout  the  thickness  of  a 
film,  the  results  would  confirm  or  deny  that  shadowing  plays  a  role  in 
density  variations  of  a  film,  if  they  exist. 
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Nodular  Defects 

As  stated  earlier  in  chapter  three,  there  are  basically  two  differ¬ 
ent  ways  that  people  have  tried  to  produce  nodular  defects  in  the  film; 
namely,  through  the  use  of  substrate  imperfections  and  the  doping  of  im¬ 
purities  into  the  microstructure.  Throughout  the  course  of  this  re¬ 
search  both  of  these  methods  were  tried,  the  results  of  which  varied. 

Substrate  Imperfections.  As  can  be  seen  back  in  figures  15-20,  the 
substrate  subroutine  had  the  ability  to  produce  almost  any  number  of  sub 
strates.  Some  of  these  were  used  in  the  deposition  of  films,  namely 
those  in  figures  27,29,30,32  and  33.  Although  no  nodules  were  formed 
as  a  result  of  the  imperfections,  some  were  able  to  produce  an  amor¬ 
phous  clump  of  disks  in  the  immediate  vicinity  around  the  imperfection. 
This  can  be  seen  in  figures  29,30  and  33.  The  reason  for  the  amorphous 
nature  of  the  film  around  the  substrate  imperfections  can  only  be  attri¬ 
buted  to  the  irregularity  of  the  imperfections  themselves.  As  to  why 
the  substrates  were  not  able  to  produce  nodules,  it  seems  to  be  due  to 
the  fact  that  the  imperfections  tend  to  be  too  small  in  comparison  to 
the  size  of  the  disks  used. 

Impurities.  The  use  of  impurities  to  produce  nodules  was  by  far 
more  effective  than  using  substrate  imperfections.  As  can  be  seen  in 
figures  30  and  31  (especially  in  figure  30),  distinct  nodules  can  be 
seen  within  the  mircostructure.  Other  attempts  to  produce  nodules  by 
the  use  of  impurities  can  be  seen  in  figures  27,28,29,32,33  and  34. 
Although  these  did  not  produce  nodules,  they  did  produce  an  amorphous 
structure  immediately  in  their  vicinity,  which  can  be  attributed  to  the 
irregularity  of  the  structure,  just  as  in  the  case  of  the  substrate  im- 


perfections.  Another  characteristic  common  to  the  film  growth  around 
the  impurities  is  the  fact  that  the  nodule  or  amorphous  structure  grows 
in  a  direction  parallel  with  the  angle  of  deposition.  For  example,  in 
figure  30  the  nodule  at  the  top  of  the  film  grew  at  an  angle  of  40  de¬ 
grees,  which  happens  to  be  the  angle  of  deposition  for  that  film.  One 
last  characteristic  which  is  common  among  some  of  the  doped  films  is 
that  the  growth  of  the  film  immediately  above  the  impurity  extends  above 
the  top  of  the  film.  This  can  be  seen  in  figures  30,31  and  34. 

It  would  appear  then  that  impurities  can  be  used  to  create  nodules 
and  nodular  effects  such  as  the  production  of  an  amorphous  structure. 

It  is  also  implied  that  one  reason  for  the  creation  of  nodules  in  films 
is  due  to  shadowing,  since  no  physical  parameters  such  as  electrostatic 
charge  are  dealt  with  in  the  model. 


Multilayers 

Perhaps  the  most  interesting  results  obtained  during  this  research 
occurred  when  two  different  size  disks  were  used  to  simulate  multilay¬ 
ered  films  (figures  27,28,32,33,35,36).  In  the  films  produced  when  the 
smaller  disks  were  deposited  first,  as  in  figures  27  and  35,  there  did 
not  appear  to  be  any  outstanding  features  to  note.  The  reason  for  this 
is  that  the  boundaries  between  the  layers  appeared  to  be  continuous  in 
nature.  However,  when  the  larger  disks  were  deposited  first,  two  dis¬ 
tinct  features  were  noticed.  First,  some  of  the  smaller  disks  from  the 
top  layer  were  able  to  penetrate  through  openings  in  the  top  of  the 
lower  layer,  which  made  for  a  discontinuous  boundary.  This  is  espe¬ 
cially  evident  in  figures  32  and  36.  Even  though  one  would  expect  this 
to  occur,  one  reason  for  choosing  different  sized  disks  was  not  neces- 


Multilayered  Film  at  Normal  Incidence 
the  Large  Disks  Deposited  First 


sarily  to  simulate  molecules  of  different  size,  but  to  simulate  mole¬ 
cules  which  possessed  different  interaction  distances.  The  second  fea¬ 
ture  exhibited  by  the  films  when  the  larger  disks  were  deposited  first, 
was  that  this  tended  to  produce  pores  in  the  microstructure  which  could 
allow  for  water  absorption  (see  figures  27,28,32).  As  stated  in  chapter 
two,  one  of  the  problems  with  vapor  deposited  films  is  that  they  are  sub 
ject  to  water  absorption.  Therefore,  these  films  would  tend  to  imply 
that  the  mechanical  properties  of  film  growth,  like  shadowing,  are  one 
reason  why  water  absorption  occurs. 

Other  Characteristics 

One  result  obtained  from  the  research,  that  was  not  expected,  was 
the  fact  that  the  density  of  the  modeled  film  immediately  above  the  sub¬ 
strate  was  significantly  greater  than  the  average  density  of  the  film. 
More  specifically,  the  packing  density  of  the  area  three  to  four  unit 
cells  above  the  substrate  was  found  to  be  25  to  30  percent  greater  than 
any  other  region  in  the  film.  The  main  reason  that  the  density  of  the 
film  near  the  substrate  is  so  much  greater,  is  most  likely  due  to  the 
fact  that  the  shadowing  effect  of  the  film  is  negligible  during  the 
early  stages  of  growth.  Then  as  the  columns  in  the  film  become  larger, 
shadowing  starts  to  play  a  more  important  role  in  the  growth  of  the 
film,  until  the  time  when  the  film  is  high  enough  above  the  substrate 
that  the  density  begins  to  approach  that  of  the  average. 
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VI.  Conclusion  and  Recommendations 


Conclusion 

Based  on  the  results  obtained  in  chapter  five,  it  is  evident  that 
the  objectives  of  the  research  were  met.  That  is,  the  ability  to  pro¬ 
duce  a  computer  simulation  of  thin  film  growth  with  variable  parameters 
was  accomplished. 

Deposition  Angle.  It  was  shown  during  this  thesis  that  the  angle 
of  deposition  plays  an  important  role  in  the  growth  of  a  film.  Not  only 
do  the  width  of  the  individual  columns  making  up  the  mircostructure  de¬ 
crease  with  increasing  angle,  but  the  packing  density  of  the  film  de¬ 
creases  as  well.  The  variation  of  the  angle  of  incidence  also  played  a 
role  in  determining  the  resulting  structure.  For  relatively  small  angle 
variations  (less  than  10  degrees)  the  density  of  the  film  went  up,  but 
for  larger  angle  variations  the  density  went  down  producing  a  rather  a- 
morphous  microstructure. 

Mobility.  By  varying  the  percentage  of  disks  undergoing  extra  mo¬ 
bilities,  three  results  were  obtained.  First,  the  larger  the  percentage 
of  disks  undergoing  extra  mobility,  the  greater  the  packing  density  of 
the  film.  Second,  for  films  with  a  relatively  large  amount  of  disk 
mobility,  the  columnar  structure  all  but  disappeared;  a  result  which  is 
in  agreement  with  experimental  work.  Third,  films  deposited  with  disks 
undergoing  extra  mobility  possessed  a  more  uniform  packing  density. 

Impurities.  Although  there  were  very  few  nodules  formed  by  the  do¬ 
ping  of  impurities,  there  were  two  interesting  results.  First,  the  film 
deposited  immediately  around  the  impurity  was  extremely  more  amorphous 


than  other  sections  of  the  film  which  contained  no  impurities.  Second, 
the  angle  of  the  nodule  and/or  amorphous  growth  was  parallel  with  the 
angle  of  deposition. 

Substrates.  Out  of  the  many  substrates  produced  by  the  model  none 
of  them  produced  nodules  of  any  great  size,  if  at  all.  One  reason  for 
this  may  be  due  to  the  fact  that  the  highest  imperfection  able  to  be 
created  was  only  seven  disks  high.  However,  the  imperfections  in  the 
substrate  did  create  an  unusually  amorphous  structure  about  them,  some¬ 
what  like  the  impurities,  which  could  have  resulted  in  nodules  had  the 
film  thickness  not  been  so  great. 

Multilayers.  The  multilayered  films  produced  by  this  computer  sim¬ 
ulation  were  the  most  interesting  of  all  the  films  produced.  It  is  ap¬ 
parent  by  looking  at  the  films  that  the  boundary  layer  between  the  dif¬ 
ferent  sized  disks  appeared  to  be  discontinuous  in  nature,  an  outcome 
which  is  in  complete  agreement  with  experimental  results.  This  was 
especially  true  of  the  cases  in  which  the  large  disks  were  deposited 
first,  for  this  allowed  the  smaller  disks  to  slip  through  holes  left 
by  the  larger  ones.  Other  things  of  interest  found  in  the  multilayered 
films  were  large  cracks  or  pores  which  could  account  for  water  absorp¬ 
tion,  and  some  nodule  growth  due  to  the  discontinuity  of  the  boundary 
layer . 

Other  Characteristics.  It  was  found  that  the  density  of  the  film 
within  a  few  unit  cells  of  the  substrate  was  significantly  higher  than 
the  average  packing  density  (25  to  30  percent).  This  is  thought  to  oc¬ 
cur  as  a  result  of  shadowing  playing  a  minor  role  in  the  early  stages  of 


Ending  Remarks.  Although  more  two-dimensional  computer  models 
using  disk  mechanics  may  be  developed  in  the  future,  the  results  ob¬ 
tained  should  not  be  significantly  different  from  those  presented  in 
this  thesis.  Three-dimensional  models,  on  the  other  hand,  may  provide 
information  about  form  birefringence  and  other  properties  that  can¬ 
not  be  studied  in  two  dimensions.  However,  for  the  most  part  compu¬ 
ter  simulated  growth  using  disk  mechanics  can  only  give  so  much  infor¬ 
mation.  Not  until  the  time  when  physical  properties  such  as  electrical 
charge  and  chemical  bounding  are  routinely  used  to  construct  the  model, 
will  computer  simulations  tell  people  more  about  thin  film  growth. 

Recommendat ions 

In  order  to  improve  upon  this  new  model  of  thin  film  growth  there 
are  several  things  which  can  be  done.  First,  if  computer  space  will  al¬ 
low,  a  simulation  in  three  dimensions  using  disk  mechanics  is  probably 
the  next  logical  step.  Upon  the  completion  of  modeling  a  three  dimen¬ 
sional  film,  the  introduction  of  physical  parameters  into  the  model  will 
ultimately  make  the  model  more  useful  and  realistic.  Some  of  the  para¬ 
meters  which  could  be  introduced  are  the  use  of  chemical  bonds  in  order 
to  produce  films  with  a  crystalline  structure,  the  variation  of  the  sub¬ 
temperature  to  allow  for  different  molecule  energies,  as  well  as  using 
electrical  charge  on  individual  molecules  to  determine  the  microstruc- 


This  appendix  is  divided  into  two  sections.  Its  purpose  is  to  show 
the  mathematical  derivations  used  to  develop  the  code  in  appendix  B. 

The  first  section  gives  the  mathematical  derivations  needed  to  find  the 
incident  disk  coordinates  after  a  collision  has  occurred.  The  equations 
for  the  final  position  of  the  incident  disk  are  derived  in  the  second 
section. 

Incident  Disk  Coordinates  After  Collision 

The  equation  of  the  trajectory  taken  by  the  center  of  the  incident 
disk  is  a  line.  Rewriting  the  slope  intercept  form  of  a  line  by  solving 
for  x  has  the  form: 

X  -  (l/m)*Y  +  a  (10) 


where 


m  »  slope  of  the  line 
a  »  X  intercept 


The  equation  of  a  contact  circle  discussed  in  chapter  two  can  be 
writen  in  the  form: 


r2  -  (X  -  h)2  +  (Y  -  k)2 


(11) 


where 


h  ■  X  coordinate  of  the  center  of  the  contact  circle 
k  «  Y  coordinate  of  the  center  of  the  contact  circle 
r  »  radius  of  the  contact  circle 


In  order  to  solve  for  the  value  of  y,  the  value  of  x  in  the  first 


equation  is  substituted  for  x  in  the  second  equation  giving  the  follow¬ 
ing: 

r2  -  [(Y/m)  +  a  -  h]2  +  [Y  -  k]2  (12) 

Expanding  the  squares  and  putting  the  equation  into  quadratic  form 
gives: 

[(1/m2)  +  1]*Y2  +  2*[(a  -  h)/(m  -  k)]*Y  +  (a  -  h)2  +  k2  -  r2  =  0  (13) 

Solving  the  above  equation  for  y  gives: 

y  ,  -(f(a  -  h)/m1  -  k)  +  (r2»r(l/m2)  +  11  -  f(a  -  h)  +  (k/m)12)1/2  (3) 

[(1/m2)  +  1] 

The  square  root  is  added,  and  not  subtracted  since  the  collision 
occurs  at  the  larger  y  value.  In  order  to  solve  for  the  value  of  x,  the 
value  of  y  is  simply  put  back  into  the  equation  of  the  trajectory.  Be¬ 
cause  the  angle  of  incidence  is  measured  from  the  substrate  normal  the 
value  of  [(1/m  )  +  1]  is  actually  the  secant  squared  of  the  incident  an¬ 
gle,  and  is  represented  by  sc2ang  in  the  program.  The  value  of  [(a-h)  + 
2 

(k/m)]  is  represented  by  g  and  the  value  of  [(a-h)/m  -  k]  by  f. 

Incident  Disk  Rest  Point  Coordinates 

Given  the  equations  of  the  collision  and  rest  disk  contact  circles 
respectively  as 


2  2  2 
rlZ  -  X  +  Y^ 

(14) 

(X  -  h)2  +  (Y  -  k)2 

(15) 

rl  is  the  radius  of  the  collision  disk  contact  circle,  r2  the  radius  of 
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the  rest  disk  contact  circle,  h  the  x  seperation  between  the  two  disks 
and  k  the  y  seperation  betweenthe  two  disks. 

It  is  important  to  note  that  the  equations  above  are  written  for 
circles  which  have  been  translated  to  the  origin.  Solving  for  x  in 
equation  14  gives: 

I  -  +/-  (rl2  -  y2)1/2  (16) 

Substituting  this  value  for  x  into  equation  15  gives: 

r22  -  [(rl2  -  Y2)1/2  -  h]2  +  (Y  -  k)2  (17) 

After  expanding  the  squares  and  cancelling  values,  the  equation 
looks  like  the  following: 

r22  »  rl2  +  h2  -  2h(rl2  -  Y2)1/2  +  k2  -  2ky  (18) 

Through  algerbraic  manipulation  this  can  be  transformed  into  the 
the  following  quadratic  equation: 

(k2  +  h2)*Y2  +  2aky  +  (a2  -  h2*rl2)  -  0  (19) 

2  2  2  2 

where  a  =  (r2  -  rl  -  h  -  k  )/2.  Solving  this  quadratic  equation 

yields  the  two  possible  values  of  y  namely 

Y  -  [-a*k  +/-  (a2*k2  -  c*b)]/c  (5) 

2  2  2  2  2 

where  b  ■  (a  -  h  *rl  )  and  c  »  (k  +  h  ).  The  two  values  of  x  can  be 
solved  for  in  a  likewise  manner  yielding  the  following  equation: 


m 


-A 


x 


[-a*h  +/-  (a2*h2  -  c*b)]/c 


(6) 


Fortran  Program  Code 

The  following  pages  of  this  appendix  contain  the  Fortran  code  used 
to  model  thin  film  growth  by  vapor  deposition.  More  specifically  the 
program  code  is  composed  of  Fortran  77  and  contains  a  fortran  statement 
using  IMSL  (International  Mathematical  and  Statistical  Libraries),  a 
library  of  subroutines  used  for  statistical  purposes.  In  this  program 
the  IMSL  statement  (ggubfs)  is  used  as  a  random  number  generator  for 
evenly  distributing  the  simulated  evaporate  over  the  substrate. 

Questions  concerning  the  use  of  particular  variables  within  the  pro¬ 
gram  can  be  addressed  by  referring  to  appendix  C,  which  contains  a  list 
of  the  variables  along  with  a  brief  discription  of  their  purpose. 
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•THIN  FILM  GROWTH  SIMULATOR  -  AN  IMPROVED  MODEL* 

VERSION  -  2 

WRITTEN  8 Y  OAVIO  J.  OORYLANO 
30  SEP  35 


LIKE  ITS  PREDECESSOR »  THE  PURPOSE  OF  THIS  PROGRAM  IS  TQ  SIMULATE 
THE  VAPOR  DEPOSITION  OP  THIN  FILMS.  THE  PROGRAM  ITSELF  IS  01- 
VIOEO  INTC  4  MAJOR  AREAS  WITH  EACH  CONTAINING  A  38I5P  DESCRIP¬ 
TION  OF  ITS  FUNCTION.  THE  REASON  BEHIND  THIS  IS  THAT  IT  PROVIDES 
NOT  ONLY  AN  EASIER  WAY  TO  CE3UG  THE  PROGRAM  BUT  IT  ALSO  HELPS  OR¬ 
GANIZE  IT.  THIS  -PROGRAM  IS  WRITTEN  IN  FORTRAN  77  ANO  REQUIRES  A 
WORKING  KNOWLEDGE  OF  IMSL  ANO  S  (FOR  FUTURE  PLOTTING).  IT  WAS  DE¬ 
SIGNED  TO  RUN  ON  A  VAX/VMS  COMPUTER. 


*  * 

*  MAIN  PROGRAM  START-UP  * 

*  ANO  INTERFACE  LOOP  * 

*  * 


THE  PURPOSE  OF  THIS  SECTION  IS  TO  OECLARE  ANO  INITIALIZE  THE  MAIN 
VARIABLES  USEO  THROUGHOUT  THE  PROGRAM.  IT  ALSO  INCLUDES  A  SIGN 
ON  MESSAGE,  A  HELP  MESSAGE,  THE  READING  OF  THE  DATE  AND  THE  MAIN 
INTERFACE  LOOP. 


aaaasaaaa*  VARIABLE  DECLARATION  aaaaaaPaaa 

REAL  OATE,  XCOCRQ(300,200),  YCOCROC300 ,200) 

INTEGER  FILMC300,200)  ,  IVAR 


aaaaaaaaaa  SIGN  ON  aaaaaaaaaa 


WRITE  (6,0C10) 

0010  FORMAT  (10(2X/),15X, 'THIN  FILM  GROWTH  SIMULATOR  -  AN 
CL*/, 32X, 'VERSION  -  2 '/, 3 3X ,•»»»**■■*«■•/, 25 X  ,' WR ITTS N 
COORYLANO'/,33X,'30  SEP  3 5 ' / , 2 ( 2 X/ )  ,  1 9X , • THIN  FILM  RES 
CAT  WAY  OF  LIPE'/,7(2X/)) 

C 

C 

C  aaaaawaaaa  READING  OF  THE  DATE  «aaa»**a*a 

C 


IMPROVED  MOOE 
BY  OAVIO  J. 
EARCH  -  A  GRE 


WRITS  (6,0015) 

0015  FORMAT  (IX, 'ENTER  CATE  ANO  TIMS  (MMQOHHMM)  • ) 

REAO  (5,0016)  OATS- 
0016  FORMAT  (1=3.0) 

C 

C 

C  aaaaaaaoaa  HELP  MESSAGE  «a»aaaaaa« 

C 

WRITE  (6,0017) 

0017  FORMAT  (10(1X/),10X,'TO  OPERATE  THIS  PROGRAM,  ANSWER  THE  3U  E  S  T IONS 
C  ASKED  3  Y  USING'/, 10X, 'TH=  KEYBaO  LCCATSO  CN  T-E  RIGHT  SIDE  OF  Twl 
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cs  Terminal,  at  the*/, iox, 'conclusion  of  each  numbered  section,  pr 
CESS  THE  CARRIAGE  RETURN •/, 10X X E Y  TO  CONTINUE  WITH  THE  PROGRAM. 
CIN  ORDER  TC  EXIT  THE  PROGR AM  •  / , 1 0 X , ' AT  ANY  TIME,  PRESS  CONTROL  C  ( 
CCTRL/C).  TO  ENTER  THE  PROGR AM , • / , 1  OX , • PR  ESS  THE  CARRIAGE  RETURN.' 
C,7(1X/)) 

READ  (5,0021) 


aaaaaaaaa*  MAIN  INTERFACE  LOOP  aaaaaaaaaa 


0019  WRITE  (6,0020) 

0020  FORMAT  ( 15 ( 2X/) , 29X  ,  '  M A  IN  INTERFACE  LOOP '///, 2QX ,' ENT ER  APPRQPRIAT 
CE  NUMBER  TO  CONT INUE ' / ,  2( 2X/ ) , 26 X , • 1  *  SUBSTRATE  SUB  ROUT INE '/, 26X , 
C ' 2  -  DEPOSITION  SUBROUT INE ' / ,  26X ,  •  3  -  ANALYSIS  SUBROUTINE'/, 26X,  '4 
C  -  MOVING  SUBROUTINE'/, 26X, 'CTRL/C  -  EXIT  TO  OPERATING  SYSTEM'/, 6C 
C2X/) ) 


READ  (5,0021)  IVAR 

FORMAT  (ID 

IF  (IVAR  .EQ.  1)  THEN 

CALL  SUB (FILM,XCOORO,YCOORO, HEIGHT) 

GOTO  0019 

ELSE  IF  (IVAR  .EQ.  2)  THEN 

CALL  C£P0(FILM,XC00R0,YC00R0, HEIGHT, DATE) 

GOTO  0019 

ELSE  IF  (IVAR  .EQ.  3)  THEN 

CALL  ANAL(FlLM,XCOORO,YCOORC) 

GOTO  0019 

ELSE  IF  (IVAR  .EQ.  4)  THEN 

CALL  MCVE(FILM,XCOORO,YCOORO) 

GOTO  0019 

ENCIF 

WRITE  (6,0030) 

FORMAT  ( 20X ,' INCORRECT  RESPONSE  OUMMYN  -  TRY  AGAIN'///) 

GOTO  0031 


aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
a  a 
*  SUBSTRATE  SUBROUTINE  a 
a  a 
aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


this  subroutine  CONSIST  of  5  areas:  that  is  THE  initializing,  pro 
OUCING,  STORING,  TRANSFERING  ANC  RECALLING  QF  SUBSTRATES.  EACH 
SECTION  CONTAINS  A  BRIEF  DESCRIPTION  OF  ITS  PURPOSE  ANC/OR  ITS  AP 
PRCACH. 


SUBROUTINE  SUB  ( FILM,  XCOORO, YC009D, HEIGHT) 

LOCAL  VARIABLES 

real  XCOLO,  HGT,  XPOS,  YCCLD,  HEIGHT,  YC(300,10),  XC(  300  ,  10  ) 

REAL  XCOORCC  300, 200  ,  YC 00 R 0 (  300 , 2 00  ) 

INTEGER  SUEBUFUOO.IO),  COL,  XC3L ,  ANS,  IVAR,  FILM(  300 , 200) ,  YHGT 


A-'Ij*  L: 
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C  **$*3$$*$#  SUBSTRATE  INTERFACE  LOOP  4444444444 

C 

0040  WRITE  <6,0041) 

0041  FORMAT  C15C2X/),25X, 'SUBSTRATE  INTERFACE  LOOP ' /// , 20 X , • ENTE R  APPRO 
CPRIATE  NUMBER  TO  CONTINUE ' / , 2 C 2 X/ )  ,  2 6X  ,  •  1  -  INITIALIZE  SUBSTRATE*/ 
C  ,  26X  ,  1  2  -  CREATE  SUBSTRATE */,  26X  3  -  STORE  SUB  ST R AT E • / , 26X , • 4  -  R 

CECALl  SUBSTRATE* /,26X, *5  -  MOVE  SUBSTRATE  ( TQ  FILM  MATRIX) '/ ,26X  ,  • 
C 9  -  EXIT  TO  MAIN  INTERFACE  LOOP • / , 26X ,  •  C TRL/C  -  EXIT  TO  OPERATING 
CSYSTEM'/,5C2X/)) 

C 

0049  READ  (5,0050)  IVAR 
0050  FORMAT  <I1) 

C 

IF  (IVAR  .EQ.  1)  THEN 
GOTO  0060 

ELSE  IF  (IVAR  .EQ.  2)  THEN 
GOTO  OCTO 

ELSE  IF  (IVAR  .EQ.  3)  THEN 
GOTO  0030 

ELSE  IF  (IVAR  .EQ.  4)  THEN 
GOTO  0090 

ELSE  IF  (IVAR  .EQ.  5)  THEN 
GOTO  0100 

ELSE  IF  (IVAR  .EQ.  9)  THEN 
RETURN 
ENOIF 

WRITE  (6,0051) 

0051  FORMAT  ( 21X , *  YOU  BLEW  IT  3L0CXHEAD 1  I  -  TRY  AGAIN'///) 

GOTO  0049 


4444444444  SUBSTRATE  INITIALIZATION  4*44444444 

THIS  PART  OF  THE  SUBROUTINE  INITIALIZES  SUB9UF,  XC,  YC  TO  ZERO  SO  THAT 
A  NEW  SUBSTRATE  CAN  BE  CREATED.  UPON  ENTERING  THIS  SECTION  THE  INTIAL- 
IZATION  DOES  NCT  AFFECT  THE  FILM  MATRICES  ( F I LM , XCOGR D  ,  YCOC RO  )  . 

0060  WRITE  (6,0064) 

0064  FORMAT  (lOX.'UPON  ENTERING  YOU  MAY  LOSE  THE  CONTENTS  OF  THE  SUBROU 
CTINc  BUFFER !’/, 10X, ' 00  YOU  WISH  TO  SAVE  THE  OLO  CONTENT S ?'//,  1 1 X  ,  ' 
CIr  YES  -  ENTER  1 */,  1 IX ,' OTHERWISE  -  HIT  RETURN'//) 

C 

READ  (5,0050)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  OOBO 
ENOIF 
C 

00  0061  X  a  1,300,1 
OC  0062  Y  a  1,10,1 
SUB3Ue( X  ,  Y  )  »  0 
XCCX.Y)  a  o.o 
YC ( X , Y )  a  0.0 
0062  CONTINUE 
0061  CONTINUE 
C 

WRITE  (6,0063) 

0063  FORMATdOX,' ALL  VALUES  OP  TmE  SUBROUTINE  BUFFER  HAVE  BEEN  SET  TO  Z 
CERO  '  //) 

*»*  HIT  RETURN' 

RSAO  (5,0050) 
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GOTO  0040 


****4#««4#  SUBSTRATE  CREATION  44444444*4 

THIS  PART  OF  THE  SUBROUTINE  ALLOWS  THE  PROGRAMMER  TO  CREATE  A  SUBSTRATE. 

THE  FIRST  STEP  IS  TO  LAY  OOWN  A  SINGLE  LATER  OP  DISK.  AFTER  THE  SINGLE 
LAYER  HAS  BEEN  LAYEO  OOWN,  THEN  THIS  PART  OF  THE  SUBROUTINE  HELPS  IN  CREAT¬ 
ING  THE  REMAINCER  OF  THE  SUBSTRATE.  THIS  IS  DONE  BY  NOTING  THAT  EACH 
DISK  IS  THE  SQUARE  ROOT  OF  TWO  IN  OIAMETER  AND  THAT  THESE  CAN  BE  STACKED 
ONE  UPON  ANOTHER. 

0070  WRITE  (6,00720 

0072  FORMAT  (IOX, 'JUST  A  REMINDER.  UPON  ENTERING  INTO  THIS  PART  OF  THE 
C'/, IOX, • SUBROUTINE  YOU  MAY  BE  DESTROYING  A  PREVIOUS  SUBSTRATE.  1 
COX, *00  YOU  WISH  TO  SAVE  THE  OLO  SUBSTRATE? '//,  UX  ,' IF  YES  -  ENTER 
C  1'/,11X, 'OTHERWISE  -  HIT  RETURN'//) 

C 

READ  (5,0050)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  0080 
ENOIF 
C 

C  CREATES  A  HORIZONTAL  SUBSTRATE  WITH  NO  DEFECTS 
C 

XC(l.l)  *  1.0 
YCC1.1)  *  1.0 
SU33UFC1.1)  *  1 
XCCLO  «  1.0 

0071  IF  (XCOLO  .LT.-300)  THEN 

XCCLO  »  XCOLO  ♦  1.4142136 
X  a  AINT (XCOLO ) 

XCCX.l)  *  XCOLO 
YC ( X  ,  1)  =  1.0 
SUBBUF(X.l)  ■  1 
GOTO  0071 
SNCIF 

HEIGHT  *  1.0 
C 

WRITE  (6,0073) 

0073  FORMAT  (5X,'00  YOU  WISH  TO  CREATE  MORE  THAN  JUST  A  HORIZONTAL  LINE 
C  FOR  THE  SUBSTRATE?' //,11X, ' IF  YES  -  ENTER  1  •/,  1 IX ,' OTHERWISE  -  HI 
CT  RETURN'//) 

C 

RSAO  (5,0050)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  0074 

ELSE 

GOTO  0040 
ENOIF 
C 

0074  WRITE  (6,1073) 

1073  format  (iox, 'there  are  213  columns  across  the  bcttcm  of  the  substr 

CATE. '/, 10X, ' EACH  ONE  OF  THESE  MUST  HAVE  A  HEIGHT  SPECIFIED  TO  IT.' 

C/« IOX , ' THE  FOLLOWING  STEPS  WILL  HELP  IN  THAT  PROCESS.  '//) 

C 

COL  *  0 
HEIGHT  *  0 

1074  IF  (COL  .EC.  213)  THEN 

GOTO  0079 
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ELSE 

COL  *  COL  ♦  1 
WRITE  (6,0075)  CCL 

0075  FORMAT  (lOX.'YOU  ARE  IN  COLUMN' ,I4/,10X, 'ENTER  HEIGHT  (MUST  BE 

♦  7  OR  LESS)'//) 

C 

1080  RSAO  <5,0050)  YHGT 

IF  (YHGT  .EQ.  1)  THEN 
GOTO  0076 

ELSE  IF  CYHGT  .GT.  7)  THEN 
WRITE  (6,1075) 

1075  FORMAT  (10X, 'HEIGHT  TOO  HIGH  -  TRY  AGAIN l • //) 

GOTO  1030 

ELSE 

C 

C  CREATES  THE  HEIGHT  SPECIFIED  FOR  AN  INDIVIDUAL  COLUMN 
C 

YCCLD  a  1.0 

XPCS  *  1.0  ♦  (COL  -  1)*1. 4142136 
X  =  AINT(XPOS) 

HGT  *  YHGT  -  1 
00  0077  I  »  1 , HGT , 1 

YCOLD  =  Y COLO  ♦  1.4142136 
IF  CYCOLO  .GT.  HEIGHT)  THEN 
HEIGHT  »  YCOLD 
SNOIF 

Y  *  AINT(YCOLO) 

XC(X,Y)  *  XPOS 
YC ( X  ,  Y )  a  YCQLO 
SUBBUFCX.Y)  a  1 
0077  CONTINUE 

C 

0076  WRITE  (6,0078)  YHGT 

0078  FORMAT  (lOX.'OO  YOU  WISH  THIS  HEIGHT  (',11,')  TO  CONTINUE? 

C'//,11X,'IF  YES  -  ENTER  1'/,11X» 'OTHERWISE  -  HIT  RETURN'//) 

C 

REAO  (5,0050)  ANS 
IF  (ANS  .EQ.  1)  THEN 
WRITE  (6,1070) 

1070  FORMAT  ( 10X ,  'TOO  WHAT  COLUMN  SHOULD  THIS  CONTINUE?*/ *1 

C8X, ' (MUST  BE  213  OR  LESS)'//) 

1077  READ  (5,1050)  XCOL 

1050  FORMAT  (13) 

IF  (XCOL  .GT.  213)  THEN 
WRITE  (6,1078) 

1078  FORMAT  (10X, 'COLUMN  NUMBER  TOO  LARGE  -  TRY  AGAIN'//) 

GOTO  1077 

SNOIF 

FILLS  IN  THE  SU3STRATE  MATRICES  TO  THE  COLUMN  SPECIFIEC  WITH  THE  LAST 
ENTERED  HEIGHT. 


COL  »  CCL  ♦  1 

IF  (YHGT  . EQ.  1)  GOTO  1079 
00  1071  I  »  COL , XCCL ,  1 

XPOS  a  1.0  ♦  (I  -  1)#1. 4142136 
X  *  AINT(XPOS) 

YCQLO  »  1.0 
00  1072  J  *  l.HGT.l 

YCCLO  «  YCQLO  ♦  1.4142136 


•  in* 


Y  »  AINT( YCCLQ) 

XC(X,Y)  «  XP0S 
YC(X,Y)  *  YCOLD 
SUB3UF(X, Y)  »  1 
1072  CONTINUE 

1071  CONTINUE 

1079  COL  »  XCCL 

GOTO  1074 

ELSE 

GOTO  1074 
ENOIF 
ENOIF 
ENOIF 

0079  WRITE  (6,1076) 

1076  FORMAT  (lOX.'THE  SUBSTRATE  IS  NOW  COMPLETE !!'/, 10X ,' R E AD Y  FOR  TRAN 
CSFER  TO  FILM  MATRICES.'///) 


CSFER  TO  FILM  I 
PRINT  *, ' 

R c AO  (5,0050) 
GOTO  0040 


HIT  RETURN' 


SUBSTRATE  STORAGE 


THIS  PART  Or  THE  SUBROUTINE  ALLOWS  ONE  TO  TEMPORARILY  STORE  THc  CONTENTS 
OF  SU33UF,  XC  AND  YC  IN  THE  file  CALLED  FOROH.OAT.  IT  DOES  NOT  DESTROY 
CONTENTS  OF  THESE  MATRICES,  aUT  MERELY  COPIES  THEM  FOR  FUTURE  REFERENCE. 

0080  WRITE  (6,0081) 

0081  FORMAT  CIOX.'UPON  ENTERING  THIS  PART  OF  THE  SUBROUTINE,  THE  SUBROU 
CTINS'/,10X, 'BUFFER  WILL  BE  STGRSO  IN  A  FILE  CALLED  FOROU.CAT.  »/,  1 
COX, 'IS  THIS  YOUR  DESIRE?'//, 11X, 'IF  YES  -  ENTER  1'/,11X, 'OTHERWISE 
C  -  HIT  RETURN'//) 

READ  (5,0050)  TVAR 
IF  (IVAR  .EO.  1)  THEN 
PRINT  *,'  ' 

PRINT  6,  'WORKING' 

GOTO  0088 

ELSE 

GOTO  0040 
ENCIF 


0088  OPEN  (UNIT  a  U,FILE  »  '  FOR  0  1 1 .  C  A  T  '  ,  S  T  ATUS  *  'NEW') 

OQ  0082  X  *  1,300,1 
00  0083  Y  *  1,10,1 

WRITE  (11,0084)  SUB9UF(X,Y),  XC(X,Y),  YC(X,Y> 
0094  FORMAT  (  1  X  ,  I  2 , 3 X ,F3 . 4 , 3X , F 3 . 4 ) 

0033  CONTINUE 
0082  CONTINUE 

CLOSE  (UNIT  =  10) 


WRITE  (6,0036) 

FORMAT  (  1 0  X ,  ' THE  VALUES  OF  THE  SUBROUTINE  BUFFER  HAVE  NOW  ! 
10X, 'COPIED  INTO  THE  PILE  FORO ll.OAT' //) 

PRINT  *, •  HIT 

ASAO  (5,0050) 

GOTO  0040 


RETURN' 
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C 

c  *#$:**$**6*  SUBSTRATE  RECALLING  $*$*$*#*#* 

C 

C  THIS  PART  OF  THE  SUBROUTINE  ALLOWS  ONE  TO  ACCESS  ANO  RECALL  THE  CQN- 
C  TENTS  OF  A  FILE  CALLEO  TAKE. OAT.  IN  DOING  SO,  THE  VALUES  QF  A  PRE- 
C  VICUSLY  OEVELOPEO  SUBSTRATE  WILL  BE  ENTERED  INTO  THE  MATRICES  SUBBUF, 

C  XC  ANO  YC. 

C 

0090  WRITE  (6,0091) 

0091  FORMAT  (10X,'JUST  A  REMINDER.  UPON  ENTERING  THIS  PART  OF  THE  SUBR 
COUTINE'/.IOX.'THE  CURRENT  VALUES  OF  THE  SUBROUTINE  SUFFER  MAY  BE  L 
COST. '/,10X, 'THE  REASON  FOR  THIS  IS  THAT  NEW  VALUES  FROM  TAKE. OAT  W 
CILL  BE  REAC  IN. » / , 10X , • 00  YOU  WISH  TO  SAVE  THE  CURRENT  VALUES  OF  T 
CHE  BUFFER?*//, 11X, 'IF  YES  -  ENTER  1 •/,  1 1 X OTHERWISE  TO  CONTINUE  - 
C  HIT  RETURN'//) 

C 

READ  (5,0050)  IVAR 
IF  (IVAR  . EQ.  1)  THEN 
GOTO  0080 
ENOIF 

PRINT  *, 'WORKING' 

C 

OPEN  (UNIT  »  11, FILE  »  • TAKE .OAT ', STATUS  *  'CLO') 

REWIND  (UNIT  *  11) 

C 

00  0095  X  *  1,300,1 
DO  0096  Y  a  1,10,1 

READ  (11,0097)  SUBBUF ( X , Y ) ,  XC(X,Y),  YC(X,Y) 

0097  FORMAT  ( 1 X  ,  12 , 3X , F8 . 4 , 3X , F8 . 4) 

0096  CONTINUE 
0095  CONTINUE 
C 

CLOSE  (UNIT  *  11) 

C 

HEIGHT  a  10 
C 

WRITE  (6,0098) 

00 98  FORMAT  (lOX.'THE  NEW  VALUES  OF  THE  SUBROUTINE  BUFFER  HAVE  NOW  BEEN 
C'/.lOX,' REAO  IN  FROM  THE  FILS  TAKE. OAT'//) 

PRINT  *,'  HIT  RETURN' 

REAO  (5,0050) 

GOTO  0040 
C 
C 

c  SU3STRATS  MOVING  as******** 

C 

C  THIS  PART  Or  THE  SUBROUTINE  COPIES  THE  CURRENT  VALUES  OF  SUBBUF,  XC  ANC 
C  YC  INTO  THE  MATRICES  FILM,  XCOCRO  ANO  YCOORO  RESPECTIVELY.  THE  REASON 
C  FOR  THIS,  IS  SC  THAT  THE  DEPOSITION  PROCESS  CAN  BE  INVCKEO  WITH  A  SUB- 

c  strate  in  place. 
c 

0100  WRITE  (6,0101) 

0101  FORMAT  (  1  OX , ' UPON  ENTERING  THIS  PART  OF  THE  SUBROUTINE  THE  CURRENT 
C  VALUES '/, 10X , 'OF  THE  SUBROUTINE  BUFFER  WILL  BE  COPIED  INTO  Th£  FI 
CLM  MATRICES. '/,10X, 'IS  THIS  TOUR  CES IRE ? ' // , 1 1 X , • IF  YES  -  ENTER  1 
C'/,11X, 'OTHERWISE  -  HIT  RETURN'//) 

C 

REAO  (5,0050)  IVAR 
IF  (IVAR  .EO.  1)  THEN 
GOTO  0103 


■'.-X 
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'  y.vr,  7  V  y  ? f_*  7  !A  ^  w* »  ■  i  ■  l  *  i  ■  h  p  ■  g  ■  l"p  m  w  n 


ELSE 

GOTO  0040 
SNOIF 
C 

0103  00  0104  X  «  1,300,1 

00  0105  Y  «  1,10,1 

XCQQROC  X  » Y  3  «  XC ( X  ,  Y  ) 

YCQCROCX,  Y)  »  YC<  X  » Y  3 
PILM<X,Y3  *  SU83UFCX  ,  Y3 
0105  CONTINUE 
0104  CONTINUE 
C 

WRITE  <6,01063 

0106  FORMAT  (lOX.'THE  SUBROUTINE  VALUES  HAVE  BEEN  TRANSFEREO ' //) 

WRITE  <6,01073 

0107  FORMAT  <10X,'00  YOU  WISH  TO  RETURN  TO  THE  MAIN  INTERFACE  LOOP?'//, 
CllX.'IF  YES  -  ENTER  1 ' /* 11X ,  'OTHERWISE  -  HIT  RETURN'//3 
C 

READ  <5,00503  IVAR 
IF  < IVAR  .EQ.  13  THEN 
RETURN 

ELSE 

GOTO  0040 
ENOIF 
ENO 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 


*  * 

*  DEPOSITION  SUBROUTINE  * 

*  * 

THIS  SUBROUTINE  CONSIST  OF  3  AREAS:  THAT  IS  THE  INITIALIZING  OF 
THE  FILM  MATRICES,  THE  SETTING  OF  ThE  DEPOSITION  VARIABLES  AND 
THE  DEPOSITION  OF  THE  FILM  ON  A  SUBSTRATE.  EACH  SECTION  CONTAINS 
A  BRIEF  DESCRIPTION  OF  ITS  PURPOSE  ANO/QR  ITS  APPROACH. 

SUBROUTINE  OEPO< FILM , XCOORO , YCOORO , HE IGHT , C ATE  3 

LOCAL  VARIABLES 

REAL  DATE 

REAL  ANG,  ANG VAR, RANG, OLOANG, ANGINC.HE IGHT, R ANNUM, YTEST.YBNC <63 
REAL  XC30RC< 30 0,200 3, YCOORO <300, 2 00 3 ,  H , K , A , 3 ,C , R3NC < 63 ,X3NC<63 
DOUBLE  PRECISION  OSE ED  ,  R E STY  ,  R ES TX  ,  R 1 , R 2 , OR E STX 
DOUBLE  PRECISION  Te ST1 , T E ST  2 ,T E ST  3 , T E ST4 , TS S T Y , OR EST Y 
INTEGER  PILM<300,2003,IVAR,NOOSA,BNCVAR,MGTLAY,LAYER,SIZE,BNCCTR 
INTEGER  YSUB,  YFILM.nO EP.OSPCTR, XCCLL.Y COLL, HCTR, IMP X.3CTR 

«»««*«»» DEPOSITION  INTERFACE  LOOP  «*«««#*## 

0110  WRITE  <6,01113 

0111  FORMAT  <15<2X/3,26X, ’DEPOSITION  INTERFACE  LOOP • /// , 2 0 X , • ENTE R  APPR 
C3PRIATE  NUMBER  TO  CONT INUE ' / , 2< 2 X/ )  , 2 5 X , • 1  -  INITIALIZE  FILM  M  A  T  R I 


OPERATING  SYSTEM'/, 6<2X/3) 


C 

C 

C 

C 
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0114 

REAO  (5,0115) 

IVAR 

0115 

FORMAT  (ID 

L 

IF  ( IVAR  .EQ. 

1)  THEN 

GOTO  0120 

ELSE  IF  (IVAR 

.EQ. 

2) 

THEN 

GOTO  0130 

ELSE  IF  (IVAR 

.EQ. 

3) 

THEN 

GOTO  0160 

ELSE  IF  (IVAR 

.EQ. 

9) 

THEN 

RETURN 

ENOIF 

WRITE  (6,0116) 

0116  FORMAT  C23X, 'WRONG  SCREWBALL!!  -  TRY  AGAIN'///) 

GOTO  0114 
C 
C 

C  *$$*#$*$$#  DEPOSITION  INITIALIZATION  a#**#**### 

c 

C  THIS  PART  OF  THE  SUBROUTINE  INITIALIZES  ALL  THE  VALUES  OF  FILM,  XC00R0 
C  ANO  YC00R0  -  FROM  T  *  11  TO  r  *  200  -  TO  ZERO.  THIS  IS  DONE  TO  ENSURE 
C  THAT  ANY  PREVIOUS  DEPOSITION  TRIALS  ARE  CLEARED  OUT. 

C 

0120  WRITE  (6,0121) 

0121  FORMAT  ( 10X , ' UPON  ENTERING  THIS  SECTION,  ALL  VALUES  OF  THE  FILM  MA 
CTRICES'/,10X,'WILL  BE  0ESTR0YE0.  IN  ORDER  TO  SAVE  THESE  VALUES  YO 
CU  MUST'/, 10X, 'EXIT  THIS  SUSROUT INE . • / , 10 X , • 00  YOU  WISH  TO  SAVE  THE 
CSE  VALUES?'//, 11X, 'IF  YES  -  ENTER  1 •/, 11 X ,' OTHERWISE  -  HIT  RETURN* 
C//) 

c 

REAO  (5,0115)  IV AR 
IF  ( I V A R  .SQ.  1)  THEN 
RETURN 

SNOIF 

PRINT  *, 'WORKING' 

C 

00  0125  X  »  1,300,1 

00  0126  Y  «  11,200,1 
FILM(X.Y)  *  0 
XCCORO ( X  ,  Y )  »  0.0 
YCOORO(X,Y)  »  0.0 
0126  CONTINUE 
0125  CONTINUE 
C 

WRITE  (6,0127) 

0127  FQfiMAT(2(lX/),10X,  'THE  VALUES  CF  THE  FILM  MATRICES  HAVE  BEEN  SET  T 
CO  ZERO. '/,10X, 'DEPOSITION  VARIABLES  ARE  READY  TO  BE  SET!'//) 

PRINT  HIT  RETURN* 

REAO  (5,0115) 

GOTO  0110 


READING  OF  DEPOSITION  VARIABLES  *n****innn> 

THIS  PART  OF  THE  SUBROUTINE  RSACS  IN  SUCH  VARIA3LS  AS:  THE  NUMBER  OF 
DISKS  TO  BE  0EP0SITS0,  THE  ANGLE  OF  INCIDENCE,  ANY  VARIATION  IN  THE 
ANGUS  CF  INCIDENCE,  WHETHER  CR  NOT  BOUNCING  WILL  OCCUR,  WHETHER  AN  IM¬ 
PURITY  IS  TO  BE  DEPOSITED  AND  Th£  SIZE  OF  THE  DISKS  T 3  BE  USED  IN  CON¬ 
JUNCTION  with  the  depth  of  the  laters.  this  is  accomplished  by  asking 


Si 


►  ’^3 


78 


C  A  NUMBER  OF  QUESTIONS. 

C 

0130  WRITS  (6,0131) 

0131  FORMAT  C10X,'3EF0RE  GOING  FARTHER,  IT  IS  ASSUMED  THAT  A  SUBSTRATE 
CHAS ’ / , 1 0 X , 'BEEN  PfiCOUCSO  OR  RECALLED  AND  THEN  TRANSFEREE  TO  THE  FI 
CLM'/.lOX, 'MATRICES.  IT  IS  ALSO  ASSUMED  THAT  THE  TCP  PORTION  OF  TH 
CS '/« 10X, ' FILM  MATRICES  HAVE  BEEN  INITIALIZED. '/,1QX»'IS  THIS  CORRE 
CCT? ' // , 1 IX , ' IP  YES  -  ENTER  1'/,11X,  'OTHERWISE  -  HIT  RETURN*//) 

C 

REAO  (5,0115)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  0132 

ELSE 

RETURN 

ENOIF 

C 

C  INITIALIZATION  OF  QUESTION  VARIABLES 
C 

•  0132  NQOSK  *  0 

OLCANG  a  0.0 
ANGVAR  a  0.0 
ANGINC  a  0 
BNCVAR  a  0 
NOMOB  «  0 
IMPUR  *  0 
LAYER  *  0 
HGTLAY  a  0 
SIZE  a  X 
C 

C  QUESTION  PERTAINING  TO  THE  NUMBER  OF  DISK 
C 

0133  WRITE  (6,0134) 

0134  FORMAT ( 10X , ' HOW  MANY  OISKS  00  YOU  WISH  DEPOSITED? '/, 10X  ,' ENTER  VAL 
CUE  IN  INTEGER  FORMAT.'//) 

REAO  # , NOCSK 

WRITE  (6,0135)  NOOSK 

0135  FORMAT  ( 2( IX/) , 10X , 1 5 , '  DISK  WILL  BE  DEPOSITED,  UNLESS  FILM  MATRIC 
CES  BECOME  FULL'//) 

C 

C  QUESTIONS  PERTAINING  TO  THE  ANGLE  OF  DEPOSITION 
C 

WRITS  (6,0136) 

0136  FORMAT  ( 1 0 X , • AT  what  ANGLE  DO  YOU  WISH  THE  FILM  TO  BE  DEPOSITED?'/ 

c.iox, 'please  enter  angle  in  degrees  -  enter  in  real  format.’//) 

REAC  #, OLCANG 

c 

WRITE  (6,0137) 

0137  FORMAT  (lOX.'OO  YOU  WISH  TO  VARY  THE  ANGLE  OF  INC ICENCE? ' // . 1 IX , • I 
CF  YES  -  ENTER  1 '/, 1 1 X ,' OTHER Wl S £  -  HIT  RETURN'//) 

C 

REAO  (5,0115)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  0133 

ELSE 

GOTO  0141 

ENOIF 

C 

0133  WRITE  (6.0135) 

0139  FORMAT  (lOX.'HCW  MANY  DEGRESS  (FULL  ANGLE)  00  YOU  WISH  TO  VARY  TM£ 
C  ANGLE?'/, 10X, 'ENTER  DEGREE  S  IN  REAL  FORMAT.'//) 


o  o  o 


READ  *,  ANG VAR 
WRITE  <6,0140) 

0140  FORMAT  < 2< IX/ ) , 1  OX , • WHAT  SIZE  ANGLE  00  YOU  WISH  TO  INCREMENT  THE  A 
CNGLE  OF  INCIOSNCE?*/,10X, 'ENTER  IN  REAL  FORMAT.'//) 

REAO  $ , ANGINC 

•QUESTIONS  PERTAINING  TO  THE  MOBILITY  OF  THE  DISKS 
0141  WRITE  <6,0142) 

0142  FORMAT  < 2< IX/) , 10X , • DO  YOU  WISH  ANYTHING  OTHER  THAN  NORMAL  MOBILIT 
C Y ? ' //  ,  11X » • IP  YES  -  ENTER  1 • / , 1 IX , ' OTHE R WISE  -  HIT  RETURN'//) 

c 

REAO  <5,0115)  IVAR 
IF  < I V AR  ,EQ.  1)  THEN 
GOTO  1143 

ELSE 

GOTO  0155 

ENCIF 

C 

1143  WRITE  <6,1144) 

1144  FORMAT  <10X,'IF  YOU  WISH  NO  MOBILITY  -  ENTER  l'/,iOX,'IF  YOU  WISH 
CEXTRA  MOBILITY  -  ENTER  2'//) 

C 

1145  READ  <5,0115)  IVAR 

IF  < I V A R  .  EQ.  1)  THEN 
NOMOB  *  1 
GOTO  0155 

ELSE  IF  <IV AR  .EQ.  2)  THEN 
GOTO  0143 
ELSE 

PRINT*, 'TRY  AGAIN* 

GOTO  1145 
ENOIF 
C 

0143  WRITE  <6,0144) 

0144  FORMAT  <10X, 'ENTER  THE  NUMBER  OF  QISKS  TO  BE  DEPOSITED  BEFORE  AN  E 
CXTRA'/.lOX, 'M03ILITY  ITERATION  TAXES  PLACE  -  ENTER  IN  INTEGER  FORM 
CAT.'//) 

READ  * , 3 NC V A R 
C 

C  QUESTION  PERTAINING  TO  WHETHER  OR  NOT  AN  IMPURITY  SHOULD  BE  DEPOSITED. 
C 

0155  WRITE  <6,0156) 

0156  FORMAT  <2<1X/),10X,'00  YOU  WISH  AN  IMPURITY  TO  BE  CEPOSITEO? • // ,  1 1 
C  X , • I F  YES  -  ENTER  1 ' / , 1 1 X , • OTHE R HI SE  -  HIT  RETURN'//) 

C 

READ  <5,0115)  IVAR 
IF  < I V A R  .EQ.  1)  THEN 
IMPUR  «  1 
ELSE 

GOTO  0146 
ENCIF 
C 

C  QUESTIONS  PERTAINING  TO  The  size  of  disk  deposited 
C 

0146  WRITE  <6,0147) 

0147  FORM A 7  < 2< IX/) , 10X, '00  YOU  WISH  TO  USE  TWO  DIFFERENT  SIZE  CIS*?'// 
C.ilX.'IP  YES  -  ENTER  1 • / , 1 1 X  ,  • 0  THE R W  I S 6  -  HIT  RETURN'//) 
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IP  (IVAR  .EQ.  1)  THEN 
LAYER  *  1 
GOTO  0148 

ELSE 

GOTO  1150 

ENOIF 


0148  WRITE  (6,1149) 

1149  FORMAT  CIQX.'OO  YOU  WISH  THE  LARGE  DISK  TO  BE  OEPQSITSO  FIRST?'//, 
Cl IX ,' IF  YES  -  ENTER  1'/,11X, 'OTHERWISE  -  HIT  RETURN'//) 

REAO  (5,0115)  IV AR 
IF  (IVAR  .EQ.  1)  THEN 
SIZE  *  2 
ENOIF 
C 

WRITE  (6,0149) 

0149  FORMAT  (lOX.'HOW  MANY  ROWS  HIGH  (IN  TERMS  OF  FILM  MATRIX)  00  YOU  W 
CISH'/,10X, ' THE  LAYERS  TC  BE?  -  ENTER  IN  INTEGER  FORMAT.'//) 

REAO  #,HGTLAY 
C 

C  PRINTING  OUT  QUESTION  VARIABLES 
C 

1150  WRITE  (6,0150) 

0150  F0RM4T  (12(2X/)) 

PRINT  4, 'THc  NUMBER  OF  DISK  TO  8E  OEPOSITEO  »',NGD5K 

PRINT  4, 'THE  ANGLE  OF  DEPOSITION  WILL  BE  * ' , CLDANG , • D EGRE ES • 

C 

IF  (ANGVAR  .EQ.  0.0)  THEN 

PRINT  4, 'NO  ANGLE  VARIATION  WILL  TAKE  PLACE.' 

ELSE 

PRINT  4,  'THE  AMOUNT  OF  ANGLE  VARIATION  ANGVAR ,• DEGREES . ' 
PRINT  4 , • THE  ANGLE  INCREMENT  «», ANGINC ,' DEGREES . • 

ENOIF 


IF  (NOMOB  .EQ.  1)  THEN 

PRINT*, 'DEPOSITION  WILL  OCCUR  WITH  NO  RELAXATION.' 

ELSE 

IF  (3NCVAR  .EQ.  0)  THEN 

PRINT  *, 'DEPOSITION  WILL  OCCUR  WITH  NORMAL  RELAXATION.' 

ELSE 

PRINT  * , ' AN  EXTRA  MOBILITY  WILL  BE  GIVEN  AFTER  EVERY  •  .3NCVAR 
C,'  ITERATIONS. ' 

ENOIF 

ENOIF 

IF  ( IMPUR  .EQ.  0)  THEN 

PRINT  *, 'NO  IMPURITY  WILL  BE  OEPOSITEO.' 

ELSE 

PRINT  4, 'ONE  IMPURITY  WILL  BE  OEPOSITEC  RANOOMLY.' 

ENOIF 


IF  (LAYER  .EQ.  0)  THEN 

PRINT  4, 'ONLY  CNS  SIZE  DISK  WILL  BE  OEPOSITEO.' 

Else 

IF  (SIZE  .EQ.  1)  THEN 

PRINT  4, 'TWO  OIPFSRSNT  SIZE  DISK  WILL  BE  OEPOSITEO  -  THE  PI 
CAST  BEING  SMALL' 

ELSE 

PRINT  4,'TWC  DIFFERENT  SIZE  DISK  WILL  BE  DEPOSITED  -  THE  FI 
CAST  3EING  LARGS' 
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ENOIF 

PRINT  *  » 1  THE  06PTH  OF  EACH  LAYER  WILL  BE  ',HGTLAY 
SNOIF 
C 

C  LAST  CHECK  BEFORE  DEPOSITION 
C 

WRITS  (6,0151) 

0151  FORMAT  (  2( IX/ ) , 10 X , • 00  YOU  WISH  TO  CHANGE  ANY  OF  THE  ABOVE  VARIABL 
CES? ' // »1 IX  » ' IF  YES  -  ENTER  1'/*11X, 'OTHERWISE  -  HIT  RETURN'//) 

C 

REAO  (5,0115)  IVAR 
IF  (IVAR  .10.  1)  THEN 
GOTO  0132 

ENOIF 

C 

WRITE  (6,0153) 

0153  FORMAT  (10X, '3EP0SITI0N  VARIABLES  HAVE  NOW  BEEN  SET. • / , 10X , • OEPOSI 
CTION  IS  REAOY  TO  OCCUR! I'//) 

PRINT  *,'  HIT  RETURN' 

REAO  (5,0115) 

GOTO  0110 
C 
C 

C  **********  DEPOSITION  OF  FILM  ********** 

C 

C  THIS  PART  OF  THE  SUBROUTINE  IS  THfi  HEART  OF  THE  ENTIRE  PROGRAM.  ITS 
C  PURPOSE  IS  TO  CEPOSIT  THE  "FILM"  ONTO  THE  FABRICATED  "SUBSTRATE". 

C  ALL  NECCESSARY  VARIABLES  TO  RUN  THIS  SECTION,  SHOULO  HAVE  BEEN  READ  IN 
C  BY  EARLIER  SUBROUTINES  OR  SECTIONS.  FOR  MORE  INFORMATION  ON  HOW  THIS 
C  PART  OF  THE  SUBROUTINE  WORKS  CONSULT  THE'  THESIS  -  AN  IMPROVED  MODEL  OF 
C  THIN  FILM  GROWTH  (AFIT  -  SCHOOL  OF  ENGINEERING) . 

C 

0160  WRITE  (6,0161) 

0161  FORMAT  (lOX.'HAVE  ALL  THE  DEPOSITION  VARIABLES  BEEN  REAO  IN  ANO'/, 
Cl  OX ,' THE  FILM  MATRICES  IN  IT  I ALI 2E0?  •  // , 1 IX , • IF  YES  -  ENTER  1'/,11X 
C, 'OTHERWISE  -  MIT  RETURN'//) 

C 

REAO  (5,0115)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  0162 

ELSE 

GOTO  0110 

ENOIF 

C 

C  BEGINNING  OF  DEPOSITION 
C 

0162  WRITE  (6,0163) 

0163  FORMAT  (22X, 'DEPOSITION  PROCESS  HAS  NOW  BEGUN!!*///) 

C 

C  SETTING  OF  THE  SEEO  FOR  RANDOM  NUMBER  GENERATION. 

C 

OSEEO  «  123457. CO ' 

OSEEC  *  (OSccO*(OAT£/iOOO.O)> 

C 

C  INITIALIZATION  OF  DEPOSITION  VARIABLES 
C 

3NCCTR  «  0 
3CTR  «  0 
IMPSET  »  0 
HGTIMP  «  o 
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SET  *  0 
HOSIZE  »  SIZE 
YCK  «  0 
IMPHGT  a  0 
IMPX  *  0 

YHGT  *  AINT ( HEIGHT)  ♦  3 
Y FILM  a  AINT (HEIGHT) 

ANGMAX  >  OLDANG  *  ANGVAR/2 
ANGMIN  a  OLOANG  -  ANGVAR/2 
ANG  >  OLOANG 
NO  *  0 
P  *  1 
HCTR  a  0- 
C 

C  DETERMINATION  OF  WHEN  IMPURITY  HILL  BE  DEPOSITED 
C 

IF  CIMPUR  .EQ.  1)  THEN 
RANNUM  a  GGUBFSCOSEEO) 

RANDY  a  (RANNUM  *  179)  ♦  YHGT 
HGTIMP  a  AINTCRANOY) 

C  HGTIMP  a  io 

SNOIF 


c  *#*»*COLLISION  POINT  DETERMINATION*#*#* 

C 

C  RANDOM  X  DETERMINATION 
C 

DO  0170  OEPCTR  a  1,N0QSK,1 

IF  C YHGT  .GT.  199)  GOTO  0201 
RANNUM  a  GGUBFSCOSEEO) 

RANOX  a  (RANNUM  *  300)  *  1 
C 

C  ANGLE  DIRECTION  CHANGE  SECTION 
C 

IF  (ANGVAR  .SC.  0.0)  THEN 
GOTO  0175 

ELSE 

ANG  >  ANG  ♦  P# ANGINC 
IF  (ANG  .GT.  ANGMAX)  THEN 

?  a  -1 

ANG  a  ANGMAX  -  ANGINC 
ELSE  IF  (ANG  .LT.  ANGMIN)  THEN 
P  »  1 

ANG  a  ANGMIN  ♦  ANGINC 

END  I  F 

ENOIF 

0175  RANG  a  ANG/57 .2957751 

CCS ANG  a  COSC RANG) 

TANANG  «  T  AN( RANG) 

SC2ANG  *  TANANG**2  ♦  1 
C 

C  SIZE  AND  COUNTER  MODIFICATION  FOR  IMPURITIES 
C 

IF  (IHPUR  .EQ.  1)  THEN 

IF  (SET  .EQ.  1)  THEN 
SIZE  »  HOSIZS 
GOTO  0177 
ENOIF 

IF  (YCK  .EQ.  0)  THEN 
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IF  cyhgt  .GE.  HGTIMP)  then 
tck  »  i 

IMPSET  *  1 
SET  *  1 
HOSIZE  a  SIZE 
SIZE  »  3 
CTR  »  15 
GOTO  0180 
ENOIP 
ENOIF 
ENOIF 
C 

C  DETERMINATION  OF  SCAN  COUNTER 
C 

0177  IF  (IMPSET  .EC.  0)  THEN 

IF  (LAYER  «EQ.  0)  THEN 
CTR  »  3 
GOTO  0130 

ELSE 

CTR  a  6 
GOTO  0179 
ENCIP 
ELSE 

IF  (LAYER  .£Q.  0)  THEN 
CTR  a  14 
GOTO  0180 
ELSE 

CTR  a  15 
GOTO  0179 
ENOIF 
END  IF 
C 

C  SIZE  CHANGING  SECTION  FOR  DIFFERENT  SIZE  DISK 
C 

0179  IF  (((YHGT-1)  -  YFILM)  .GE.  HGTLAY)  THEN 

YFILM  a  YMGT 
IF  (SIZE  .EQ.  1)  THEN 
SIZE  ■  2 

ELSE  IF  (SIZE  .SO.  2)  THEN 
SIZE  «  1 
ENCIF 

HOSIZE  a  SIZE 
ENOIF 
C 

0130  Y2  a  0 

X2  »  0 

0*INT  »  0 
0 Y I  NT  a  o 
XCCLL  »  0 
YCOLL  ■  0 
A2  a  0,0 

00  0131  STRCTR  a  QfCTR,l 
C 

C  DETERMINATION  CF  X  INTERCEPT 
C 

IP  (IMPSET  .EC.  1)  Then 
IP  (SIZE  .EG.  1)  THEN 

XINCT  a  RANCX  ♦  (  6 . 3 6 3 96  1  /CC S ANG ) 
ELSE 

XINCT  a  RANOX  ♦  (  7 . 0 7  1  067 3/CC S ANG ) 
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\  ENCIF 

V  ELSE 

XINCT  *  RANOX  ♦  (SIZE  $  1.4142136)/C0SANG 
ENOIF 
C 

c  determination  of  streamer  intercept 
c 

SINCT  s  RANOX  ♦  .942309416  *  (STRCTR/CQSANG) 

c 

XPOS  *  YHGT  *  TANANG  ♦  SINCT  -  1 
OXPOS  *  XPOS 
COLL  *  0 

c 

C  BEGINNING  OF  STREAMER  SEARCH 
C 

DO  0190  Y  *  YHGT, 1,-1 

XPOS  a  Y  3  TANANG  +  CSINCT  -  1) 

ISHFT  »  CCAINTCXPOS/30D)  *  300) 

NEWX  *  AINTCXPOS  -  ISHFT) 

CLDX  a  AINT(OXPOS  -  ISHFT) 

CC  0195  X  a  CLOX , 1 *-l 

IF  (FILM(X.Y)  .NE.  0)  THEN 

IF  ((X  .EC.  OXINT) . ANO.C Y  .EQ.  0 Y INT )  )  GOTO  0194 
NXINCT  a  XINCT  -  ISHFT 

G  *  (NXINCT-XCOQRO(X,Y))+(TANANG  *  YCCORO ( X  ,  Y ) ) 

C 

C  COLLISION  CIRCLE  RADIUS  DETERMINATION 
C 

IF  (FILM(X.Y)  .EG.  1)  THEN 
IF  (SIZE  .EQ.  1)  THEN 
RA02  a  2 

ELSE  IF  (SIZE  .EQ.  2)  THEN 
R AC 2  a  4.5 

ELSE  IF  (SIZE  .EQ.  8)  THEN 
RA02  a  40.5 
ENOIF 

ELSE  IF  ( F I LM ( X , Y )  .EG.  2)  THEN 
IF  (SIZE  .EQ.  1)  THEN 
RA02  a  4.5 

ELSE  IF  (SIZE  .EQ.  2)  THEN 
RAC2  a  s 

ELSE  IF  (SIZE  .EG.  3)  THEN 
R  AO  2  a  50 
END  IF 

ELSE  IF  (FILM(X,Y)  .EG.  3)  THEN 
IF  (SIZE  .EQ.  1)  THEN 
R AD 2  a  40.5 

ELSE  IF  (SIZE  .EQ.  2)  THEN 
R  AC  2  *  50 
ENOIF 
.  ENOIF 
C 

C  TESTING  eOR  PCSSI3LE  COLLISION  PARTNER 
C 

T  E  S  T 1  •  ( R  AD  2  a  SC2ANG  -  G**2) 

IF  (TEST1  .LT.  0.0)  THEN 
GCTO  0194 
ELSE 

F  a  YCCCRCCX.Y)  -  TANANG*(NXlNCT-XCOORO(X,Y)) 
YI  a  (F  ♦  SGRT(TEST1))/SC2ANG 
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QXINT  *  X 
QYINT  a  Y 

IF  C Y I  .ST.  Y2  )  THEN 
Y 2  a  YI 
A2  a  NXINCT 
XCCLL  *  X 
YCCLL  »  Y 
5N0IF 
COLL  =  1 
ENDIF 
ENOIF 


0194 

IF  CX  .EQ.  NEWX) 

GOTO 

0195 

CONTINUE 

0196 

CXPOS  »  XPOS 

IF  (COLL  .EQ.  1)  GOTO 

0131 

0190 

CONTINUE 

0131 

CONTINUE 

C 

C  DETERMINATION  CF  X  COLLSION  POSITION 
C 

X2  a  Y2  #  TANANG  ♦  A2 
C 

C  TEST  FOR  NO  MOBILITY 
C 

IF  CNGMG3  .EQ.  1)  THEN 
RESTX  a  X2 
RESTY  a  Y2 
GOTO  0221 
ENOIF 
C 

C  *****REST  POINT  DETERMINATION***** 

C 

C 

0205  SHORT  a  10000 

MIN  a  o 

MAX  a  o 

c 

C  INTIALIZATION  FOR  MOBILITY  ROUTINE 
C 

IF  ( 3NC VAR  .NE.  0)  THEN 
00  0207  B  a  1,6,1 
RSNCCB)  a  10000 
XBNCC3)  a  0.0 
YSNCC3)  »  0.0 
0207  CONTINUE 

ENCIF 
C 

C  DETERMINATION  CF  THE  SIZE  OF  THE  ARRAY  SEARCH 
C 

IF  (IMPSET  .EQ.  0)  THEN 
IF  (LAYER  .EQ.  0)  THEN 
MIN  a  -3 

MAX  a  3 

GCT3  0210 
ELSE 

WIDTH  a  FILMCXCOLL, YCCLL)  ♦  SIZE 
IF  (WIOTH  .EC.  2)  THEN 
min  *  -A 

MAX  a  4 

ELSE  IF  (WIDTH  .EC.  3)  THEN 


o  o  o 


MIN  a  -5 
MAX  *  5 

ELSE  IF  (WIDTH  .EQ.  4)  THEN 
MIN  «  -6 
MAX  a  6 
END  IF 
ENOIF 
ELSE 

IF  (LAYER  .EQ.  0)  THEN 
MIN  a  -3 
MAX  a  3 
GOTO  0210 
ELSE 

WIDTH  a  FILM(XCQLL.TCQLL)  ♦  SIZE 
IF  (WIDTH  .EQ.  2)  THEN 
MIN  a  -3 
MAX  a  3 

ELSE  IF  (WIDTH  .EC.  3)  THEN 
MIN  a  -10 
MAX  a  10 

ELSE  IF  (WIOTH  .EC.  4)  THEN 
MIN  a  -10 
MAX  a  10 

ELSE  IF  (WIDTH  .EC.  9)  THEN 
MIN  a  -9 
MAX  a  9 

ELSE  IF  (WIOTH  .EQ.  10)  THEN 
MIN  *  -10 
MAX  a  10 
ENOIF 
ENOIF 
ENOIF 

BEGINNING  OF  SEARCH  ARCUNO  COLLISION  PARTNER 

0210  CO  0220  J  »  MIN, MAX, l 

Y  a  0 

Y  a  j  ♦  YCOLL 

IF  ((Y  .LT.  1)  .OR.  (Y  .GT.  YHGT) )  GOTO  0220 
00  0230  I  a  MIN, MAX, 1 

X  a  0 

X  a  I  ♦  XCOLL 
3HFT  a  0 

IF  (X  .LT.  1)  THEN 
X  a  x  ♦  300 
SHFT  a  1 

ELSE  IF  (X  .GT.  300)  Then 
X  a  x  -  300 
SHFT  a  300 
ENOIF 
C 

C  INITIALIZATION  OF  REST  POINT  VAR IA9LES 
C 

IF  <film(x,y)  .ne.  0)  then 

A  a  0.0 
3  »  0.0 
C  a  0.0 
0  a  0.0 
R 1  a  0.0 

R 2  a  0.0 
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H  »  0.0 
K  »  0.0 
TEST2  «  0.0 
RINC  >  0.0 
C 

C  SETTING  OP  REST  CIRCLE  RAOII 
C 

IP  <  SIZE  .60.  1)  THEN 
RINC  ■  .707136781 
ELSE  IP  (SIZE  . EQ .  2)  THEN 
RINC  >  1.414213562 
ELSE  IP  (SIZE  -EQ.  8)  THEN 
RINC  *  5.6568542 
eNOIP 
C 

c  DETERMINATION  op  rest  point  variables 
c 

R1  «  (PILMCXCOLLiYCOLL)  *  .707106781)  ♦  RINC 
R2  *  <FILM(X,V)  *  .707106781)  ♦  RINC 
C 

IP  (SHPT  .S3.  0)  THEN 

H  «  XC00R0C  X , Y)  -  XC00RCCXCCLL, YCOLL) 

ELSE  IP  (SHPT  .EQ.  1)  THEN 

H  «  XCQORQCX.Y)  -  300.0  -  XC0QR0C XCCLL, YCOLL) 
ELSE  IP  < SHPT  .EQ.  300)  THEN 

H  *  XCOOROCX.Y)  ♦  300.0  -  XCOORDC XCCLL . YCOLL ) 
ENOIP 
C 

K  *  YCCQROC  X  »  Y )  -  YC OORO C XCOLL , YCCLL ) 

A  *  (R2*R2  -  R1PR1  -  H9H  -  K*K)/2 
8  «  A9A  -  CH*M)*CR1*R1) 

C  »  <***  ♦  H*H) 

IF  CC  .£3.  0.0)  GOTO  0230 
T6ST2  «  ( A3 A6K4K  -  <C*8)> 

C 

C  TESTING  FOR  POSSIBLE  REST  PARTNERS 
C 

IF  CTEST2  .LT.  0.0)  THEN 
GOTO  0230 
C 

C  DETERMINATION  OF  POSSIBLE  REST  POINTS 
C 

ELSE 

00  0250  L  ■  -1,1,2 

Y 1  »  (-CA3K)  ♦  L»S3RTCTSST2))/C 
CO  0260  M  »  -1,1,2 

o  ■  A 4 A  -  (X3X)#(R1*RJ  ) 

TEST3  *  { H9H3A4A  -  ( C30  )  ) 

IF  CTEST3  .LT.  0.0)  GOTO  0250 
XI  «  C-CM*A)  ♦  M*SCRT(TEST3))/C 
TEST4  *  <(X1  -  H)#$2  ♦  (Tl  -  K>#*2> 

TCL  ■  ASSCTESTA  -  R2«RZ) 

IF  (TCL  .LT.  .005)  THEN 

XOIST  «  XI  ♦  XCCCRCC XCOLL, YCOLL)  -  X2 
YOIST  *  Y1  ♦  YCOORCC XCOLL,  YCOLL)  -  Y2 
R  *  (XCIST4«2  ♦  Y0IST**2) 

C 

C  STORING  OP  VALUES  fqR  M08ILITY  ROUTINE 
C 

IF  (3NCVAR  .NS.  0)  Then 
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OQ  02S5  SNC  *  1,6,1 

IF  (R  .LT.  R  SNC ( SNC ) )  THEN 

IF  C3NC  .63.  6)  SOTO  02ST 
BNC1  «  SNC  ♦  1 
00  02S6  81  >  5  ,  SNC 1,-1 
R3NCCB1)  »  R3NCCE1-1) 

XBNCC31)  *  X3NCCB1-1) 

Y3NCCB1)  «  YBNCCB1-1) 

0256  CONTINUE 

0257  R3NCC3NC)  •  R 

XBNCC8NC)  >  X1+XC0CR0CXC0LL,YCCLL) 
YBNCC3NC)  •  Y1*YC0CR0(XC0LL, YCOLL) 
GOTO  0258 

ELSE 

GOTO  0255 
ENCIF 

0255  CONTINUE 

ENOIF 

C 

C  TEMPORARY  STORAGE  OF  REST  POINT  VALUES 
C 

0258  IF  (R  .LT.  SHORT)  THEN 

TShORT  »  SHORT 
SHORT  «  R 
TRESTX  *  0.0 
TRESTY  •  0.0 
TRESTX  «  RESTX 
TRESTY  «  RSSTY 
RESTX  *  0.0 
RESTY  »  0.0 
TESTY  *  0.0 

TESTY  «  Y1  ♦  YCOOROCXCOLL, YCOLL) 

IF  (TESTY  .LT.  2.0)  GOTO  0260 
RSSTY  »  Yi  ♦  YCOOROCXCOLL, YCOLL) 

RESTX  *  XI  ♦  XCCQRD(XCOLLfYCCLL) 

C 

C  MOVEMENT  OF  OISK  IF  NOT  WITHIN  80UNCARIES 
C 

IF  (RESTX  .GE.  301.0)  THEN 
RESTX  »  RESTX  -  300.0 
ELSE  IF  (RESTX  .LT.  1.0)  THEN 
RESTX  ■  RESTX  ♦  300.0 
ENOIF 
C 

C  TESTING  OF  POSITION  FOR  OCCUPATION 
C 

TPX3  »  OINT (RESTX) 

TPY3  «  OINT(RESTY) 

IF  (FILMCTPX3 ,TPY3)  .N6.  0)  THEN 
IF  (TRESTX  ,N6.  0.0)  THEN 
SHORT  »  TSHORT 
RESTX  »  TRESTX 
RSSTY  «  TRESTY 
ENOIF 
ENOIF 

ORE  ST  X  »  RESTX 
CRSSTY  »  RESTY 


C 


ENOIF 

ELSE 
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************************************ ***************************** 
EXTRA  MOBILITY  ROUTINE 


r.V 


PRINT*,'  • 

PRINT*, *RX  RY' ,R£STX,RESTY 
IF  C8NCVAR  .5Q.  0)  THEN 
GOTO  0221 
ELSE 


TEST  TO  MAKS  SURE  THE  FILM  IS  HIGH  ENOUGH 


Y 3  «  OINTCRESTY) 

IF  (SHE  .EQ.  1)  THEN 

IF  CY3  .LE.  (1+HEIGHT) )  GOTO  0221 
ELSE  IF  (SIZE  .EQ.  2)  THEN 

IF  <Y3  .LE.  (2+HEIGHT))  GOTO  0221 
ELSE  IF  (SIZE  .EQ.  8)  THEN 

IF  ( Y 3  .LE.  (*-*HE!GHT)>  GOTO  0221 
6N0IF 
C 

BNCCTR  «  BNCCTR  ♦  1 
C  PRINT*, BNCCTR 

IF  (BNCCTR  .GE.  SNCVAR)  THEN 
OC  0235  B  >  1,6,1 
C 

IF  (XBNC(B)  .GT.  300.0)  THEN 
X3NC(B)  »  XBNC(B)  -  300.0 
ELSE  IF  (X3NC( B)  .LT.  1.0)  THEN 
X3NC ( 3 )  *  X3NC(3)  ♦  300.0 
ENOIF 
C 

C  INTI AL  CHECK  FCR  REOUNOANCY 
C 

IF  (CX3NC(8>.NE.0RS$TX).AND.CY3NC(3).NE.CRESTY))  THEN 
IF  ((XSNC(3).NE.O.0).ANO.(YBNC(B).N:.0.O))  THEN 
TXBNC  »  A I NT( X  3NC ( 3  )  ) 

TY3NC  «  AINT(Y3NC(3)) 

IF  (FILM(TXSNC,TY3NC)  .EQ.  0)  THEN 
C 

C  BEGINNING  OF  OISTANCE  AREA  CHECK 
C 


CO  0236  JJ  *  MIN, MAX, 1 
T3CK  «  JJ  »  TT3NC 
00  0237  II  *  MIN, MAX, 1 
XSCK  ■  II  ♦  TXSNC 
83MFT  »  0 


C 


IF  (X3CK  .GT.  300)  Then 
X3CK  »  X3CK  -  300 
3  ShF T  •  300 
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ELSE  IF  (X3CK  .LT.  I)  THEN 
XBCK  •  XBCK  ♦  300 
8SMPT  ■  I 
ENOIF 

if  cfilmcxbck.ybcio  .ne.  o)  then 

TBOIST  ■  YCOCROCXaCX.YSCK)  -  Y3NCC3) 

IF  (BSHFT  .EG.  0)  THEN 

XBDIST  a  XC0QRDCX3CK, Y3CK)  -  XBNCC3) 
ELSE  IP  (3SHFT  .EQ.  300)  THEN 

X30IST  >  XC00R0CX3CX, Y3CK ) ♦ 30 0-X3NC CS ) 
ELSE  IP  (3SHFT  .EQ.  1)  THEN 

X30IST  »  XCG0R0CXaCK,Y8CX)-300-X8NCCB) 
ENOIF 

TESTS  ■  (XBDIST**2  ♦  Y3DIST**2) 


c  determining  the  tolerance  OISTANCS 
c 

IF  <FILMCX3CK(Y3CK)  .60.  I)  THEN 
I?  (SIZE  .EC.  1)  THEN 
BTCL  >  1.98 

ELSE  IF  (SIZE  .EQ.  2)  THEN 
3T0L  >  *.45 

ELSE  IF  (SIZE  .EQ.  8)  THEN 
BTCL  »  *0.4 
ENOIF 

ELSE  IF  (FILN(X3CA,YBCK)  .EQ.  2)  THEN 
IF  (SIZE  .EQ.  1)  THEN 
BTCL  *  4.45 

ELSE  IF  (SIZE  .EQ.  2)  THEN 
BTCL  »  7.9 

ELSE  IF  (SIZE  .EQ.  8)  THEN 
8TCL  »  49.0 
ENOIF 

ELSE  IF  (FILM(X8CK«Y3CK)  .EQ.  3)  THEN 
IF  (SIZE  .EC.  1)  THEN 
BTCL  •  40.4 

ELSE  IF  (SIZE  .EQ.  2)  THEN 
STCL  ■  49.0 
ENQIF 
ENOIF 
C 

C  TESTING  FOR  OISTANCE  TOLERANCE 
C 

IF  (TESTS  .LT.  3TCL)  GOTO  0235 
ENOIF 

0237  CONTINUE 

0236  CONTINUE 

RESTX  «  X8NCC3) 

RESTY  »  Y3NCC3) 

.YTEST  *  ORESTY  ♦  .3 
IF  (RESTY  .GT.  YTEST)  GOTO  0235 
ENCCTR  ■  0 
8CTR  *  3CTS  ♦  1 
GOTO  0219 
ENOIF 
ENOIF 
ENOIF 

0235  CONTINUE 

GOTO  0219 
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ENOIF 

ENOIF 

C 

C  LAST  CHECK  ANO  OEPAULT  STRUCTURE 
C 

0218  IF  CRESTY  .LT.  1.0)  THEN 

8CTR  »  9CTR  -  1 
0219  RESTX  »  ORESTX 

RESTY  *  ORESTY 
8NCCTR  »  CSNCVAR  -  1) 

ENOIF 

C  PRINT#, 1 RX  RY* , RESTX, RESTY 

C 

C  t##*#*##*##*####*#**#*###*#*#**#**#*#*####*****##*#*##**#****** 
C 

0221  X3  *  DINT(RESTX) 

Y 3  «  OINTCRESTY) 

C 

C  TESTING  FOR  OVERWRITES 
C 

IF  (FILMCX3.Y3)  .NE.  0)  THEN 
HCTR  »  HCTR  ♦  1 
ENOIF 
C 

C  FIXING  OF  OISK  POSITION 
C 

XCQ0R0(X3,Y3)  «  RESTX 
YC00RCCX3, Y3)  «  RESTY 
FILM(X3,Y3)  «  SIZE 
C 

C  SETTING  OF  THE  STARTING  HEIGHT  FOR  NEXT  DISK 
C 

IF  (YCK  .EQ.  1)  THEN 
YCK  *  0 

TYHGT  »  YHGT  -  CY3  *  6) 

IF  CTYHGT  .GT.  0)  THEN 
YHGT  *  YHGT 
ELSE 

YHGT  ■  Y3  ♦  9 
ENOIF 
ELSE 

IF  (Y3  .GE.  YHGT)  THEN 
YHGT  *  Y3  ♦  3 
ENOIF 
ENOIF 
C 

C  INCREMENTING  THE  OISK  COUNTER 
C 

NO  «  NO  ♦  1 
IF  (NO  .E3.  500)  THEN 
WRITE  (6,0200)  OEPCTR 

0200  FORMAT  (28X,I5,*  QISKS  OEPOSITEO’) 

NC  ■  0 
ENOIF 
C 

C  RESETTING  OF  PARAMETERS  IN  THE  CASE  OF  IMPURITIES 
C 

IF  ( IMP  SET  .50.  1)  THEN 

IF  (YHGT  .GT.  (HGTIMP*50) )  THEN 
IHPSET  «  0 


•  %  "V  ■-  *_■  V." 


ENOIP 

ENOIP 


C  FINOING  THE  POSITION  OP  IMPURITY 
C 

IP  (SIZE  .EQ.  3)  THEN 
IMPHGT  *  Y 3 
IMPX  *  X3 
ENOIP 
C 

C  MESSAGE  CENTER  FOR  ENO  OP  DEPOSITION 
C 

0170  CONTINUE 

0201  WRITE  (6*0202)  (OEPCTR  -  1) 

0202  FORMAT  ( 2C1X/ ) , 1 SX , • OEPOSIT ION  COMPLETEO  WITH  •  ,  15  »  •  OISKS  06P0SIT 
CEO! •///) 

C 

WRITE  (6*0225)  HCTR 

0225  FORMAT  (23X, 'THERE  WERE  *,I2,*  OVERLAPPING  DISKS'//) 

C 

IP  (3NCVAR  .NE.  0)  THEN 
WRITE  (6*0229)  BCTR 

0229  FORMAT  (13X,I5*'  OISKS  WERE  GIVEN  AN  EXTRA  MOBILITY'//) 

ENOIP 

C 

IF  (IMPHGT  .EQ.  0)  THEN 
GOTO  0226 
ELSE 

WRITE  (6,0227)  IMPHGT 

0227  FORMAT  (22X,'THE  Y  HEIGHT  OP  THE  IMPURITY  IS  »,I3) 

WRITE  (6,0223)  IMPX 

0228  FORMAT  (21X,'THE  X  POSITION  OF  THE  IMPURITY  IS  ',13//) 

ENOIP 

C 

0226  WRITE  (6,0222) 

0222  P0RMATC9X, 'DEPOSITION  VARIABLES  ARE  RE AO Y  TO  BE  ANALYZED  ANO/OR  ST 
CORED. '//) 

PRINT  3,'  HIT  RETURN* 

REAO  (5,0115) 

RETURN 


333*3333333333333333333333333333 

*  3 

3  ANALYSIS  SUBROUTINE  3 

*  3 
3*333333**3333333333333333333333 

this  subroutine  consist  OF  2  AREAS:  That  is  The  calculating  of  the 
DENSITY  of  The  film  matrices  ANC  The  DETERMINATION  Or  THE  ANGLE  OF 
The  OSPOSITEO  FILM.  EACH  SECTION  CONTAINS  A  BRIEF  DESCRIPTION  OF 
ITS  PURPOSE  AnC/OR  ITS  APPROACH. 


SUBROUTINE  ANAL(FILM,XCOORC,YCOORO) 
LOCAL  VARIABLES 


***  «  ’  l**#  k  ’A  , 
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REAL  XCOORC(300*200),YCOORD(30Q,200),RSMLON»RLRGON,T0TOEN,CORR(90) 
REAL  LOCGEN,LLRGON*LSMLQN»SLNOSN,LLNDEN,CCRFAC»ISHFT ,XP0S,QXP0S 
REAL  ANG, R ANG, T AN ANG.TLNOEN, LARGE  1 ,L ARG E2 , LARGE3 , AVGANG 
INTEGER  FILM(300,200),lVAR,$MLCTR,LRGCTR,TQTCTR,SONCTR,L0NCTR 
INTEGER  L0CCTR,30THGT,T0PHGT,SUB,CELCTR,LDKCTR,SDKCTR,MAXCNT,A 
INTEGER  TEST.ANGLElt ANGLE 2 , ANGLE3 , OENCK , LEFT .RIGHT 
C 
C 

c  **********  ANALYSIS  INTERFACE  LOOP  ********** 

C 

0400  WRITE  (6,0401) 

0401  FORMAT  ( 19( 2X/ ) , 27X , • ANALY S IS  INTERFACE  LOCF •///, 20X, • ENTER  AFFROP 
CRIATEO  NUM6ER  TO  CONTINUE '/, 2(2X/)  ,25X  ,»  1  -  CALCULATION  OF  OENSITY 
C ' /, 2  5X  , '2  -  CALCULATION  OF  ANGLE • / ,2SX  ,  • 9  -  EXIT  TO  MAIN  INTERFACE 
C  LG OP* /*25X, 'CTRL/C  -  EXIT  TO  OPERATING  SY STEM ' / , 6( 2X/ ) ) 

C 

0405  REAO  (5,0406)  IVAR 
0436  FORMAT  (ID 
C 

IF  (IVAR  .EQ.  1)  THEN 
GOTO  0410 

ELSE  IF  (IVAR  .EQ.  2)  THEN 
GOTO  0430 

ELSE  IF  (IVAR  .EQ.  9)  THEN 
RETURN 
ENOIF 

WRITE  (6,0407) 

0407  F0RMAK25X,  ‘SORRY  CHARLIE!!  -  TRY  AGAIN'///) 

GOTO  0405 


$ 
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**********  CALCULATION  OF  OENSITY  ********** 

THIS  PART  OF  THE  SUBROUTINE  IS  LSEO  TO  CALCULATE  THE  OENSITY  OF  THE  OE- 
P0SIT6C  FILM.  THIS  IS  ACCOMPLISHED  3Y  SCANNING  THE  ENTIRE  FILM  MATRIX 
ANC  CHECKING  FOR  OCCUPATION  BY  EITHER  SMALL, LARGE  QR  IMPURITY  OISKS. 
AFTER  COUNTING  THE  NUMBER  OF  OISK  THE  AREA  TAKEN  UP  3Y  THE  DISK  IS  COM¬ 
PARED  TO  THE  TOTAL  POSSIBLE  AREA  THUS  CREATING  THE  PACKING  DENSITY.  A 
LOCAL  OENSITY  IS  ALSO  CALCULATED  TO  STUOY  THE  AFFECTS  CF  CHANGING  CER- 
T I  AN  PARAMETERS. 

0410  WRITE  (6,0411) 

0411  FORMAT  ( 10X , • IT  IS  ASSUMEO  UPON  ENTERING  THIS  SECTION  OF  THE  SU3R0 
CUTINE  THAT'/, 10X, ' YOU  WISH  TO  FINO  THE  DENSITY  OF  THE  CURRENT  FILM 
C  MATRICES. '/,  10X  ,' IS  THIS  CORRECT?'//, UX, 'IF  YES  -  ENTER  1'/,11X, 
C’OTHERWISE  -  HIT  RETURN'//) 

C 

REAO  (5,0406)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
ANGCK  >  0 
GOTO  0412 
ELSE 

GOTO  0400 
ENOIF 

c 

0412  WRITE  (6,0413) 

0413  FORMAT  ( 1  OX , ' IF  YOU  WISH  TO  LOOK  AT  VARIATIONS  IN  Y  -  ENTER  l'/, 10 
CX.'IF  YOU  WISH  TO  look  AT  VARIATIONS  IN  X  -  ENTER  2'/,10X,'IF  YOU 
CWISH  TO  LOCK  AT  VARIATIONS  IN  3CTH  K  ANO  Y  -  ENTER  3'//) 
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RE AQ  (5,0490)  OENCK 
0490  FORMAT  (ID 
C 

IF  (OENCK  .EQ.  1)  THEN 


GOTO 

0415 

ELSE  IF 

(OENCK  .EQ.  2) 

THEN 

GOTO 

0430 

ELSE  IF 

(OENCK  .EQ.  3) 

THEN 

GOTO 

0415 

ELSE 

PRINT*, *  TRY  AGAIN' 

GOTO  0412 
ENOIF 
C 

0415  WRITE  (6,0416) 

0416  FORMAT  (lOX.'IN  ORCeR  TO  CHECK  FOR  VERTICAL  VARIATIONS  IN  DENSITY 
CTWO  HEIGHTS'/, 10X, 'NEED  TO  Be  SPeCIFIEC.  PLEASE  ENTER  THE  Y  VALUE 
CS  IN  INTEGER  FORMAT.'//) 

C 

0434  PRINT*,  • ENTER  BOTTOM  HEIGHT* 

0413  RE AO  (5,0417)  80THGT 
0417  FORMAT  (13) 

IF  (30THGT  .LT.  1)  THEN 

PRINT*,'  BOTTOM  HEIGHT  TOO  LOW  -  TRY  AGAIN' 

GOTO  0418 
SNOIF 

PRINT*, • ENTER  TOP  HEIGHT' 

0419  READ  (5,0417)  TQPHGT 

IF  (TOPHGT  .GT.  200)  THEN 

PRINT*,'  TOP  HEIGHT  TOO  HIGH  -  TRY  AGAIN' 

GOTO  0419 
ENOIF 

PRINT#,'  ' 

C 

IF  (CANGCK  .Eg.  1). OR. (OENCK  .EG.  1))  THEN 
LEFT  »  1 
RIGHT  »  300 
PRINT*, 'WORKING' 

GOTO  0433 
ENOIF 
C 

0430  WRITE  (6,0431) 

0481  FORMAT ( 1  OX , • IN  CROER  TO  CHECK  FOR  HORIZONTAL  VARIATIONS  IN  OENSITY 
C  TwO  WlOTHS'/flOX, ' NEEO  TO  SE  SPECIFIEC.  PLEASE  ENTER  THE  X  VALUE 
CS  IN  INTEGER  FORMAT.'//) 

C 

PRINT*, • ENTER  LEFT  30UN0ARY • 

0482  RSAO  (5,0417)  LEFT 

IF  (LEFT  .LT.  1)  THEN 

PRINT*,'  LEFT  BOUNCARY  TOO  FAR  LEFT  -  TRY  AGAIN' 

GOTO  0432 
SNCIF 

PRINT#, • ENTER  RIGHT  30UNCARY* 

0483  READ  (5,0417)  RIGHT 

IF  (RIGHT  .GT.  300)  Then 

PRINT#,*  RIGHT  BOUNDARY  TOO  FAR  RIGHT  -  TRY  AGAIN' 

GOTO  0433 
ENOIF 
C 

Is  (OENCK  .SQ.  2)  Then 


8? 
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TOPHGT  a  200 
3QTHGT  «  1 
ENOIF 

PRINT#, •  • 

PRINT#, • WORKING' 

INTIALIZATION  CP  DENSITY  VARIABLES 

3433  TOTCTR  »  0 
LRGCTR  «  0 
SHLCTR  »  0 
IMPCTR  «  0 
IDNCTR  «  0 
SONCTR  «  0 
LONCTR  «  0 
LOCCTR  »  0 

COUNTING  OF  OISK  PGR  OENSITY  CALCULATIONS 

00  0420  Y  *  1,200,1 
03  0421  X  »  1,300,1 

IF  CFILMCX.Y)  .Efi.  1)  THEN 
SMLCTR  a  SMLCTR  ♦  1 
ELSE  IF  CFILMCX.Y)  .EG.  2 >  THEN 
LRGCTR  *  LRGCTR  ♦  1 
ELSE  IF  CFILM(X.Y)  . EQ .  8)  THEN 
IMPCTR  »  IMPCTR  ♦  1 
ENOIF 

IF  CCY  .GE.  30THGT). ANO.C Y  .LE.  T0PHGT3)  THEN 
IF  CCX  .Ge.  LEFT3.ANC.CX  .LS.  RIGHT))  THEN 
IF  CFILMCX.Y)  .EQ.  1)  THEN 
SDNCTR  «  SONCTR  ♦  1 
ELSE  IF  CFILMCX.Y)  .E3.  2)  THEN 
LONCTR  »  LONCTR  ♦  1 
ELSE  IF  CFILMCX.Y)  .EQ.  8)  THEN 
ICNCTR  a  IONCTR  ♦  1 
ENOIF 

LOCCTR  a  LOCCTR  ♦  1 
ENOIF 
ENOIF 

0421  CONTINUE 
0420  CONTINUE 
C 

C  FINOING  THE  TOTAL  NUMBER  OF  DISK 
C 

TOTCTR  a  SMLCTR  ♦  LRGCTR  ♦  IMPCTR 
IF  (ANGCK  .EQ.  1)  GOTO  0429 
C 

C  OUTPUT  FOR  THE  NUMBER  CP  OISK 
C 


0422 

WRITE  C  6  •  0 4 2 2 )  SMtCTR 
FORMAT  CIOX.'ThE  number 

OF 

SMALL 

OISKS 

OEPCSITSC  a  ',15) 

0423 

WRITS  C6.0423)  LRGCTR 
FORMAT  CIOX.'THE  NUMBER 

Or 

LARGE 

DISKS 

QEPC SITED  »  »,I5) 

1423 

WRITS  (6,1423)  IMPCTR 
FORMAT  CIOX.'THE  NUMBER 

OF 

IMPURITIES 

»  '  ,  11  ) 

0424 

WRITE  (6,0424)  TOTCTR 
FORMAT  CIOX.'THE  TCTAU 

NUMBER  OF 

OISK 

OE POSITED  *  ',15) 

C 
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C  CALCULATIONS  OF  THE  DENSITIES 
C 

0429  RSMLON  «  C SMLCTR41 . 5707963*100 )/60000 
RLR60N  «  (LRGCTR*6.2831853*100)/60000 
RIMPON  *  (IMPCTR*100. 53096*100/60000 
TOTOeN  >  RSMLON  ♦  RLRGON  ♦  RIMPDN 
LSMLON  »  ($ONCTR*1.7320508O3*100)/L0CCTR 
LLRGON  «  (L0NCTR*6.2331353*100)/L0CCTR 
LIMPDN  ■  (I0NCTR*1Q0.53096*100)/LQCCTR 
LOCOEN  *  LSMLON  ♦  LLRGON  ♦  LIMPDN 
C 

IP  (ANGCK  .EQ.  1)  GOTO  0435 
C 

PRINT*,'  • 

C 

C  OUTPUT  OP  THE  CENSITIES 
C 

WRITE  (6,0425)  RSMLON 

0425  PCRMAT  (10X, 'SPACE  OCCUPIED  BY  SMALL  DISK  »  ',F5.2,'*') 

WRITE  (6,0426)  RLRGON 

0426  FORMAT  (10X, 'SPACE  OCCUPIED  8Y  LARGE  OISK  *  ',P5.2,'T') 

WRITE  (6,1426)  RIMPON 

1426  FORMAT  (10X, 'SPACE  OCCUPIED  BY  IMPURITIES  *  ',P5.2,'X') 

WRITE  (6,0427)  TOTOEN 

0427  FORMAT  (IOX.'THc  TOTAL  PACKING  DENSITY  *  ',P5.2,'X') 

PRINT*,'  • 

WRITE  (6,0423)  LOCCEN 

0423  FORMAT  (lOX.'THE  LOCAL  PACKING  DENSITY  »  • , P6.2 , ' *'//) 

WRITE  (6,0404) 

0404  FORMAT  ( 10X ,  • OENSITY  CALCULATIONS  ARE  COMPLETE'//) 

PRINT*,'  MIT  RETURN’ 

REAO  C5.04C6) 

GOTO  0400 


44*«44«44*  CALCULATION  OF  ANGLE  ********** 

this  part  of  the  subroutine  is  used  to  calculate  the  angle  of  the  film 

GROWTH.  THIS  IS  ACCOMPLISHED  BY  SCANNING  THE  DEPOSITED  FILM  AT  ANGLES 
RANGING  FROM  1  TO  30  DEGREES.  AFTER  THIS  THE  LARGEST  "CORRELATION"  FAC 
TOR  IS  CHOSEN  AS  THE  ANGLE  OF  FILM  GROWTH.  THE  IDEA  HERE  IS  THAT  A 
LARGE  "CORRELATION"  WILL  OCCURR  WHEN  THE  VARYING  ANGLE  EQUALS  THE  ANGLE 
OF  FILM  GROWTH.  IT  HAS  3ESN  FCUNO  OUT  THROUGH  TRIAL  ANO  ERROR  THAT  THE 
BEST  RESULTS  ARE  03TAINE0  WHEN  THE  SEPARATION  BETWEEN  THE  TWO  HEIGHT  LE 
VELS  ARE  BETWEEN  10  ANC  20. 

0430  WRITE  (6,0431) 

0431  FORMAT  (lOX.'lN  THIS  PART  OF  THE  SU9R0UTIM6  THE  ANGLE  OF  THE  FILM 
CGROwTH'/,10X,'WILL  BE  CALCULATED  USING  A  SCANNING  TECHNIQUE .'/, 1  OX 
C.'IS  THIS  YOUR  CESIRS?'//, 11X,  • IF  YES  -  ENTER  1 • / , 1 1 X , • QThE RW I S E  - 
C  HIT  RETURN'//) 

C 

REAO  (5,0406)  IYAR 
IF  ( I V A R  .SO.  1)  THEN 
ANGCK  »  1 
GOTO  0459 
ELSE 

GOTO  0400 
ENCIF 


0459  WRITE  (6,0460) 

0460  FORMAT  CIOX.'IN  QROER  TO  SPEED  UP  THE  ANGLE  CALCULATIONS  TWO  HEIGH 
CT  LEVELS  NEED'/,I0X, 'TO  BE  SPECIFIED.  PLEASE  ENTER  THE  Y  VALUES  I 
CN  INTEGER  FORMAT.'//) 

C 

GOTO  0434 
C 

c  beginning  of  angle  analysis 
c 

0435  MAXCNT  *  299+T0PHGT-30THGT 
00  0437  C  «  0,89,1 
CORRCC)  «  0.0 
0437  CONTINUE 
C 

C  BEGINNING  OF  1N0IVICUAL  ANGLE  CHECKS 
C 

00  0439  A  «  1,80,1 
ANG  a  FLOAT (A  ) 

WRITE  (6,0455)  A 

0455  FORMAT  (23X ,'  ANALYZING  FILM  ANGLE  AT  ',12,'  DEGREES') 

RANG  =  ANG/<57. 2957751) 

TANANG  a  TANC  RANG ) 

CRSTRM  a  o.O 
00  0439  SUB  >  1,300,1 
C 

C  INTIALI2ATI0N  OF  ANGLE  ANALYSIS  VARIABLES 
C 

SDKCTR  a  0 
LOKCTR  a  o 
IOKCTR  a  0 
CSLCTR  a  0 
SLNOEN  a  0.0 
LLNOEN  a  0.0 
ILNOEN  a  0.0 
TLNOEN  a  0.0 

XPOS  a  CTANANG  *  TCPHGT)  ♦  SUB 
OXPOS  a  XPOS 
C 

C  START  OF  THE  ARRAY  SEARCH 
C 

00  0440  Y  a  TO PHG T ,1,-1 

XPOS  a  (TANANG  *  Y)  ♦  SUB 
ISHFT  a  (CAINTCXP0S/301))  *  300) 

NEWX  a  AINTCXPOS  -  ISHFT) 

OLOX  a  AINTCOXPOS  -  I S fl F T ) 

00  0441  X  a  OLOX, 1,-1 
C 

C  COUNTING  OF  THE  OISKS 
C 

IF  ( F ILM( X , Y )  .NE.  0)  THEN 

IF. (FILMCX.Y)  .EC.  1)  THEN 
SCKCTR  a  SOKCTR  ♦  1 
ELSE  IF  (FILMCX.Y)  .£0.  2)  TrtEN 
LOXCTR  a  LCKCTR  ♦  1 
ELSE  IF  (FILMCX.Y)  .ED.  8)  TrtEN 
IOKCTR  a  IOKCTR  ♦  1 
ENOIF 
ENOIF 


C 


X2  a  X  ♦  1 


I 
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C 

IF  CX2  .GT.  300)  THEN 

X2  •  X2  -  300 

ENOIP 

IF  ( PILMC  X2  »  Y )  .NE.  0)  THEN 

IF  ( PILM(  X2  *  T )  .EQ.  1)  THEN 

SOKCTR  a  SOKCTR  ♦  1 
ELSE  IF  (FILHCX2.Y)  . EQ. 

2)  THEN 

LOKCTR  a  LOKCTR  ♦  1 

ELSE  IF  CFILMCX2.Y)  .50. 

8)  THEN 

c 

IOKCTR  a  ICKCTR  ♦  1 
ENOIF 

ENOIF 

CELCTR  a  CELCTR  ♦  2 

IF  (CELCTR  .Ge.  MAXCNT)  GOTO 

0442 

0441 

IF  (X  .EQ.  NEWX  )  GOTO  0443 
CONTINUE 

0443 

OXPOS  a  XPOS 

0440 

IF  ( Y  .LE.  30THGT )  GOTO  0442 
CONTINUE 

C 

c  calculation  of  the  line  oensities 
c 

0442  SLNOEN  «  C C SCKCTR*1 . 570796 3 >/C ELCTR > 

LLNOEN  »  (CLCKCTR*6. 2831853) /CELCTR) 

ILNOEN  «  CCIDKCTR*100.53096)/CELCTR) 

TLNOEN  a  SLNCEN  ♦  LLNOEN  ♦  ILNOEN 
C 

IF  CTLNOEN  .LS.  <LOC06N/100))  then 

CCRFAC  *  C ( (LCC06N/1Q0)  -  TLNCEN)/( LQCOEN/lOO) )**2 
ELSE 

CCRFAC  «  C (TLNOEN  -  CLCCOEN/1 0 0 ) ) / C 1-C LOCC EN/1 00 ) ) )**2 
6N0IF 
C 

CRSTRM  »  CRSTRH  ♦  CORFAC 
0439  CONTINUE 
C 

C  DETERMINING  THE  CORRELATION  NUMBER 

c 


CORR(A)  =  CRSTRM/ 300 
0438  CONTINUE 
C 

C  FINDING  The  THREE  LARGEST  CORRELATION  NUMBERS 
C 

ANGLE1  *  0 

ANGLS2  «  0 

ANGLE  3  »  0 

LARGEl  «  0.0 

LARGE2  «  0.0 

LARGE  3  «  0.0 

00  0470  TEST  «  1,30,1 

IF  CCORRC  TEST )  .GT.  LARGEl)  THEN 
ANGLE 3  *  angle 2 
ANGLE  2  *  ANGLS1 
ANGLE1  a  TEST 
LARGE3  «  LARGE  2 
LARGE2  «  LARGEl 
LARGEl  «  CORR(TEST) 

ELSE  IF  (CGRRCTEST)  .GT.  LARGS2)  THEN 
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ANGLE3  >  ANGLE  2 
ANGLS2  »  TEST 
LARGE3  «  LARGE2 
LARGE2  *  CQRR(TEST) 

GOTO  0470 

ELSE  IF  (CORR(TEST)  .GT.  LARGE3)  THEN 
ANGLE3  »  TEST 
LARGE3  «  CORR(TEST) 

ENDIr 

0470  CONTINUE 
C 

C  CALCULATION  OF  THE  AVERAGE  OF  THE  TCP  THREE  VALUES 
C 

AVGANG  «  FLOAT (( ANGLE1  ♦  ANGLE2  ♦  ANGLE3)/3) 

C 

C  OUTPUT  OF  THE  ANGLE  ANALYSIS 
C 

PRINT#,'  • 

PRINT#,'  • 

WRITE  C6.04S1)  ANGLE!,  ANGLE2 , ANGLS3 

0451  FORMAT  C13X,'THE  TOP  THREE  CORRELATED  ANGLES  IN  OECSNOING  0R0ER  AR 
CE'/,T30,I2,aX,I2,8X,I2/> 

WRITE  (6,0453)  AVGANG 

0453  FORMAT  C 1 8X , 'THE  AVERAGE  CORRELATED  ANGLE  IS  ',F4.1,'  DEGREES'//) 
WRITE  (6,0452) 

0452  FORMAT  (20X,'THE  FILM  ANGLE  CALCULATIONS  ARE  COMPLETE'//) 

PRINT#,'  HIT  RETURN* 

RE AO  (5,0406) 

GOTO  0400 
ENO 


ft###**########*####**############ 

#  * 

*  MOVING  SUBROUTINE  # 

*  * 

«***#9#«99# ft***##*##*#*****#**### 

This  SUBROUTINE  IS  0IVI0E0  INTO  3  SECTIONS.  THE  FIRST  2  SECTIONS 
MOVE  the  oepositeo  film  to  a  file  for  storage  ano  back  again,  the 

LAST  SECTION  HELPS  IN  PREPARING  THE  DEPOSITED  FILM  FOR  PLOTTING. 
EACH  SECTION  CONTAINS  A  BRIEF  DESCRIPTION  OF  ITS  PURPOSE  ANO/OR 
ITS  APPROACH. 

SUBROUTINE  MOVE  (FILM,XCOORC,TCOORD) 

LOCAL  VARIABLES 

real  XC00RC(3OO,200),YCQORO(300,20O),XCOQR,YCO0R 
INTEGER  FILM(  300,2  00)  ,  IV AR,FIL,CTR,RECCTR 


#####«####  MOVING  interface  loop  ########## 

0300  WRITE  (6,0201) 

0  301  FORMAT  ( 18(2X/), 27X,  'MOVING  INTERFACE  LOOP ' /// ,  2 0 X , • E NT £R  APPROPRI 
CATS  NUMBER  TO  CON T  INUE  '  /  ,  2( 2 X/ ) , 2 5 X , • 1  -  STORE  FILM  M ATR I CES ' / , 2 SX 
C , ' 2  -  RECALL  FILM  M A TR IC S S ' / ,  2 S X ,  '  3  -  STORE  FILM  FCR  PLOTTING'/, 25 
C X  , ' A  -  RECALL  PLOTTING  V ALUE S ' / , 2 S X , • 9  -  EXIT  TO  MAIN  INTERFACE  LO 
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C0P'/,2SX,'CTRL/C  -  EXIT  TO  OPERATING  S Y STEP • / , 5 ( 2 X/ ) ) 

0305  READ  <5,0306)  IVAR 
0306  P0RMAT  CI1) 


IF  (IVAR 

.EQ. 

1)  then 

GOTO 

0310 

ELSE  IF 

(IVAR 

.EQ. 

2) 

THEN 

GOTO 

0330 

ELSE  IF 

(IVAR 

.EQ. 

3) 

THEN 

GOTO 

0350 

ELSE  IF 

(IVAR 

.EQ. 

4) 

THEN 

GOTO 

0380 

ELSE  IF 

(IVAR 

.EQ. 

8) 

THEN 

RETURN 

ENOIF 

WRITE  (6,0308) 

0308  FQRMAT(21X« ' YOU  MESSEO  UPH  -  TRY  AGAIN'///) 
GOTO  0305 


****«*a«3*  STORING  OP  FILM  MATRICES  ********** 

THIS  PART  OF  THE  SUBROUTINE  STORES  THE  VALUE  OF  FILM,  XCOORD,  YCOORD 
INTO  THE  OATA  FILE  F0R012.0AT.  THIS  IS  DONS  SO  THAT  FINISHED  FILMS 
CAN  BE  STCRSO  FOR  FUTURE  REFERENCE  AND  RECALLED  FOR  ANALYZING  IF  DE- 
SIREO.  THIS  METHOO  OF  STORAGE  IS  SCMEWH AT  INEFFICIENT  IN  THAT  UN¬ 
OCCUPIED  STORAGE  CELLS  ARE  ALSO  READ  INTO  THE  FILE.  THIS  METHOO  OF 
STORAGE  SHOULD  ONLY  BE  USEC  WHEN  THE  OPERATOR  HAS  ALOT  OF  SPACE  IN 
HIS  DIRECTORY. 

0310  WRITE  (6,0311) 

0311  FORMAT  (lOX.'IN  THIS  PART  OF  THE  SUBROUTINE  THE  CURRENT  VALUES  OF 
CTHfi  FILM'/, 10X, 'MATRICES  WILL  BE  STORED  IN  A  FILE  CALLED  FOR012.OA 
CT  FOR  FUTURE'/, 10X, 'REFERENCE. '/,10X, 'IS  THIS  YOU  CES IRE? ' // , 11X  ,  • 
CIF  YES  -  ENTER  1'/,11X, 'OTHERWISE  -  HIT  RETURN'//) 

C 

READ  (5,0306)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
PRINT*,'  ' 

PRINT*, 'WORKING' 

GOTO  031S 
ELSE 

GOTO  0300 
SNCIF 
C 

0315  OPEN  (UNIT  «  12,  FILS  *  'F0R012.DAT',  STATUS  *  'NEW') 

C 

00  0320  Y  >  1,200,1 
DO  0325  X  >  1,300,1 

WRITE  (  12,0326)  F I LM< X  ,  Y )  ,  XCOO RO ( X , Y ) , Y COQRO ( X , Y ) 

0326  FORMAT  (1 X , 1 2 , 3 X , F 3 . 4 , 3X  ,  F 8  .  A ) 

0325  CONTINUE 
0320  CONTINUE 
C 

CLOSE  (UNIT  »  12) 

C 

WRITS  (6,032?) 

0327  FORMAT ( 1  OX , ' THE  VALUES  OF  THE  FILM  MATRICES  HAVE  NCW  BEEN  COPIED'/ 
c*iox, ' into  the  fils  foroi2.dat. • //) 


oooooooooo 


PRINT*  » ' 

REAO  (5,0306) 
GOTO  0300 


MIT  RETURN* 


**********  RECALLING  OP  FILM  MATRICES  ********** 

THIS  PART  OF  THE  SUBROUTINE  RECALLS  THE  VALUES  OF  THE  FILE  F0R013.DAT 
SO  THAT  THEY  MAY  BE  ANALYZEO  ANO/OR  READIED  FOR  PLOTTING.  THE  VALUES 
REAO  IN  ARE  PLACED  IN  THE  MATRICS  FILM,  XCQQRQ,  YCOORO.  IN  DOING  SO 
THE  VALUES  STORED  BY  THE  PREVIOUS  SECTION  ARE  REAO  IN,  THEREFORE  THE 
VALUES  NEED  TO  OF  BEEN  COPIED  INTO  A  FILE  CALLEO  F0R013.DAT. 

0330  WRITE  (6,0331) 

0331  FORMAT  (10X,'JUST  A  REMINDER.  UPON  ENTERING  INTO  THIS  PART  OF  THE 
CSUBROUTINE'/.lOX, ‘THE  CURRENT  VALUES  OF  THE  FILM  MATRICES  WILL  BE 
CLOST.  THE  REA-'/, 10X, 'SON  FOR  THIS  IS  THAT  THE  NEW  VALUES  FROM  TH 
CE  FILE  FOR013.0AT*/,10X, 'WILL  BE  REAO  IN . • / , 10X ,  • DC  YOU  WISH  TO  SA 
CVE  THE  CURRENT  VALLES  OF  THE  FILM  MATRIC ES? * // ,  1 IX  ,  *  IF  YES  -  ENTER 
C  1*/,11X, 'OTHERWISE  -  HIT  RETURN*//) 

C 

REAO  (5,0306)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 
GOTO  0310 
ENOIF 

PRINT*, 'WORKING* 

C 

OPEN  (UNIT  »  13, FILE  *  • FORO 13. OAT* , STATUS  •  'OLD') 

REWINO  (UNIT  «  13) 

C 

DO  0335  Y  *  1,200,1 
00  0336  X  *  1,300,1 

REAO  (12,0337)  FILM( X , Y) , XC00RD( X , Y> , YCOORD( X , Y) 

0337  FORMAT  (1X,I2,3X,F8.*,3X,F8.*) 

0336  CONTINUE 
0335  CONTINUE 
C 

CLOSE  (UNIT  »  13) 

C 

WRITE  (6,03*0) 

03*0  FORMAT  (lOX.'THE  NEW  VALUES  OF  THE  FILM  MATRICES  HAVE  NOW  BEEN*/,1 
COX, 'REAO  IN  FORM  THE  FILE  FORO 1 3 . 0 AT . ' // ) 

PRINT*,'  HIT  RETURN* 

REAO  (5,0306) 

GOTO  0300 
C 
C 

c  **********  storing  of  film  for  PLOTTING  ********** 

c 

c  this  part  of  the  subroutine  stores  the  non-ierc  VALUES  of  The  MATRICES 
C  FILM,  XCOCRO,  ANO  YC0CR0  INTO  THE  FILES  FQR007.0AT  ANO  PCR014.0AT  FOR 
C  FUTURE  PLOTTING.  THIS'  IS  ACCOMPLISHED  BY  CHECKING  EACH  AOORESS  OP  THE 
C  TWO  OIMENSIONAL  AREAS  FOR  OCCUPATION  ANO  THEN  TRANSFSRING  ONLY  THOSE 
C  STORAGE  CELLS  THAT  HAVE  NONE  ZERO  INFORMATION. 

C 

0350  WRITE  (6,0351) 

0351  F0RMAT(10X,*WARNINSI  UPON  ENTERING  INTO  THIS  PART  OF  TmE  SUSROUTI 
CNS  THE'/, 10X, 'CURRENT  VALUES  OF  THE  FILM  MATRICES  WILL  BE  COPIEO  I 
CNTG  THE'/, 10X, 'FILES  P0R007.CAT  ANC  FOROi*.OAT.  IN  OOlNG  SO,  TME 
CAOCRESS  OP'/.lOX, 'The  INDIVIDUAL  POINTS  WILL  BE  LOST  SO  THAT  NO  AN 


oooooonooooo 
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CALYSIS  CAN  3S'/,10X,  ‘EVALUATED  OIRECTLY.  */,10X  ,  •  IS  THIS  YOUR  OESIR 
CE? ' // » 11X , • IP  YES  -  ENTER  1 •/, 11 X ,• OTHERWISE  -  HIT  RETURN'//) 

C 

REAO  (5,0306)  I  VAR 
IF  (IVAR  .EQ.  1)  THEN 
PRINT*, 'WORKING* 

GOTO  0355 
ELSE 

GOTO  0300 
ENOIF 
C 

0355  OPEN  (UNIT  ■  7, PILE  «  • P0R007 . OAT' , STATUS  »  'NEW') 

OPEN  (UNIT  »  14, PILE  ■  'F0R014. OAT', STATUS  «  'NEW') 

.  C 

DO  0357  Y  *  1,200,1 
00  0358  X  -  1,300,1 

IF  (FILM(X,Y)  .NE.  0)  THEN 

WRITE  (7,0359)  FILM(X,Y),XCOORO(X,Y),YCOORO(X,Y) 

0359  FORMAT  ( IX , 12 , 3X , F 8 . 4, 3X  , F8 . 4 ) 

IF  ( PILM( X , Y )  .EQ.  2)  THEN 

WRITE  (14,0356)  FILM( X , Y ) , XCQORO( X , Y) , YCOQRO CX , Y ) 

0356  FORMAT  ( IX  ,  1 2 , 3 X , F8 . 4, 3X , PB . 4 ) 

ENOIF 

ELSE 

GOTO  0358 
ENOIF 

0358  CONTINUE 
0357  CONTINUE 
C 

CLOSE  (UNIT  »  14) 

CLOSE  (UNIT  »  7) 

C 

WRITE  (6,0360) 

0360  FORMAT  ( 1  OX , ' THE  VALUES  OF  THE  FILM  MATRICES  HAVE  BEEN  COPIED  TO'/ 
C  « 10X , ' THE  FILES  CALLEO  F0R007.0AT  AND  FORO 14 . DAT .'//, 1  OX ,' PLOTTING 
C  OF  THE  FILM  MATRICES  CAN  NOW  OCCUR!'//) 

PRINT*,'  HIT  RETURN* 

READ  (5,0306) 

GOTO  0300 


*****«***«  RECALLING  OF  FILM  FOR  ANALYZING  ********** 

this  last  part  OF  The  subroutine  takes  The  values  film,  xcooro,  ycooro 

from  THE  file  FQR008.0AT  50  THAT  THEY  CAN  BE  ANALYZED.  THIS  IS  ACCOM¬ 
PLISHED  3t  LOOKING  AT  THE  VALUES  OF  XCOORO  ANO  YCOORO  ANO  TRUNCATING  THE 
VALUES.  AFTER  TRUNCATION  ThE  FILM  MATRIX  IS  THEN  GIVEN  THE  PROPER  VALUES. 
THIS  SECTION  IS  THE  COUNTERPART  OF  THE  PREVIOUS  SECTION  IMPLYING  THAT  THE 
VALUES  MUST  8E  COPIED  INTO  A  FILE  F0R003.0AT  BEFORE  ATTEMPTING  TO  READ 
FROM  THE  FILS. 

0380  WRITE  (6,0331) 

0381  FORMAT  ( 1 0 X , • C AUT ICN I  UPON  ENTERING  INTO  THIS  SECTION  OF  THE  SU8R0 
CUTINE,  THE'/, 10X, 'CURRENT  VALUES  OF  F0R0C8.DAT  WILL  BE  REAO  INTO  T 
CHE  FILM  MATRICES ' /,10x , ' DESTROY InG  ThE  CLO  V ALUE S . • / , 1 0 X , ' H A VE  THE 
C  FILM  MATRICES  BEEN  INITIALIZED  to  ZERO?'//, UX,  •  if  yes  -  ENTER  1* 
C/.UX,  'OTHERWISE  -  HIT  RETURN'//) 

C 

REAO  (5.03C6)  IVAR 
IF  (IVAR  .EQ.  1)  THEN 


•*.*  va-  *v*\ ■"«' tj  w;  vj v j v. »v ipj w .  pj r>  /;  r* r_  i-  ™> t> ".■ 
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SOTO  0332 
ELSE 

SOTO  0386 
ENOIF 
C 

0386  WRITE  (6.0339) 

0389  FORMAT  (lOX.'OO  YOU  WISH  TO  INITIALIZE  THE  FILM  MATRICES? •//,  11X , ' 
CIr  YES  -  ENTER  1 •/, 11X, • OTHERWISE  -  HIT  RETURN*//) 

C 

REAO  (5*0306)  IV AR 
IF  (IVA'R  .EQ.  1)  THEN 
PRINT*** WORKING* 

GOTO  0397 
ELSE 

PRINT*,'  • 

GOTO  0382 
ENOIF 
C 

0397  00  0398  X  «  1,300,1 

00  0399  Y  *  1,200,1 
FILMCX.Y)  *  0 
0399  CONTINUE 
0398  CONTINUE 
C 

0382  WRITE  (6,0383) 

0383  FORMAT  (10X,*IN  OROER  TO  COMPLETE  THIS  PART  OF  THE  SUBROUTINE  ONE 
CMUST  KNOW*/, 10X, 'THE  NUMBER  OF  RECORDS  IN  THE  FILE  TO  BE  REAO.  IF 
C  YOU  KNOW  THE'/, 10X, 'CORRECT  NUMBER  ENTER  THAT  BELOW.  IF  NOT,  EXI 
CT  THE  PROGRAM  ANO ' / , 10X , ' FIND  OUT ! I  -  USE  CTRL/C*//) 

C 

PRINT*, 'ENTER  THE  NUMBER  OF  RECORDS' 

REAO  (5,0385)  NUMREC 
0335  FORMAT  (15) 

C 

PRINT*,*  ' 

PRINT*, 'WORKING* 

C 

OPEN  (UNIT  »  8, FILE  »  • FOR008. CAT* .STATUS  •  *OLD*) 

REWINO  (UNIT  *  8) 

C 

RECCTR  «  0 

00  0390  CTR  «  1, NUMREC, 1 

REAO  (8,0391)  FIL,  XCOOR,  YCOOR 
0391  FORMAT  (  1  X  ,  12 , 3 X  ,  F8 . 4 , 3 X , F8 . 4 ) 

X  *  AINT(XCOOR) 

Y  »  AINT(YCQOR) 

FILM(X.Y)  «  FIL 
RECCTR  «  RECCTR  ♦  1 
IF  (RECCTR  .60.  1000)  THEN 
WRITE  (6,0388)  CTR 

0338  FORMAT  C23X.I5,'  RECCROS  REAO  IN') 

RECCTR  ■  0 
ENOIF 
0390  CONTINUE 
C 

CLOSE  (UNIT  »  15) 

C 

PRINT*,'  • 

WRITS  (6,3396)  (CTR-1) 

0396  FORMAT  (21X,'A  TOTAL  OF  ',15,'  RECCROS  wSRE  READ  IN'//) 
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WRITE  (6*0395) 

0395  FORMAT  (lOX.'THE  NEW  VALUES  OF  THE  FILM  MATRICES  HAVE  SEEN  READ  IN 
C  FROM  THE*/, 10X, 'FILS  FGR008.0AT  AND  ARE  READY  TO  3E  ANALYZED.'//) 
PRINT*,'  HIT  RETURN* 

REAO  (5,0306) 


RETURN 


Program  Documentation 

The  purpose  of  this  appendix  is  to  provide  some  documentation  of 
the  program  listed  in  appendix  B.  This  is  accomplished  by  listing  the 
variables  used  in  the  program  and  giving  a  short  description  of  their 
purpose. 


Introduction 

This  program  is  the  second  in  a  series  of  programs  written  by  stu¬ 
dents  attending  the  Air  Force  Institute  of  Technology  located  at  Wright 
Patterson  AFB,  Ohio.  The  purpose  of  this  program,  like  its  predecessor, 
is  to  simulate  the  vapor  deposition  of  thin  films.  It  is  written  in  For¬ 
tran  77  and  contains  a  statement  using  IMSL  (International  Mathematical 
and  Statistical  Libraries),  a  library  of  intrinsic  functions  used  for 
statistical  purposes.  It  was  designed  to  run  on  a  VAX/VMS  computer  sys¬ 
tem. 

The  program  is  organized  into  five  main  sections  which  includes  the 
the  main  start-up  and  four  subroutines.  The  main  start-up  consist  of  a 
help  message  in  order  to  get  the  user  aquainted  with  the  program  and  an 
interface  loop.  The  substrate  subroutine  allows  the  user  to  intialize, 
create,  store,  recall  and  move  substrates.  The  deposition  subroutine  is 
used  to  initialize  the  main  arrays,  read  in  deposition  parameters  and  de¬ 
posit  the  film  onto  the  substrate.  The  analysis  subroutine,  provides 
the  user  the  ability  to  check  the  density  of  the  film  and  the  angle  of 
growth  without  having  to  print  out  the  film.  The  moving  subroutine 
allows  the  user  to  store  and  recall  films  for  plotting  and  analysis. 


Variable  Listing 

The  following  is  a  listing  of  the  variables  used  in  the  program 
listed  in  appendix  B  along  with  a  short  description  of  their  purpose.  A 
qualifier  on  the  end  of  the  description,  tells  in  what  subroutine  of  the 
program  the  variable  was  used.  The  listing  is  in  alphabetical  order. 


A  -  A  variable  derived  in  appendix  A  used  to  find  the  rest  point 
of  the  incident  disk.  -  depo 

ANDVAR  -  The  number  of  degrees  (full  angle)  that  the  angle  of  in¬ 
cidence  should  vary  about  the  main  angle  of  deposition.  -  depo 

ANG  -  The  main  angle  of  deposition  measured  from  the  subsrate  nor¬ 
mal.  -  depo, anal 

ANGCK  -  This  variable  is  used  to  help  in  the  angle  analysis  of  the 
film,  especially  in  the  density  check.  -  anal 

ANGINC  -  The  amount  of  angle  to  be  incremented  or  decremented  from 
the  old  angle  of  deposition.  -  depo 

ANGLE1,2,3  -  Three  variables  used  to  hold  the  angles  corresponding 
to  the  largest  correlation  numbers.  -  anal 

ANGMAX  -  A  variable  which  contains  the  largest  angle,  the  angle  of 
incidence  can  have  when  angle  variation  takes  place.  -  depo 

ANGMIN  -  A  variable  which  contains  the  smallest  angle,  the  angle  of 
incidence  can  have  when  angle  variation  takes  place.  -  depo 

AVGANG  -  The  average  of  the  angles  with  the  three  largest  correla¬ 
tion  factors  used  in  finding  the  angle  of  growth  by  the  film.  -  anal 

A2  -  The  x  intercept  of  the  trajectory  for  the  center  of  the  inci¬ 
dent  disk.  -  depo 

B  -  A  variable  derived  in  appendix  A  used  to  find  the  rest  point  of 
an  incident  disk.  -  depo 

BNCCTR  -  A  counting  variable  used  to  monitor  if  an  extra  mobility 
is  to  occur.  -  depo 

BNCVAR  -  The  number  specifying  (indirectly)  the  percentage  of  disks 
to  undergo  an  extra  mobility.  -  depo 

BOTHGT  -  Specifies  the  bottom  height  to  be  checked  in  the  density 
and  angle  analysis  routines.  -  anal 


BSHFT  -  This  variable  used  in  the  mobility  routine  helps  maintain 
the  periodic  nature  of  the  film  and  is  like  the  other  shift  variables,- 
depo 

BTOL  -  This  variable  is  just  like  TOL  but  is  used  in  the  mobility 
routine.  -  depo 

C  -  A  variable  derived  in  appendix  A  used  to  find  the  rest  point  of 
an  incident  disk.  -  depo 

CELCTR  -  Counts  the  number  of  cells  during  a  streamer  check  in  the 
the  angle  analysis  routine.  -  anal 

COLL  -  A  variable  used  to  keep  track  of  whether  or  not  a  collision 
has  occured  during  checking  of  the  collision  corridor  by  the  streamers. 

-  depo 

CORFAC  -  The  correlation  number  of  an  individual  streamer  used  to 
find  the  angle  of  growth  by  the  film.  -  anal 

CORR  -  An  array  used  to  hold  the  correlation  numbers  for  finding 
the  angle  of  growth  by  the  film.  -  anal 

COSANG  -  The  cosine  of  RANG.  -  depo 

CRSTRM  -  The  summation  of  CORFAC  across  the  bottom  of  the  sub¬ 
strate.  -  anal 

CTR  -  The  variable  used  to  specify  the  number  of  streamers  in  the 
scanning  routine.  -  depo 

D  -  A  variable  defined  from  the  variables  A,  K  and  Rl.  -  depo 

DATE  -  The  date  and  time  entered  are  used  to  modify  the  seed  of  the 
random  number  generator  in  the  deposition  subroutine.  -  depo 

DENCK  -  Used  to  signal  that  the  density  algorithm  is  to  be  used.  - 

anal 

DEPCTR  -  The  main  deposition  counter  which  keeps  track  of  how  many 
disks  are  deposited.  -  depo 

DRESTX  -  This  default  variable  stores  the  value  of  RESTX  and  is 
used  in  case  a  test  is  failed  in  the  mobility  routine.  -  depo 

DRESTY  -  The  counter  part  of  DRESTX  but  stores  the  value  of  RESTY.- 

depo 

DSEED  -  The  seed  numbers  used  to  determine  the  random  numbers  in 
deposition.  -  depo 

F  -  A  variable  used  to  find  the  y  coordinate  of  the  collision  point 
and  is  derived  in  appendix  A.  -  depo. 
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FIL  -  A  counterpart  to  FILM  used  to  read  back  in  the  film  for  anal¬ 
ysis.  -  move 


FILM  -  A  300  by  200  array  used  to  keep  track  of  whether  a  cell  is 
occupied  by  a  disk.  If  its  value  is  zero  the  cell  is  unoccupied.  -  sub, 
depo ,anal , move 

G  -  A  variable  used  to  find  the  y  coordinate  of  the  collision  point 
and  is  derived  in  appendix  A.  -  depo 

H  -  A  variable  defined  in  appendix  A  and  is  used  to  find  the  rest 
point  of  the  incident  disk.  -  depo 


HCTR  -  This  counting  variable  keeps  track  of  how  many  overlapping 
disks  occur  during  the  deposition  process.  -  depo 

HDSIZE  -  This  variable  holds  the  last  used  value  of  SIZE  when  an 
impurity  is  deposited.  -  depo 

HEIGHT  -  A  variable  used  to  store  the  largest  height  of  the  sub¬ 
strate.  -  sub, depo 

HGT  -  A  variable  read  in  during  the  creation  of  the  substrate.  Its 
value  gives  the  height  or  the  number  of  disks  to  be  placed  in  a  column.  - 
sub 

HGTIMP  -  A  number  determined  randomly.  If  the  value  of  YHGT  equals 
or  is  greater  than  this  value  an  impurity  is  deposited  if  this  is  de¬ 
sired.  -  depo 

HGTLAY  -  Its  value  specifies  the  height  of  each  layer  in  terms  of 
the  spatial  arrays.  -  depo 

IDKCTR  -  Counts  the  number  of  impurities  in  a  streamer  search  to 
find  the  angle  of  growth  by  the  film.  -  anal 

IDNCTR  -  Counts  the  number  of  impurities  in  the  local  density  var¬ 
iation  check.  -  anal 

ILNDEN  -  The  line  density  of  streamer  for  impurities  in  checking 
for  angle  of  growth  by  the  film.  -  anal 

IMPCTR  -  Counts  the  total  number  of  impurities  in  the  film.  -  anal 

IMPHGT  -  This  variable  stores  the  truncated  value  of  the  y  posi¬ 
tion  of  the  impurity  after  it  is  deposited.  -  depo 

IMPSET  -  If  equal  to  one  an  impurity  is  about  to  be  or  has  been 
deposited  into  the  film.  Its  value  may  change  to  zero  if  the  impurity 
is  covered.  -  depo 
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IMPUR  -  If  this  has  the  value  of  one  then  an  impurity  will  be  de¬ 
posited  randomly  into  the  mircrostructure.  -  depo 

IMPX  -  This  variable  stores  the  truncated  value  of  the  x  position 
of  the  impurity  after  it  is  deposited.  -  depo 

ISHFT  -  The  truncated  value  of  XPOS  used  to  help  maintain  periodic 
boundaries.  -  depo 

IVAR  -  This  interface  variable  used  throughout  the  program  is  re¬ 
sponsible  for  reading  in  answers  to  questions  asked  by  the  program.  - 
sub , depo , anal , move 

K  -  A  variable  defined  in  appendix  A  and  is  used  to  find  the  rest 
point  of  the  incident  disk.  -  depo 

LARGE1,2,3  -  Three  values  which  store  the  largest  correlation  num¬ 
bers  for  finding  the  angle  of  growth  by  the  film.  -  anal 

LAYER  -  If  this  has  the  value  of  one  then  a  multiple  layered  film 
will  be  produced.  -  depo 

LDNCTR  -  Counts  the  number  of  large  disks  in  the  local  density  var¬ 
iation  check.  -  anal 

LDKCTR  -  A  counter  of  large  disks  in  the  streamer  search  for  find¬ 
ing  the  angle  of  growth  by  the  film.  -  anal 

LEFT  -  The  default  boundary  condition  on  the  left  side  for  vertical 
variations.  -  anal 

LIMPDN  -  The  impurity  packing  density  in  the  local  area  check  - 


LLNDEN  -  The  line  density  of  the  large  disks  during  a  streamer 
search  to  find  the  angle  of  growth  by  the  film.  -  anal 

LLRGDN  -  The  large  packing  density  in  the  local  area  check.  -  anal 

LOCCTR  -  Counts  the  number  of  cells  in  the  local  density  variation 
check.  -  anal 


LOCDEN  -  The  total  local  packing  density.  -  anal 

LRGCTR  -  Counts  the  total  number  of  large  disks  in  the  film.  -  anal 

LSMLDN  -  The  small  packing  density  in  the  local  area  check.  -  anal 

MAX  -  The  value  of  this  variable  determines  how  far  to  the  right 
and  up,  the  area  search  is  to  take  place  around  the  collision  disk 
looking  for  possible  rest  disk.  -  depo 
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MIN  -  The  value  of  this  variable  determines  how  far  to  the  left 
and  down,  the  area  search  is  to  take  place  around  the  collision  disk 
looking  for  possible  rest  disk.  -  depo 

ND  -  A  counting  variable  used  to  keep  track  of  the  number  of  disks 
that  have  been  deposited.  -  depo 

NEWX  -  The  truncated  value  of  XPOS  minus  ISHFT  and  is  used  to  de¬ 
termine  to  what  address  the  search  is  to  take  place  along  an  individ- 
al  row  during  the  scanning  routine.  -  depo 

NODSK  -  The  number  of  disk  to  be  deposited  in  the  deposition  of  the 
film  if  the  spatial  arrays  do  not  become  full.  -  depo 

NOMOB  -  A  variable  used  in  determining  the  mobility  of  the  disk. 

If  equal  to  one  no  mobility  will  occur.  -  depo 

NUMREC  -  The  number  of  records  to  be  read  back  in  for  analysis.  - 

move 

NXINT  -  The  value  of  XINT  minus  ISHFT  and  is  used  to  account  for 
the  perodic  nature  of  the  film.  -  depo 

OLDANG  -  The  main  angle  of  deposition  measured  from  the  substrate 
normal.  -  depo 

OLDX  -  The  truncated  value  of  OXPOS  minus  ISHFT  and  is  used  to  de¬ 
termine  from  what  address  the  search  began  along  an  individual  row  dur¬ 
ing  the  scanning  routine.  -  depo 

OXINT  -  A  variable  used  to  minimize  the  calculations  made.  It  con¬ 
tains  the  truncated  x  position  of  the  last  streamer  intercept.  -  depo 

OXPOS  -  This  variable  contains  the  last  new  value  of  XPOS  for  an 
individual  row  during  the  scanning  by  a  streamer.  -  depo 

OYINT  -  A  variable  used  to  minimize  the  calculations  made.  It  con¬ 
tains  the  truncated  y  position  of  the  last  streamer  intercept.  -  depo 

P  -  This  variable  has  the  value  of  1  or  -1  and  is  used  in  incre¬ 
menting  or  decrementing  the  angle  of  incidence.  -  depo 

R  -  This  variable  gives  the  distance  squared  that  a  disk  has  to 
move  before  coming  to  rest.  -  depo 

RAD2  -  The  radius  of  the  collision  disk  contact  circle  squared  and 
is  used  to  find  the  y  coordinate  of  the  collision  point.  -  depo 

RANDX  -  The  random  number  between  I  and  301  used  to  deposit  the 
disks  evenly  across  the  substrate.  Its  value  is  the  x  intercept  of 
the  left  most  streamer.  -  depo 


RANG  -  The  angle  of  deposition  in  radians.  -  depo 


RANDY  -  The  random  number  between  1  and  180  produced  in  order  to  de¬ 
cide  where  the  impurity  is  to  be  deposited  in  the  mircrostructure.  - 
depo 


RANNUM  -  The  randum  number  between  0  and  1  produced  by  the  IMSL 
subroutine  named  ggubfs.  -  depo 

RBNC  -  An  array  seven  in  length  and  is  used  like  SHORT  to  keep 
track  of  how  far  a  disk  has  to  move  before  coming  to  rest.  It  is  used 
in  the  mobility  routine.  -  depo 

RECCTR  -  A  counting  variable  used  to  keep  track  of  how  many  records 
are  read  back  in  for  analysis.  -  move 

RESTX  -  The  x  coordinate  of  the  final  position  of  an  incident  disk. 

-  depo 

RESTY  -  The  y  coordinate  of  the  final  position  of  an  incident  disk. 

-  depo 

RIGHT  -  The  default  boundary  on  the  right  side  used  for  vertical 
variation  checks.  -  anal 


RIMPDN  -  The  total  packing  density  due  to  impurities.  -  anal 

RINC  -  The  incremental  radius  used  to  calculate  the  values  of  R1 

and  R2  and  is  determined  by  the  size  of  the  incident  disk.  -  depo 

RLRGDN  -  The  total  packing  density  due  to  large  disks.  -  anal 

RSMLDN  -  The  total  packing  density  due  to  small  disks.  -  anal 

R1  -  The  radius  of  the  collision  disk  contact  circle  and  is  used  in 
finding  the  rest  point  of  the  incident  disk.  -  depo 

R2  -  The  radius  of  a  possible  rest  disk  contact  circle  and  is  used 
in  finding  the  rest  point  of  the  incident  disk.  -  depo 

SC2ANG  -  The  secant  squared  of  RANG.  -  depo 


SDKCTR  -  Counts  the  number  of  small  disks  in  the  streamer  search 

for  finding  the  angle  of  growth  by  the  film.  -  anal 

SDNCTR  -  Counts  the  number  of  small  disks  in  the  local  density  var¬ 

iation  check.  -  anal 


SET  -  If  equal  to  one  an  impurity  is  about  to  be  or  has  been  de¬ 
posited.  -  depo 


SHFT  -  A  variable  like  ISHFT  used  to  keep  the  film  periodic  but  is 
used  in  the  rest  point  determination.  -  depo 

SHORT  -  A  variable  used  to  keep  track  of  how  far  an  incident  disk 
has  to  move  before  coming  to  rest.  The  smaller  the  value  the  shorter 
the  distance.  -  depo 

SINCT  -  The  x  intercept  of  the  current  streamer.  -  depo 

SIZE  -  This  variable  contains  information  pertaining  to  the  dia¬ 
meter  of  the  disk  in  question.  If  equal  to  one  the  diameter  equals  the 
square  root  of  two,  if  equal  to  two  the  diameter  is  two  times  the  square 
root  of  two,  it  equal  to  eight  the  diameter  equals  eight  times  the 
square  root  of  two.  -  depo, anal, move 

SLNDEN  -  The  line  density  due  to  the  small  disks  and  is  used  in 
finding  the  angle  of  growth  by  the  film.  -  anal 

SMLCTR  -  Counts  the  total  number  of  small  disks  in  the  film.  -  anal 

STRCTR  -  A  variable  which  contains  the  current  value  of  incrementa¬ 
tion  of  the  number  of  streamers  used  to  search  the  collision  corridor, 
depo 

SUBBUF  -  A  300  by  10  array  used  to  store  the  occupancy  data  of  the 
substrate  in  its  development  before  being  transfered  to  FILM.  -  sub 

TANANG  -  The  tangent  of  RANG.  -  depo 

TESTY  -  A  variable  used  to  make  sure  the  disk  is  not  below  the  film 
array.  -  depo 

TEST1  -  A  variable  used  for  a  number  of  test.  This  is  used  to  make 
sure  a  collision  is  possible  between  two  disks.  -  depo 

TEST2  -  A  variable  like  TEST1.  This  is  used  to  make  sure  the  two 
contact  circles  do  indeed  intercept  one  another.  -  depo 

TEST3  -  A  variable  almost  exactly  like  TEST2  but  this  one  uses  the 
x  values  to  check  for  intersection  between  the  circles.  -  depo 

TEST4  -  A  variable  used  in  conjuction  with  TOL  to  make  sure  that 
the  values  of  XI  and  Y1  go  together.  -  depo 

TLNDEN  -  The  total  line  density  of  the  streamer  used  to  find  the 
angle  of  growth  by  the  film.  -  anal 

TOL  -  A  variable  used  in  conjunctin  with  TEST4  to  make  sure  that 
the  values  of  XI  and  Y1  go  together.  -  depo 

TOPHGT  -  This  is  used  to  specify  the  top  height  in  the  density  and 
angle  analysis  routines.  -  anal 


TOTCTR  -  Counts  the  total  number  of  disks  in  the  film.  -  anal 

TOTDEN  -  The  total  packing  density  of  the  film.  -  anal 

TPX3  -  The  truncated  value  of  RESTX,  it  is  used  to  check  for  prior 
occupancy  by  another  disk.  -  depo 

TPY3  -  The  truncated  value  of  RESTY,  it  is  used  to  check  for  prior 
occupancy  by  another  disk.  -  depo 

TRESTX  -  A  temporary  storage  place  used  to  store  the  last  value 
of  RESTX.  -  depo 

TRESTY  -  A  temporary  storage  place  used  to  store  the  last  value 
of  RESTY.  -  depo 

TSHORT  -  A  temporary  storage  place  used  to  store  the  last  value 
of  SHORT.  -  depo 

TXBNC  -  The  truncated  value  of  XBNC  and  is  used  to  check  for  oc¬ 
cupancy  in  the  mobility  routine.  -  depo 

TYBNC  -  The  truncated  value  of  YBNC  and  is  used  to  check  for  oc¬ 
cupancy  in  the  mobility  routine.  -  depo 

WIDTH  -  A  variable  used  to  determine  the  values  of  MIN  and  MAX.  - 

depo 

XBNC  -  An  array  seven  in  length  and  stores  possible  x  values  of 
rest  points  for  the  mobility  routine.  -  depo 

XBCK  -  This  is  used  in  the  area  search  of  the  mobility  routine.  - 

depo 

XC  -  A  300  by  10  array  used  to  store  the  y  coordinates  of  the  disks 
in  the  development  of  the  substrate  before  being  transfered  to  YC00RD.  - 
sub 

XCOL  -  Its  value  gives  the  current  column  that  the  substrate  sub¬ 
routine  is  working  on  in  the  creation  of  a  substrate.  -  sub 

XCOLD  -  A  variable  used  to  store  the  old  values  of  XC  in  the  manu¬ 
facturing  of  the  substrate.  -  sub 

XCOLL  -  The  integer  value  which  contains  the  x  address  of  the  col¬ 
lision  disk.  -  depo 

XC00R  -  A  counterpart  of  XC00RD  used  for  reading  back  in  films  for 
analysis.  -  move 

XC00RD  -  A  300  by  200  array  used  to  store  the  x  coordinates  of  the 
disks.  -  sub, depo, anal, move 


XDIST  -  A  variable  used  in  conjunction  with  YDIST  to  find  how  far  a 
disk  has  to  move  before  coming  to  rest.  -  depo 

XINCT  -  This  variable  contains  the  x  intercept  of  the  center  of  the 
incident  disks  trajectory.  -  depo 

XPOS  -  The  address  of  the  cell  used  to  start  the  search  of  the  col¬ 
lision  corridor.  -  depo 

XI  -  A  variable  used  to  calculate  the  two  possible  x  values  of  the 
rest  point.  -  depo 

X2  -  The  x  position  of  the  incident  disk  immediately  upon  collision 
with  the  collision  disk.  -  depo 

YBCK  -  A  variable  used  in  the  area  search  of  the  mobility  routine.  - 

depo 

YBNC  -  An  array  seven  in  length  and  stores  possible  y  values  of 
rest  points  for  the  mobility  rountine.  -  depo 

YC  -  A  300  by  10  array  used  to  store  the  y  coordinates  of  the  disks 
in  the  development  of  the  substrate  before  being  transferred  to  XCOORD.  - 
sub 

YCOLD  -  A  variable  used  to  store  the  old  values  of  YC  in  the  manu¬ 
facturing  of  a  substrate.  -  sub 

YCK  -  A  variable  used  to  recalculate  the  variable  YGHT  when  an  im¬ 
purity  is  deposited.  -  depo 

YCOLL  -  The  integer  value  which  contains  the  y  address  of  the  col¬ 
lision  disk.  -  depo 

YC00R  -  The  counterpart  of  YC00RD  used  to  read  back  in  the  film  for 
analysis.  -  move 

YCOORD  -  A  300  by  200  array  used  to  store  the  y  coordinates  of  the 
disks.  -  sub. depo, anal, move 

YDIST  -  A  variable  used  in  conjunction  with  XDIST  to  find  out  how 
far  a  disk  has  to  move  before  coming  to  rest.  -  depo 

YFILM  -  A  variable  used  to  store  the  largest  height  of  the  sub¬ 
strate  in  conjuntion  with  determing  the  thickness  of  the  multilayers.  - 
depo 

YHGT  -  This  variable  stores  the  y  starting  point  for  the  scanning 
routine  for  depositing  the  film.  -  depo 

YI  -  A  variable  used  to  calculate  the  new  possible  y  coordinate  of 
the  collision  point.  -  depo 


YTEST  -  A  variable  uaed  in  the  mobility  routine  to  make  sure  that 
the  film  does  not  become  less  dense.  -  depo 

Y1  -  A  variable  used  to  calculate  the  two  possible  y  values  of  the 
rest  point.  -  depo 

Y2  -  The  y  position  of  the  incident  disk  immediately  upon  collision 
with  the  collision  disk. 

Y3  -  The  truncated  value  of  RESTY  and  is  used  to  make  sure  that  mo¬ 
bility  is  possible  near  the  substrate.  -  depo 
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A  VAX-11/785  computer  was  used  to  simulate  the  two-dimensional  .>;v‘ 
growth  of  thin  films  produced  by  vapor  deposition.  In  this  model 
molecules  and  impurities  were  represented  by  three  different  sized 
disks.  In  order  to  simulate  varying  deposition  conditions  and  evapor- 
ants,  several  variable  parameters  were  introduced.  Among  these  para¬ 
meters  were  the  variation  of  the  deposition  angle  about  some  main  an¬ 
gle,  the  mobility  of  the  disks  upon  collision,  the  ability  to  intro¬ 
duce  impurities  into  the  mircrostructure,  the  simulation  of  multi¬ 
layered  coatings  and  the  ability  to  introduce  imperfections  into  the 
substrate. 

The  results  obtained  by  this  model  show  that  disks  can  be  used 
to  simulate  some  of  the  main  features  exhibited  by  vapor  deposited 
films.  Among  these  features  are  the  formation  of  columns  and  their 
compliance  with  the  "tangent  rule*^  and  the  disappearance  of  this 
structure  in  the  case  of  large  disk  mobility.  Another  feature  found 
to  be  exhibited  in  the  modeled  films  is  that  under  certian  conditions, 
impurities  and  substrate  imperfections  can  produce  large  voids  and/or 
nodules.  Other  characteristics  found  in  the  simulated  films  include 
pores  which  could  allow  water  absorption,  and  increased  packing  den¬ 
sity  for  films  produced  with  angle  variations  along  with  a  moderate 
amount  of  disk  mobility.  (  ■"  •  ■ 
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